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Chapter  I 
INTRODUCTION 

Based  on  numerous  studies  of  laser  system  applications  in  recent 
years,  important  applications  of  pulsed  laser  operation  of  all  types  of 
lasers  have  been  apparently  identified.  The  pulsed  operation  of  a high 
energy  gas  laser  presents  special  problems  in  the  area  of  shock  wave 
generation  due  to  the  sudden  energy  deposition  and  the  interaction  of 
these  waves  with  other  components  of  the  laser  system.  These  problems 
are  present  in  both  closed  and  open  cycle  laser  systems  and  in  all  gas 
lasers  such  as  electric  discharge,  chemical,  and  gas  dynamic  lasers  when 
operated  in  a pulsed  mode.  In  addition  to  shock  wave  generation,  other 
transient  gas  flow  thermodynamic  problems  are  involved  during  the  start 
up  of  a high  energy  gas  laser  system  from  a room  temperature,  zero  flow 
initial  condition  until  the  establishment  of  "steady"  flow  conditions 
corresponding  to  sustained  pulsed  operation.  These  start  up  transients 
are  characteristic  of  all  high  energy  gas  laser  systems,  both  pulsed  and 
CW  and  both  open  and  closed  cycle  to  varying  degrees.  The  purpose  of  the 
work  reported  here  is  to  investigate  some  of  the  major  transient  fluid 
and  thermal  characteristics  of  high  energy  pulsed  laser  operation  and  to 
study  methods  of  reducing  or  eliminating  the  problems  associated  with  such 
systems.  The  pulsed  closed  cycle  electric  discharge  laser  (EDL)  was  used 
as  the  basis  for  study  in  this  work  since  it  has  most  of  the  major  fluid 
and  thermal  operating  characteristics  associated  with  all  high  energy  gas 
laser  systems.  The  following  major  technical  issues  were  determined  to  be 
associated  with  the  development  of  the  pulsed  closed  cycle  high  energy 
laser  system: 


A 
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A.  Transients 

1.  Start  up  recirculator 

2.  Start  up  laser  input  power 

3.  Acoustic  waves/Attenuation 

B.  Boundary  Layer  Interaction 

C.  Heat  Exchanger  Performance 

D.  Medium  Quality,  Clearing  Time 

E.  Plasma  - Chemistry  Stability 

To  answer  some  of  these  issues,  fundamental  knowledge  of  fluid  and  thermal 
characteristics  of  the  recirculating  gas  flow  in  the  closed  cycle  circula- 
tor during  the  laser  operation  is  considered  essential  and  necessary  as 
one  of  basic  engineering  design  tools. 

To  the  best  of  our  knowledge,  experimental  data  and  theoretical 
understanding  directly  attributable  to  the  gain  of  fundamental  knowledge 
and  the  substantiation  of  technical  feasibility  are  lacking  for  the  re- 
circulating gas  flow  peculiar  to  the  pulsed  laser  operation.  At  present, 
some  of  the  major  technical  problems  associated  with  the  laser  develop- 
ment of  this  kind  are  related  to  the  questions  of  flow  and  thermal  stabi- 
lities of  the  recirculating  gas,  plasma  formation,  pressure  and  thermal 
waves  attenuations.  For  these  reasons,  a program  was  initiated  at  the 
Army  Missile  Command  to  design  and  fabricate  a unique  small  scale  closed 
cycle  gas  circulator  for  repetitively  pulsed  electric  discharge  lasers 
at  200°K  temperature  and  high  pressures  in  order  that  these  potential 
problem  areas  can  be  identified  and  investigated  at  relatively  low  cost 
and  risk.  Results  of  this  small  scale  program  may,  through  comparison 
with  other  programs  of  larger  scale  circulators,  contribute  to  the  estab- 
lishment of  a realistic  similarity  law  for  the  design  of  the  full-scale 
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circulator  at  a saving  of  weight  and  volume. 


As  a part  of  the  overall  team  effort  in  the  program  under  the  di- 
rection of  the  Army  High  Energy  Laser  Directorate,  the  UAH  team  has  ex- 
pended best-possible  efforts  in  the  following  areas  of  concern. 

(1)  Design  and  fabricate  a subscale  closed  cycle  laser  gas  cir- 
culator to  be  incorporated  with  the  Army's  single  pulsed  E-beam 
controlled  laser  to  conduct  simulation  experiments. 

(2)  Design  and  fabricate  an  acoustic  attenuator  (or  muffler)  and 
acoustic  reflectors  to  be  installed  in  the  closed  cycle  circulator 

in  order  to  study  the  wave  attenuation  and  interaction  characteristics. 

(3)  Develop  an  instrumentation  system  including  capabilities  to 
measure  and  record  pressure,  temperature  and  velocities  of  the  flow. 

(4)  Experimentally  investigate  fluid  and  thermal  characteristics 

of  the  recirculating  flow  associated  with  the  pulsed  wave  propagation. 

(5)  Experimentally  determine  the  effects  of  muffler  and  reflectors 
upon  the  wave  attenuation  of  both  acoustic  and  thermal  types. 

(6)  Study  experimental ly  the  effectiveness  of  the  blower  and 
simulating  heat  exchangers  as  interacting  elements  with  acoustic 
and  thermal  waves. 

(7)  Develop  computer  models  to  numerically  simulate  the  recirculating 
gas  flows  in  laser  operation,  including  steady  flow  case,  unsteady 
flow  case,  and  one-dimensional  wave  analysis. 

(8)  Assist  in  the  design  of  a laser  cavity  to  be  incorporated  with 
the  Rocketdyne-bui 1 t circulator  from  the  standpoints  of  gas  dynamics 
and  thermodynamics. 

To  the  best  of  our  knowledge,  the  experimental  results  reported  here 

are  the  first  to  show  the  transient  fluid  and  thermal  characteristics  of 
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a pulsed  electric  discharge  in  a functioning  closed  cycle  gas 
circulator.  These  results  should  be  of  interest  to  those  concerned 
with  the  design  and  operation  of  all  types  of  pulsed  gas  lasers 
either  closed  or  open  cycle  and  including  the  electric  discharge 
and  chemical  laser  systems. 
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Chapter  II 

NUMERICAL  SIMULATION  OF  THE  RECIRCULATING  AND  UNSTEADY 
GAS  FLOW  IN  LASER  OPERATION 

The  experimental  program  described  in  ChapterHI  was  supplemented 
by  a numerical  investigation  of  the  fluid-thermal  behavior  of  the  gas. 

The  numerical  simulation  work  involved  three  computer  programs.  The 
first  to  be  described  is  the  steady  state  program  which  is  a modifica- 
tion of  a computer  program  which  was  available  at  UAH  but  needed  modi- 
fications to  accommodate  the  low  temperature  operating  range.  This 
computer  program  is  designed  to  predict  the  steady  state  performance  of 
the  system  which  means  that  the  pulsed  operation  of  the  system  could 
only  be  simulated  by  assuming  an  average  power  input. 

The  second  computer  program  to  be  described  is  designed  to  simulate 
the  unsteady  flow  resulting  from  the  sudden  deposition  of  energy  in  the 
cavity.  This  program  is  capable  of  predicting  the  pressure,  temperature, 
and  velocity  waves  generated  at  the  cavity. 

The  third  computer  program  is  the  result  of  a one-dimensional 
analysis  of  the  unsteady  flow  generated  at  the  cavity.  This  program  was 
used  to  provide  estimates  of  the  pressure  and  temperature  environments 
expected  in  the  cavity  region. 

A.  Steady  State  Computer  Program 

The  basic  steady  state  computer  program  was  first  developed  by 
AVCO  for  application  to  the  MTU.  This  program  was  modified  at  UAH  under 
a previous  contract  in  order  to  accommodate  a different  heat  exchanger 
design  than  originally  employed  in  the  program.  Under  the  present  con- 
tract, this  program  was  modified  further  to  simulate  the  200  K gas  system 


under  consideration  in  this  work. 
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The  basic  computer  program  is  adequately  described  in  Reference  1 
and  will  not  be  covered  in  this  report.  A listing  of  the  program  with 
the  current  modifications  is  given  in  Appendix  A.  The  results  for  a 
typical  run  is  given  in  the  tables  which  follow.  Table  1 shows  the 
input  list  used  to  generate  the  results.  The  system  dimensions  are  those 
which  correspond  to  the  system  being  built  by  Rocketdyne.  The  computer 
program  calculates  the  fluid-thermal  characteristics  at  eleven  points 
around  the  system.  The  stations  are  located  as  follows: 


Station 

Location 

1 

Cavity  Entrance 

2 

Cavity  Exit 

3 

Diffuser  Exit 

4 

Constant  Area  Duct 

5 

Heat  Exchanger  (Room  Tet 

6 

Compressor  Inlet 

7 

Compressor  Exit 

8 

Constant  Area  Duct 

9 

Constant  Area  Duct 

10 

Heat  Exchanger  (200  K) 

11 

Heat  Exchanger  Exit 

The  gas  is  loaded  with  energy  at  the  cavity  at  the  rate  of  450  KW.  This 

corresponds  to  about  900  J/ e.  on  a continuous  basis.  The  results  for 

this  run  are  presented  in  Table  2. 

The  results  show  that  the  pressure  in  the  system  doesn't  vary  a 

great  deal  and  the  thermal  load  can  be  accommodated  on  a steady  state 

basis  with  the  heat  exchanger  loads  shown. 
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The  gas  properties  for  the  range  of  temperature  and  pressure  en- 
countered by  the  gas  in  traveling  about  the  system  were  investigated  care- 
fully to  ensure  that  large  errors  were  not  made.  No  major  difficulties 
were  encountered  since  the  pressure  remains  near  atmospheric  and  the 
temperature  is  above  200  K.  This  investigation  also  revealed  that  the 
system  could  likely  be  operated  at  even  a lower  temperature  without  danger 
of  CO^  snowing  out  of  the  gas.  This  can  be  seen  from  the  phase  diagrams 
for  the  laser  gas  mixtures  given  in  Appendix  B.  This  work  was  performed 
by  Mr.  David  Washington  as  an  exercise  in  an  advanced  thermodynamics 
course  taken  at  UAH. 

B.  Unsteady  Flow  Case 

1.  Theoretical  Considerations 

Fluid  and  thermal  characteristics  are  to  be  studied  along  the 
recirculating  gas  flow  path  as  a function  of  time  and  space.  As  the 
laser  pulses,  energy  will  be  input  into  the  laser  gas.  The  effect  of 
the  pulse  will  be  felt  along  the  fluid  loop.  Disturbances  in  the  form 
of  shock  waves  will  travel  out  from  the  cavity,  as  well  as  other  thermal 
changes  that  travel  with  the  fluid  particles. 

The  basic  equations  which  are  combined  into  the  governing  equations 
are  one-dimensional  fluid  mechanics  and  thermodynamics.  Together  with 
the  accompanying  boundary  conditions  and  initial  conditions,  these  equa- 
tions form  a mathematical  model  of  the  fluid  and  thermal  characteristics 


of  the  recirculating  laser  gas  flow. 

Based  on  the  principles  of  similitude,  the  pertinent  quantities  of 
the  flow  are  in  non-dimensional ized  form  and  presented  as  follows.  The 
continuity  equation  is  written  in  the  form 


1 9p  + 9p  + 9u  + _u_  9a 

p 9t  9x  9x  a 9x 


(1) 
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and  the  dynamical  equation  of  motion  is  given  as 


3u  , 3u  1 3P  fu|ul 

3t  u 3x  p 3x  ' 2D 

Combining  the  first  and  second  laws  of  thermodynamics  gives 

dS  = C , - R 

v T p 

The  rate  of  change  of  entropy  may  be  a function  of  position  and  time 
as  shown  implicitly  below. 


(2) 


(3) 


-jyp  = f (x,  t,  a,  u,  S)  (4) 

Normalizing  and  rearranging  the  above  equations  yield  a set  of  partial 
differential  equations  governing  the  transient  flow  in  one-dimension. 


= f U,  x.  A,  U,  S) 


(5) 


6t 


A+  uj  = - 


AU  + A |=£  + (y-1) 


DS 

Dt 


where 


sf—  * i A) 


3 

9 £ 


(6) 

(7) 


Since  two  of  the  governing  partial  differential  equations  are  hyper- 
bolic, tiey  possess  "characteristics,"  invariants  in  time.  At  each  point 
in  time  and  space,  there  are  two  invariants  moving  away  from  the  point. 

A numerical  solution  can  take  advantage  of  this  very  fact  and  is  thus 
called,  the  method  of  characteristics.  It  uses  invariants  from  nearby 
known  points  to  calculate  other  points.  These  characteristics  are  called 
P and  Q. 
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FIGURE  1 . DETERMINATION  OF  FLOW  CONDITIONS  AT  "LATER"  OR  AT 
"EARLIER"  POINT  FROM  GIVEN  INITIAL  CONDITIONS. 


FIGURE  2.  CHARACTERISTIC  DIAGRAM  FOR  THE  FLOW  ANALYSIS. 


Q 


A + U 


A - U 


(P  - Q) 


^ (P  + Q) 


= reference  length  (ft) 

= reference  speed  of  propagation  (ft/sec) 

= one-dimersional  space  coordinate  (ft) 

= time  (sec) 

= entropy  (Btu/lbm  - 'R) 

= flow  velocity  (ft/sec) 

= gas  constant 
= ratio  os  specific  heats 
= speed  of  sonic  wave  propagation  (ft/sec) 
= heat  transfer  (Btu/lbm) 

= pressure  (lbf/ft  ) 

= density  (slug/ft^) 

= temperature  (°R) 

= (y  R T)4 
= area  ratio 
= friction  factor 
= speed  of  sound 
= fluid  particle  velocity 


S = fluid  particle  entropy 

Q‘  = fluid  particle  heating 

£ = coordinate  distance 

t = coordinate  time 

P = right  moving  characteristic  of  governing  equation 
Q = left  moving  characteristic  of  governing  equation 


Thei  r 
Let  I 


respective  mathematical  definitions  are: 

■ yT  A±  U;  Na  . - AU  ^ ; Nb  . A;  N£ 


(9-DA  {£);  N 


D 


f Lq  U|U| 
2D 


ir  (I>  * na  + nb  -§r  <s>  * "c  + nd  <8> 

-S  (I)  = (Na  + Nc  + Nd)  <ir  + 6 (S)  Nb  (9) 

In  finite  difference  form,  Eq.  (9)  becomes 

A ' ■ 'VV  ND>a»e  *',4S  <NB>a,e  <’°> 

Let  the  characteristics  be  given  as 

P = — — y A + U (11) 

Y-l 

0 = A - U 02) 

Now  the  governing  equations  can  be  applied  to  the  numerical  solution  of 
the  unsteady  one-dimensional  flow  in  the  circulator. 


.J 
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2.  Numerical  Analysis 

To  apply  a numerical  scheme,  the  governing  equations  are 
written  in  the  finite  difference  form.  Equation  (5)  then  becomes  (13) 
and  Equation  (6)  produces  Equations  (14)  and  (15). 


DS 

Dt 


= a U.  i . A,  U,  S) 


03) 


(St 


= - AU 


9ilna  , „ 6+S 
9 r 6r 


+ (y-1)  A 


DS 

Dt 


- f 


Lo  UIUI 
2D 


(14) 


6 Q ~0„  6 S 

-f—  = - AU  + A -jf- 

6t  9f,  <5t 


DS  Lo  UIUI 
+ (y-1)  A -j£-  - f 


Dt 


2D 


(15) 


where 


P = A + U moving  at  | df~)  = ^ + A 


Q = A * u at  (-gf~) 


= U - A 


and  a particle  moves  at  = U. 

Two  known  points  are  used  to  calculate  a third  new  point,  located 
at  ( f; ^ » *3)*  The  p characteristic  from  (f^,  t-|)  intersects  with  the  Q 
characteristic  from  (f.^,  T2),  since  the  slopes  of  the  invariants,  called 
Riemann  variables,  are  known.  By  knowing  that  the  P and  Q are  constant, 
they  must  exist  at  (f,3,  t^).  Now  we  can  calculate  and  A^  by 

u3  ■ -r (p  - Q> 


A3  = (P  + Q). 
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Iteration  of  this  principle  continues  as  follows  with  Eq.  (16),  (17),  and 
(18)  representing  Eq.  (13),  (14),  and  (15),  respectively. 


S3  S4  + ( Dt  )43  (t3  " t4) 


(16) 


/3  an  (D,/D,)\ 

P3  = |"  A13  U13  \ E3  - ^ ) + 


(y_1)  A13  (St  )- 


. 3Lo 
13  2 °13 


U1 3 ' U1 3 I 


(c3  ' V + A1 3 ^S3  ' V + P1 


(17) 


A23  (St). 

U23  lU23  I J ^t3  " r2^  + A23  ^S3  " S2^  + 


/ 3 £n  (D,/D?) 

Q3  = |'  A23  U23  \ ^~C2 


f23Lo 

23  2 °23 


(18) 


Then, 


P3  ' Q3 


(19) 


A3  = (p3  + ^3^ 


(20) 


Test  for  iteration: 


Ui  - U3 


< 0.001 


(21) 


A3  ' A3 


0.001 


(22) 


where  the  primed  quantities  are  the  newly  calculated  ones. 

Once  the  method  of  characteristics  solution  is  completed,  the  thermo- 
dynamic quantities  can  be  calculated.  The  thermodynamics  equations  are: 

13 


d £n  p = 


da 


Y-l 


ds 


density 


(23) 


T = 


yRgr 


temperature  (24) 


P = p R T 


pressure 


(25) 


Equations  (23)  to  (25)  are  modified  according  to  the  numerical  scheme  of 
calculation  of  the  unknown  point  as 


)3  = p1  exp 


2 

y-l 


A, 


£n 


Y (S3  - 5,) 


density 


(26) 


T3  - V ' 11  "9c1 


„ . „ T / 32. 1 7 ]bm \ 

P3  - p3  R T3  \ T slug"') 


temperature  (27) 


pressure 


(28) 


Now  the  solution  scheme  moves  on  to  two  more  known  points  to  calculate 
a third.  Finally,  an  entire  new  line  is  calculated.  It,  in  turn,  is 
then  used  to  generate  another  line,  and  so  on. 

The  method  of  characteristics  only  breaks  down  when  it  must  cross  a 
shock  wave  or  contact  surface.  Any  discontinuity  requires  special 
attention.  The  solution  must  find  the  intersection  of  the  discontinuity 
with  the  characteristics  and  "jump"  across  to  a new  point  to  continue  as 
before.  This  "jump"  across  the  characteristics  is  accomplished  by  using 
Eqs.  (29)  and  (30)  which  are  derived  from  the  Rankine-Hugoniot  relations 
fcr  the  case  of  shock  waves  and  Eqs.  (31)  through  (34)  for  the  case  of 


contact  surfaces  as  illustrated  in  Figures  3 and  4. 
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lil  ill 

!il  = 1 + y-i  pi 


ill  + III 

7-1  Pj 


'II 


(29) 


II 


II 


JII 


(30) 


Q,  = Qd  + (pi  + Qtj)  tanh 


nl 


(s.  - sD) 


(31) 


PR  " PL  + QR  ‘ \ 


U2  = 


prA 


(32) 

(33) 


A'  = *ll 


T~  <PR  + V 


(34) 


Each  shock  wave  and  interface  must  be  continually  followed  and 
updated  as  the  solution  progresses.  The  "cut  off"  known  point  must  be 
replaced  by  a new  point  on  the  other  side  of  a discontinuity  to  continue 
the  method  of  characteristics. 

A flow  diagram  is  presented  in  Figure  5 to  illustrate  the  numerical 
procedure  for  the  solution. 

The  final  results  should  be  a location-time  history  of  the  velocity, 
pressure,  density,  and  temperature  of  the  unsteady  flow  in  the  circulator 
as  shown  below. 
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FIGURE  5.  FLOW  DIAGRAM  OF  THE  NUMERICAL  SOLUTION  OF  UNSTEADY  FLOWS 
USING  THE  METHOD  OF  CHARACTERISTICS. 


C.  One-Dimensional  Wave  Analysis 


Following  the  work  of  Foa  (Ref.  2),  a one-dimensional  analysis 
was  performed  on  the  waves  produced  by  the  energy  deposition  in  the 
cavity.  The  initial  conditions  are  assumed  to  result  from  an  instantaneous- 
constant-vol ume  heating  of  a region  of  gas  of  length  L in  a constant 
area  duct.  The  energy  deposited  goes  into  heating  the  gas  according  to 
the  following  equation 

■ cv 

Using  the  perfect  gas  law,  p = pRT,  and  using  the  fact  that  the  initial 
and  final  density  are  the  same,  we  obtain  the  following  relationship  for 
the  change  in  pressure 


A P _ q 
P C T 

O V 0 

where  q is  the  energy  per  unit  mass  deposited.  Thus,  the  equations  above 
can  be  used  to  find  the  initial  temperature  and  pressures  produced  by 
the  energy  pulse. 

The  over  pressure  at  the  boundaries  of  the  region  cause  a shock 
wave  to  propagate  away  from  the  cavity.  The  Mach  number  of  this  shock 
wave  is  determined  uniquely  from  the  pressure  ratio  using  normal  shock 
relations.  For  a pressure  ratio  of  P.j/P  , the  following  relationship 
holds 


P1  . 

2 k M - k + 1 

p 

0 

1 ' & (Ms  - -*;) 

2k 

(k+1 ) 

k-1 

The  value  of  M is  found  numerically. 
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The  wave  propagation  is  shown  sketched  in  Figure  6 for  flow  from  left 
to  right.  The  shock  waves  travel  in  both  directions  away  from  the  cavity. 
Behind  the  shock  is  a hot-cold  interface  which  initially  moves  in  the 
same  directions  as  the  shocks.  Behind  the  interface  is  an  expansion  fan 
which  moves  in  the  opposite  direction  of  the  shocks.  The  expansion  fans 
cross  at  the  center  of  the  cavity  (for  no  flow  case)  and  then  interact 
with  the  hot-cold  interface  which  was  generated  at  the  opposite  end  of 
the  cavity.  This  intereation  has  the  effect  of  stopping  the  outward 
movement  of  the  hot  gas. 

A computer  program  was  developed  to  calculate  the  shock  speed, 
pressure  ratios,  temperatures,  and  expansion  of  the  hot  gas  region. 

This  program  is  given  on  the  next  page.  The  results  are  given  in  the  table 
following  the  program.  The  notation  is  as  follows: 


QV 

= 

the  pulse  energy  in  J/8, 

U3MU1 

= 

velocity  in  region  3 minus  the  velocity  in  region  1 

FM 

= 

M , the  shock  Mach  number 

P1VP0 

= 

the  initial  pressure  ratio  p-j/P 

T1VT0 

= 

the  initial  temperature  ratio  Tj/T 

P3VP0 

= 

the  pressure  ratio  across  the  shock 

T3VT0 

= 

temperature  ratio  across  shock 

T2VT0 

= 

temperature  of  hot  gas  with  respect 

to  initial  temperature 

P0,P1, 

etc, 

= pressure  in  respective  region. 

N/m^ 

T0.T1, 

etc 

. = temperature  in  respective  region,  K. 

EXPAN 

= 

, a point  at  which  the  hot-cold 

interface  is  inter- 

sected  by  the  expansion  wave,  thereby  stopping  the 

expansion  of  the  hot  gas  region.  The  value  is  normalized 

with  respect  to  the  cavity  length,  L. 
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r- 


8v=n.o 

P_0»u^7 

ACK  = 0.000001 
T 0 = 2 0 0 • 

RH0C=1 .7656 
C SUBV*Ti 8 .285 
GAMMA=1 .4 
W AT  *2  2 . 4 

I F AC  = 1.15 

W9ITfc(C,100)  ACK, RH00,CSUl->V, GAMMA, WAT 
100  FORMAT (lHl,59X, ’LASER  S HO C K ’ , / , 6 0 X , ' 3 4 X , 

1 ’ A CK ' 7l2X , ' RHOO  * , 1 0 X , ' C S U 0 V ' , 1 0 X , ' G A MMA  ' ,11X,  *WAT  ' , / , / , 2 6 X , 

2 SEIS.b././././././.Wx.’QV'.llX.’UjMUl’.^X.’FM’.llx.’PlVPO’, 

3 1 OX , ' T1V  TO ' ,10X, ’ P 3 VP0 ’ ,10x,  'T3VTO’,10X, * T2  VTO ' , / , 1 2 X , ' P0  ' , 

4~  1 2 X i ' TO*  , 1 4 X , 'Pi'  ,12*  , * T 1 ' , 1 3 X , ' p 3 ' , 1 3 X , ' T 3 ' , 1 3 x , ' T 2 ' , 1 2 X , 

5 ' E X P A N ' ) 

00  2 0 J = 1 , 1 0 0 
' OW  = J *50  • 

QF  = Q V * F A C 

A0»SuRT(GAMKA*ii314.34*T0/WAT) 

PI VPO  = QF »1000./(RHOO*CSUOV*T0)*1  .0 

T 1 V TO  = P 1 VPO 

TG=2.*GAMMA 

GM1  =0AMMA-r.b 

6 P I *6  AMM  A ♦ 1 .0 

IF  M * 1 , 0 

9 FMS=$3RT(GM1  / TG*  (f.Pl /TG  ) *Tl  VTO*  ( I .0- (GM1 /GPI  > *SQRT  (1  . /T1  VTO)  * 

1 ( F M-1 ./ FM) ) ** ( T6/GM1 ) ) 

FCH=ABS(F,MS-FM) 

F M=  FMS 

I F ( FCh-ACK  >10,10,9 
10  PiVPO=(TG*FM*FM-GMl >/6P1 

T1WT2  = (PI  VP0/P3VPf)>  **  (GMl /GAMMA) 

T2VT0  = P1VP0/T1VT2 

T3VTC  = P3VPC*(1  • 0 ♦ ( 6 W 1 / G p 1 ) * P 3 V P 0 ) / (P3VPf'*GMl/GPl  ) 

A 1 « A 0 * $Q  R T (Pi  V P 0 ) 

T3=  T3VT0*TG 
P3=P3VP0*PU 
P'l  = T1VT0«P0 
T1=T1VT0*T0 
T2  * T?VTO*TQ 

u3MU1=(2.*A0/Gp1)*(FM-1./FM) 

EXPAN  = A 1 / ( A 1 -UJMOI  ) -1  . 0 

WRITb(G,101>  av,U3MUl,FM,P1VPO,TlVTO,P3VPO,T3VTO,T2vTO,PO,TO, 

II  >1  , T 1 ,P3,T3,T2, EXPAN 

101  F0HmaT(1hC,2x,<JE15.6,/,3x,8e15.6> 

20  CONTINUE 
STOP 
END 
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NORMALIZED  TIME,  T»  ta0/Lc 


FIGURE  6. 


WAVE  PROPAGATION  AND  INTERACTION  BASED  ON  ONE-DIMENSIONAL 
ANALYSIS.-1  ^ 


D.  Muffler  Design  Considerations 

Mufflers  are  described  as  devices  to  reduce  the  transmission  of 

sound  while  at  the  same  time  allowing  the  free  flow  of  gas  in  a duct  or 

pipe.  An  ideal  muffler  for  a pulsed  laser  system  would  be  a device  which 

would  reduce  considerably  the  transmission  of  the  shock  waves  produced 

by  the  high  energy  deposition  while  at  the  same  time  allow  for  the  free 

unobstructed  flow  of  the  laser  gas.  Mufflers  of  this  type  employ  holes 

in  the  side  of  the  duct  which  connect  to  a closed  region.  Through  proper 

design,  the  pressure  wave  energy  can  be  directed  from  the  gas  flow,  through 

the  side  holes,  and  into  the  dissipation  region  without  disturbing  the 

steady  flow  of  the  gas  in  the  duct.  Mufflers  of  this  type  are  called 

volume  resonators.  Parameters  of  importance  in  the  design  of  such  mufflers 

are:  the  frequency  of  waves  to  be  attenuated,  the  volume  of  the  resonator, 

the  size  of  the  holes  in  the  duct,  and  the  length  (if  other  than  zero)  of 

4 

the  tubes  which  connect  the  holes  with  the  resonator  volume.  As  is  typical 
of  muffler  design,  the  geometric  parameters  are  dependent  upon  the  fre- 
quency of  sound  to  be  attenuated.  For  example,  the  required  resonator 
volume  is  inversely  proportional  to  the  frequency.  For  this  reason,  high 
frequency  sound  can  often  be  attenuated  by  simply  making  the  duct  walls 
out  of  acoustic  tile  or  other  commonly  employed  sound  proofing  materials. 
Low  frequency  pressure  waves  require  volumes  much  larger  than  can  be 
accommodated  in  wall  coverings.  Thus,  for  low  frequency  sound,  the  duct 
must  have  openings  along  the  side  which  connect  to  the  large  volume 
resonator. 

As  discussed  in  previous  sections,  the  ideal  one-dimensional  pressure 
wave  generated  by  the  laser  pulse  will  appear  as  a square  wave  pressure 
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pulse  traveling  through  the  system.  The  square  wave  can  be  approximated 
by  an  infinite  series  of  sine  and  cosine  waves.  Thus,  the  muffler  must 
be  able  to  attenuate  waves  over  a wide  band  in  frequency  to  be  effective 
in  pulsed  laser  operation.  In  order  to  make  the  muffler  effective  over 
a wide  band,  it  was  decided  to  couple  the  duct  to  the  resonator  volume 
through  the  use  of  horn  amplifiers  rather  than  the  tubes  normally  used. 
The  straight  tube  design  would  be  effective  at  only  one  frequency  while 
the  horn  has  the  potential  for  a wide  band  of  frequencies. 

The  ideal  exponential  horn  amplifier  has  the  characteristic  of  a 
high  pass  filter  with  no  reflection.  The  geometry  of  the  area  variation 
of  the  ideal  horn  is 


A(x)  = Ao  e 

where  A(x)  is  the  area  at  a position  x away  from  the  throat,  Aq  is  the 
area  at  the  horn  throat,  and  B is  the  exponential  factor.  The  cut  off 
frequency  of  the  horn  is  determined  by  the  value  of  B and  the  total 
length  of  the  horn.  For  zero  reflection,  the  ideal  horn  must  be  in- 
finite in  length.  For  a finite  length  horn,  the  reflected  energy  is 
frequency  dependent,  with  the  reflection  increasing  for  the  lower  fre- 
quencies. For  example,  if  the  exponential  factor  B is  such  that  the 
radius  of  the  horn  of  circular  cross  section  at  the  exit,  r , were  to  be 

2 

re  B 

then  the  reflection  at  w/B  = 4 x 104  is  10%  but  at  w/B  = 2 x 10^  the 
reflection  is  30%.  If,  however,  the  relationship  were 


r 


e 


1 

B 


the  reflection  at  w/B  = 4 x 10^  is  25%  and  at  w/B  = 2 x 104  the  reflection 

is  55%.  Thus,  the  horn  geometry  for  a finite  horn  has  considerable  effect 
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on  the  filtering  and  reflection  capabilities  of  the  horn.  The  table  below 
gives  a more  complete  picture  of  the  geometric  influence  on  the  horn 
characteristics. 

For  pulsed  laser  operation  at  200  pulses  per  second  maximum,  all  of 
the  energy  in  the  resulting  pressure  waves  is  expected  to  be  in  the  fre- 
quency range  above  200  Hz.  Thus,  a horn  designed  for  200  Hz  should  re- 
move from  the  flow  all  the  pressure  wave  energy  in  the  200  Hz  and  above 
range  since  a horn  acts  as  a high  pass  filter.  Consider  a design  in 

4 

which  cj/B  = 4 x 10  and  rg  = 2/B  which  gives  a 5%  reflection  as  seen  in 
the  table.  If  this  design  were  to  be  effective  at  200  Hz, 

B = 1 257/4x1 04  = .037  cm'1 


and 


r = 63.66  cm 
e 


The  length  of  the  horn  is  given  by 

, Sin  ( Ae/Ao  ) 
L ~ B 


where  the  length  is  seen  to  be  dependent  upon  the  throat  area  Ao.  Take 
the  throat  area  to  correspond  to  a radius  of  1 cm  and  we  obtain 

L = = 276-9  cm 

which  shows  that  a horn  9 ft  long  is  needed  to  be  95%  effective  at  200  Hz. 
Ap  effective  length  of  9 ft  could  likely  be  obtained  in  a folded  horn 
design  which  would  occupy  a reasonable  volume  and  thus  be  applicable  for 
pulsed  laser  systems.  Therefore,  theoretically,  the  horn-coupled-resonator 
muffler  appears  to  have  the  potential  for  the  wide  band  performance  re- 
quired in  the  pulsed  laser  systems. 
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Table  of  Percent  Reflection  as  a Function  of  oj/B  for  Three  Exponential 

Horn  Shape  Factors 


sec  r 

\e 

1/B 

4/3B 

2/B 

2 x 10? 

55% 

40% 

33% 

4 x 10? 

25% 

12% 

5% 

6 x 10? 

15% 

8% 

3% 

8 x 10? 

11% 

5% 

2% 

16  x io? 

6% 

3% 

1% 

24  x 10? 

4% 

2% 

0% 

32  x 10q 

3% 

1% 

0% 

A horn-coupled-resonator  was  designed  and  built  for  the  purpose  of 
investigating  experimentally  the  performance  of  such  a muffler  in  the  UAH 
subscale  recirculation  system.  Available  materials  were  employed  which 
limited  the  length  of  the  horn  to  6.5  cm.  The  horn  was  designed  to  be  two 
dimensional  with  a rectangular  cross  section  rather  than  circular.  The 
width  at  the  opening  was  limited  to  11.5  cm  plus  the  throat  width.  Thus, 
with  these  geometric  constraints,  the  horn  expansion  factor  was  restricted 
to  be  approximately  .3  for  values  of  throat  width  of  about  1 to  2 cm.  For 
example,  for  a throat  width  of  1.5  cm. 


Z n Ae/Ao 
L 


2.n  (13/1.5) 
6.5 


For  this  value  of  B,  the  value  of  w/B  for  200  Hz  becomes, 
id  „ 1252  ,^4 


at  1000  Hz,  this  becomes 


.42  x 10 


2.1  x 10 


and  at  2000  Hz, 

= 4.2  x 104 

Thus,  we  expect  the  horn  described  above  to  be  effective  for  frequencies 
above  1000  Hz  based  on  the  m/B  ratio  and  reference  to  the  table  on 
percent  reflection. 


Chapter  III 

EXPERIMENTAL  INVESTIGATION  OF  SUBSCALE  CLOSED  CYCLE 
LASER  GAS  CIRCULATOR 

3 

A recirculation  system  was  designed  for  use  with  the  MICOM  S single 
pulse  laser.  This  system  was  built  for  the  purpose  of  investigating 
experimentally  the  transient  flow  and  heating  problems  characteristic  of 
a high  energy  laser  system.  The  equipment  employed  in  this  work  is  de- 
scribed in  Section  A of  the  following  which  includes  a description  of  the 
reci rculation  apparatus,  a muffler  used  in  wave  attenuation  tests, 
acoustic  reflectors,  and  the  instrumentation  employed  in  the  experimental 
program.  The  experimental  data  collected  by  this  system  will  be  described 
in  general  in  Section  B. 

A.  Description  of  Experimental  Equipment 

1.  Recirculating  Flow  Apparatus 

The  recirculation  system  consisted  of  a blower,  PVC  piping, 

3 

mounting  flanges  for  attachment  to  the  S laser,  an  E-beam  mask,  and  a 

3 

plastic  structure  which  was  mounted  inside  the  S cavity.  A schematic 

3 

of  the  system  attached  to  the  S laser  is  given  in  Figure  7.  The  flow 
through  the  cavity  region  is  from  right  to  left  in  that  figure, 
a.  Blower  Specifications 

The  blower  is  a Dayton  Model  4C108  with  a 10  5/8"  wheel. 

The  inlet  is  6V  in  diameter.  The  outlet  is  2\"  square.  The  outside 
of  the  blower  casing  is  approximately  13"  in  diameter. 

The  motor  is  connected  directly  to  the  blower  wheel.  The  motor  is 
a 1 hp  Dayton  Model  6K232  operated  at  3450  rpm. 
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b.  PVC  Piping 


The  majority  of  the  piping  consisted  of  3109662  "Certain- 
Teed  Plastics"  4"  PVC  1120  Schedule  40  220  psi  ASTMD-2665.  This  pipe 
has  4"  ID  and  4V ' OD.  The  90°  turns  were  achieved  using  standard  PVC 
fittings  which  have  a 4"  radius  of  curvature  along  the  center  line  of 
the  90°  bend.  A number  of  couplings  were  employed  throughout  the  path 
length  to  accommodate  assembly,  disassembly,  and  modification  during  the 
testing. 

Also  shown  in  Figure  7 is  a cross-over  or  by-pass  pipe  with  valve. 
This  pipe  was  J24G264  1"  PVC  1120  Schedule  80  630  psi  at  73°F  ASTMD  1785. 
The  pipe  was  1"  ID  and  IV  OD  and  was  attached  to  the  4"  PVC  through 
a standard  4"  to  2"  tee  with  a reducer  used  to  bring  the  size  down  to 
1".  A 1"  gate  valve  was  used  in  the  1"  line. 

c.  Dimensions 

The  piping  dimensions  are  shown  on  Figure  8 in  centimeters. 
Total  path  length  around  the  4"  circuit  is  approximately  1140  cm  (37  ft.) 

d.  Flanges 

Special  flanges  had  to  be  constructed  in  order  to  attach 

3 

the  PVC  pipe  to  the  blower  and  to  the  S laser.  For  the  blower  attach- 
ment, a 6"  to  4"  reducer  was  used  on  the  inlet  side  and  an  aluminum 

piece  was  machined  to  go  from  the  square  to  the  circular  cross-section. 

3 

The  flanges  for  attachment  to  the  S were  constructed  from  alum- 
inum consisting  of  4b"  diameter  holes  with  suitable  bolt  circle  and  o- 
ring  groove. 

e.  E-beam  Masks 

The  purpose  of  the  E-beam  mask  is  to  contain  the  E-beam  to 

within  a specified  region  of  the  cavity.  Figure  9 shows  the  dimensions 

of  the  three  masks  employed  in  this  work. 
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MEASUREMENTS  ON  PYC  MODEL  WITH  MUFFLER  ADDED 


f.  Cavity  Structures 


The  cavity  structures  serve  to  box-in  the  discharge  region 
in  the  vertical  and  longitudinal  directions  and  also  to  provide  a 
support  for  the  PVC  pipe.  Two  structures  were  employed  in  this  work; 
the  15  cm  structure  was  used  with  the  15  cm  mask  while  the  100  cm  structure 
was  used  with  both  the  45  cm  and  100  cm  masks. 

A sketch  of  a cavity  structure  is  shown  in  Figure  10.  Both  the  15 
and  100  cm  structure  were  constructed  in  the  same  manner  with  only  the 
dimension  of  the  opening  being  different. 

2.  Muffler 

In  consultation  with  Mr.  Cason  of  MIC0M,  a muffler  was  de- 
signed, constructed,  and  tested  in  the  recirculation  system.  The  location 
of  the  muffler  in  the  system  is  shown  on  Figure  8.  The  muffler  has  no 
obstructions  to  the  flow  and  was  constructed  to  attach  easily  to  the  4" 

PVC  system.  The  physical  dimensions  of  the  muffler  are  shown  in  Figure  11. 

3.  Reflectors 

The  reflectors  tested  are  shock  wave  attenuating  and/or 
reflecting  materials  with  frontal  configurations  of  honeycomb  and  square. 
The  honeycomb  is  made  of  aluminum  and  the  square  is  made  of  THERMAC0MB 
alumina  ceramics.  Table  I presents  the  measured  and  calculated  physical 
properties  of  the  materials. 

Table  I 


Type  of 
Openings 

Specimen 

Thickness(in) 

No.  of  Openings 
per  1 inear  in . 

Wall 

Thickness(in) 

Percent 
Open  Area 

Hyd . Rad. 
per  Open 

Honeycomb 

2.0 

4.5 

0.00075 

92.86 

0.0245 

Square  I 

0.5 

9 

0.0044 

85.14 

0.0146 

Square  II 

1.068 

7.5 

0.0090 

74.29 

0.0170 
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r 

/ Area  \ 

1 Calculated  values  of  the  hydraulic  radius  ( . — r — : — r — ) and 

yWetted  Perimeter/ 

I percent  open  area  are  based  on  an  average  of  the  appropriate  measured 

1 values.  Frontal  appearances  of  these  gratis  materials  are  shown  in 

Figure  12.  These  reflectors  are  placed  in  the  flow  apparatus  for  a series 
of  tests  to  determine  the  effectiveness  of  wave  attenuation. 
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Figure  1 


Frontal  Appearan  e 


Reflectors 


4.  Instrumentation  System 
a.  Pressure  Measurement 


The  instrument  used  for  pressure  measurements  is  a battery 
of  four  Piezotron  pressure  sensing  units  each  of  which  consists  of  a 
piezoelectric  pressure  transducer  (Type  201B5)  and  a coupler  (549R)  con- 
nected with  a 128M  cable.  The  unit  is  then  connected  to  an  oscilloscope 
for  the  readout  of  voltage  signals.  Specifications  of  the  Piezotron 
miniature  pressure  sensor  or  transducer,  the  Piezotron  coupler,  and  the 
Textronix  oscilloscope  are  presented  in  Appendix  D. 

The  pressure  of  up  to  100  psi  was  sensed  by  the  mini-gage  which  gives 
a direct,  high  level,  voltage  signal  with  less  than  100  ohms  output  im- 
pedance and  high  frequency  response  of  50  K HZ  and  low  frequency  response 
of  0.005  HZ.  The  sensor  then  converts  the  pressure  into  electrical 
voltage  with  bias  of  up  to  11  ^2  volts.  The  power  required  by  the  trans- 
ducer to  operate  is  supplied  by  the  coupler,  and  the  signal  from  the 
transducer  to  the  readout  equipment  is  transmitted  through  the  coupler 
over  a single  inexpensive  cable.  This  eliminates  all  of  the  inherent 
piezoelectric  high  impedance  problems  of  electrical  leakage,  cable  noise 
and  signal  attenuation  and  allows  the  transducers  to  be  used  in  contami- 
nated environments  and  with  long  and  moving  cables  at  low  noise  and  with- 
out use  of  charge  amplifiers. 

The  calibration  of  the  transducers  was  performed  at  the  factory, 
and  the  values  of  the  calibration  were  noted  to  be,  on  the  average  for 
all  probes,  50  mv  per  psi  for  the  pressure  measurement  up  to  100  psi. 


I 
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b.  Temperature  Measurement 

Due  to  the  extremely  transient  nature  of  temperature  varia- 
tion in  the  recirculating  flow  as  a result  of  the  pulsed  laser  operation, 
a sensor  of  high  frequency  response  in  excess  of  1 K HZ  is  considered 
necessary  for  the  temperature  measurement.  Search  of  an  adequate  sensor 
resulted  In  the  selection  of  a hot-wire  sensor  made  of  0.00015  in  diameter 
tungsten  wire  coated  with  platinum  powdery  film.  The  hot-wire  sensor  is 
connected  to  the  Temperature  and  Switching  Module  (Thermo-Systems  Model 
1040)  which  is  in  turn  connected  to  the  power  supply  (Model  1031-10A). 

The  Module  consists  of  a bridge  circuit  and  amplifier  in  an  open 
loop  configuration  so  the  hot-wire  sensor  which  is  ordinarily  used  as 
an  anemometer  probe  can  be  switched  to  function  as  a resistance  thermo- 
meter. Since  there  Is  a linear  proportional ity  between  the  voltage  out- 
put and  the  temperature,  the  calibration  can  be  simply  performed  by  ad- 
justing the  zero  and  gain  set  potentiometers  to  a desired  temperature 
range  using  the  calibrate  pots  of  two  temperatures. 

c.  Velocity  Measurement 

For  the  measurement  of  velocities,  hot-wire  probes  the  same 
as  those  used  for  the  temperature  measurement  is  applied.  The  probe  is 
connected  to  the  constant  temperature  anemometer  module  (Model  1010A). 

The  amplified  output  signal  from  the  anemometer  is  sent  to  the  Linearizer 
(Model  1005B)  so  that  the  voltage  signal  is  processed  in  such  a way  that 
it  became  linearly  related  to  velocity  of  the  gas  flow. 

The  use  of  these  modules  ensures  the  frequency  response  above  500 
K HZ  with  power  output  as  high  as  1.5  amps.  The  noise  associated  with 
the  anemometer  is  noted  to  be  less  than  0.007%  equivalent  turbulent  in- 
tensity. Frequency  response  to  the  Linearizer  is  found  to  be  up  to  400  K 

HZ  and  the  accuracy  of  linearization  can  reach  + 0.2%. 
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With  these  special  features  of  the  instrument,  it  is  able  to  measure 

both  average  velocity  and  turbulence  in  one-dimension. 

Calibration  of  the  probe  is  performed  by  using  a Thermo-Systems 

Calibrator  (Model  1125)  in  accordance  with  the  furnished  instructions. 

The  readout  system  for  both  temperature  and  velocity  is  the  Tektronix  type 

oscilloscope  (Type  564-3A74-3B3) . 

B.  Presentation  of  Experimental  Data 

An  experimental  run  typically  involved  firing  the  S laser  and 

recording  pressure,  temperature  or  velocity  at  various  locations  in  the 

recirculation  system.  Besides  the  pressure,  temperature  and  velocity 

data,  the  laser  discharge  current  was  also  monitored  in  order  to  deter- 

3 

mine  the  power  put  into  the  gas.  Each  firing  of  the  S laser  was  logged 
with  respect  to  the  sequence  number.  The  first  UAH  experiment  began  with 
sequence  number  3043  on  December  7,  1976  and  the  last  experiment  was  with 
shot  number  3238  on  February  2,  1977.  The  log  of  all  the  UAH  shots  is 
given  in  the  table. 

The  tests  were  designed  to  investigate  the  fluid-thermal  disturbances 
produced  by  the  sudden  energy  deposition  in  the  cavity  region.  In  the 
following  sections,  we  describe  the  wave  forms  produced  by  the  laser, 
the  wave  interactions  with  the  system  components,  and  the  wave  attenuations 
within  the  system.  Each  of  these  topics  will  be  divided  into  the  pressure 
wave  results  and  the  thermal-velocity  wave  results  if  appropriate. 

1.  Wave  Forms  Produced  by  the  High  Energy  Laser 

3 

The  S laser  was  operated  as  normal  except  that  no  attempt 

3 

was  made  to  produce  lasing  action.  The  S laser  merely  served  to  deposit 
energy  into  the  gas  at  the  cavity  section,  thereby  simulating  the  energy 
pulse  that  takes  place  in  such  a laser.  The  E-beam  was  operated  at  90  KV, 
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LOG  OF  UAH  SHOTS 


Reference 

Number 

Date  of 
Tent 

C.is 

Used 

H 1 ower 
Slat  us 

liner  gy 
Pulse , 

3043 

12/7/76 

N2 

Off 

3044 

M 

i I 

On 

41 

3045 

II 

•• 

" 

26.45 

3046 

If 

II 

Off 

10.  11 

304  7 

1 1 

II 

• 1 

28.22 

J048 

»» 

It 

28.27 

J053 

12/8/76 

M 

Of  f 

3054 

• I 

H 

On 

30  5 8 

12/9/76 

II 

- 

9.46 

3050 

H 

II 

On 

15.51 

3060 

M 

I# 

II 

18.45 

306) 

• 1 

II 

II 

15.81 

3062 

fl 

II 

II 

13.88 

3064 

M 

» 

II 

14.67 

3065 

• 1 

.. 

10.11 

3066 

II 

II 

II 

io.  n 

306  7 

II 

• 1 

II 

1 4 . 6 1 

3068 

• I 

II 

II 

14.67 

3069 

If 

If 

- 

1 l . 50 

30  70 

If 

" 

of  1 

1 l .00 

30/1 

If 

II 

On 

13.32 

3075 

12/10/76 

Ar 

Off 

56.65 

3076 

If 

•i 

On 

36.65 

of  Test 


Pressure  wave  propagation  with 
no  obstructions 


Pressure  wave  propagation  with 
large  pore  metal  honeycomb 
placed  at  heat  exchanger  locations 

Pressure  & temperature  wave  pro- 
pagation tests  with  large  pore 
metal  honeycomb  at  heat  exchanger 
1 oca t ions 
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Reference 

Number 

Date  nt 
Test 

(las 

Used 

It  I owe  r 
St at  us 

laierp.y  of 
Pulse,  J 

Purpose  of  Test 

3079 

It 

• • 

II 

8.20 

Pressure  wave  tests  with  two 
27  mm  sections  of  cer  honeycomb 
between  stations  3 and  4 

3080 

H 

• 1 

Off 

9.24 

II 

3082 

11 

11 

II 

7.04 

3083 

11 

'• 

II 

8.80 

• 1 

JO8/4 

H 

lie 

• 1 

8.67 

" 

? 185 

II 

II 

On 

9.35 

3086 

II 

II 

12.97 

308  7 

II 

• 1 

it 

12.97 

,, 

3088 

• 1 

II 

Off 

19.59 

3089 

It 

• 1 

II 

12.07 

3090 

II 

Ar 

On 

30.12 

1 f 

3092 

1 1 

26.60 

309  3 

II 

1 1 

Ol  f 

22.63 

II 

3094 

II 

'• 

• 1 

26.10 

II 

3095 

It 

If 

On 

27.9 

" 

3096 

12/14/76 

II 

• • 

39 . 2 

Pressure  wave?  propagation  tests 
with  4 sections  of  cer  honeycomb 
between  stations  3 and  4 

309  7 

•' 

• 1 

II 

25.66 

•i 

3098 

II 

M 

on 

17.87 

•• 

3099 

It 

•' 

H 

13.32 

it 

3100 

II 

•• 

• I 

9.75 

1 1 

3101 

12/17/76 

On 

19. 76 

Pressure  wave  propap.at  Ion  tests 
with  muffler  placed  between 
stations  3 and  4 

3102 

II 

II 

II 

21.88 

II 

. > 

i 

. — 

3103 

II 

II 

• 1 

12.00 

43 

II 

A 

Kefe rente 
Number 


blower 

Status 


Energy  of 
I’ulso,  J 


Purpose  of  Test 


rN 


rs 


Da t e of 
Test 


(Ms 

Used 


3131 

12/22/76 

Ar 

Off 

91.87 

3132 

II 

n 

• 1 

83.33 

3133 

M 

n 

II 

79.12 

3134 

II 

- 

If 

98.58 

3135 

II 

H 

On 

66.75 

3136 

II 

H 

• 1 

52.13 

3137 

II 

H 

tf 

40.13 

3138 

II 

n2 

• 1 

53.67 

3139 

II 

II 

II 

35.00 

3140 

II 

II 

Off 

48.42 

3141 

ft 

Ho 

On 

55.73 

3142 

II 

II 

II 

45.68 

3143 

II 

II 

Off 

52.70 

* 

Enlargements  made  of 

scope 

t races 

for  purpose 

Pressure  wave  tests  with  new 
probe  location  near  muffler  and 
time  delayed  signals 


data  analysis 


I 
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Kef  ercnce 

Date  of 

das 

1!  lower 

Energy 

of 

Purpose  of  Test 

Number 

Test 

Used 

Status 

Pulse , 

.1 

3144 

1/7/77 

N2 

Off 

All  four  pressure  probes  put 

3145 

M 

50.0 

at  same  location  to  test 
relative  responses. 

3146 

H 

49.9 

3147 

H 

53. 

3150 

1/25/77 

Argon 

51.3 

Installed  new  probe  at  up- 

3151 

M 

II 

32.0 

stream  of  discharge  to  test 
effect  of  corner. 

3152 

It 

II 

On 

Also,  thermal  wave  tests. 

3153 

II 

II 

II 

30.7 

3154 

II 

II 

II 

30.8 

3155 

1/26/77 

N2 

H 

4 7.6 

3156 

II 

II 

• I 

3157 

1/27/77 

II 

II 

23.1 

3158 

II 

17.8 

Corner  and  thermal  wave  test 

3159 

II 

N2 

Off 

17.5 

in  N2. 

3160 

II 

II 

II 

Thermal  wave  at  different 

3161 

It 

II 

On 

29.0 

locations  around  system. 

3162 

ft 

If 

m 

34.8 

3163 

II 

II 

it 

22.7 

3164 

II 

II 

ii 

26.9 

3165 

II 

II 

34.0 

Wave  near  muffler. 

3166 

II 

II 

On 

26.5 

3167 

II 

II 

II 

25.3 

3168 

If 

He 

II 

52.8 

Wave  propagation  and  thermal 

3169 

II 

II 

II 

18.0 

wave. 

3170 

II 

• 1 

II 

20.  7 

3171 

II 

II 

II 

38.  5 

3172 

1/28/77 

N2 

II 

62.0 

Muffler  incident  and  reflection 

3173 

II 

53.6 

test  and  thermal  wave. 

3174 

II 

N2 

Off 

39.6 

3175 

II 

ft 

31.2 

3176 

II 

II 

Off 

25.3 

A LI, 

OK  THE  ABOVE 

TESTS 

USED  A 15 

CM  CAVITY  REGION 
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:ference 

Date  of 

(Ins 

K 1 owor 

Energy  of 

Purpose  of  Test 

Number 

Test 

Used 

S tatus 

I’ulsc,  .1 

THE 

FOLLOW  1 Nti 

; TESTS  USE 

A 100  CM  CAVITY 

REGION 

3177 

2/1/77 

N2 

>>126 

Muffler  incident  reflection 

tests  and  thermal  wave  tests. 

3178 

M 

II 

On 

208 

3179 

II 

II 

II 

312 

3180 

II 

II 

182 

3181 

II 

II 

On 

186 

3182 

II 

300 

Corner  test  and  thermal  wave. 

3183 

• 1 

N2 

303 

3184 

II 

II 

308 

3185 

II 

II 

On 

280 

3186 

II 

II 

it 

288 

3187 

II 

II 

Off 

303 

Corner  test  with  probes 
reversed . 

3188 

II 

II 

On 

325 

3189 

II 

II 

it 

221 

Muffler  transmission  test. 

3190 

II 

165 

3191 

II 

N2 

On 

187 

3192 

II 

II 

it 

194 

3193 

II 

II 

H 

184 

3194 

II 

II 

H 

3195 

II 

II 

H 

294 

3196 

II 

II 

H 

128 

3198 

II 

II 

290 

Muffler  reflection  and 
incident  test 

3199 

II 

3200 

It 

N2 

On 

3201 

2/2/77 

Arson 

On 

133 

3202 

II 

It 

it 

277 
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Reference 

Date  of 

Gas 

Blower 

Energy  of 

Purpose  of  Test 

Number 

Test 

Used 

Status 

Pulse,  J 

3203 

2/2/77 

Argon 

280 

3204 

II 

II 

On 

232 

3205 

II 

II 

II 

242 

3206 

II 

II 

• 1 

294 

Muffler  transmission 

3207 

• 1 

II 

II 

239 

3208 

II 

It 

305 

3209 

II 

II 

260 

Thermal  wave  tests 

3210 

II 

400 

3211 

II 

Argon 

On 

390 

Thermal  wave  and 

3212 

II 

II 

H 

347 

corner  test 

3213 

II 

i: 

•i 

372 

Thermal  wave  and  Fan 

3214 

II 

ii 

H 

260 

Reflection  Test 

3215 

II 

ii 

•i 

371 

3216 

M 

ii 

n 

307 

3217 

II 

ii 

322 

Fan  transmission  test 

3218 

II 

ii 

300 

and  thermal  wave  test 

3219 

II 

it 

289 

3220 

II 

ii 

350 

48 
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Reference 

Date  of 

das 

11  Lower 

Energy 

of 

Purpose  of  Test 

Number 

Test 

Used 

Status 

Pulse , 

J 

3221 

2/2/77 

Argon 

146 

Fan  transmission  test  and 
thermal  test 

3222 

H 

II 

154 

3223 

II 

II 

158 

3224 

II 

II 

183 

Fan  reflection  and 
incident  test 

3225 

II 

II 

128 

3226 

H 

II 

On 

172 

Corner  test  and  thermal 

wave  tests 

3227 

• 1 

II 

II 

154 

3228 

II 

II 

II 

135 

3229 

II 

II 

II 

131 

Muffler  reflection  and 
incident  test 

3230 

II 

II 

II 

197 

3231 

II 

II 

II 

179 

3232 

II 

II 

H 

171 

3233 

It 

II 

II 

203 

Muffler  transmission  test 

3234 

II 

II 

II 

168 

3236 

II 

II 

II 

205 

* 

3237 

II 

II 

II 

104 

Corner  test 

3238 

II 

II 

II 

122 

Til K AUOVE 

TESTS  USE 

A 45  CM 

CAVITY 

REGION 

the  sustainer  voltage  was  varied  between  15-30  KV,  and  the  pulse  width 
was  varied  between  2 and  3 p sec.  As  pertains  to  the  fluid  dynamics  of 
the  cavity  flow,  the  energy  deposition  is  instantaneous  and,  therefore,  the 
process  can  be  described  as  an  instantaneous  constant  volume  heating 
process.  The  gas  pressure  and  temperature  are  instantaneously  increased 
by  the  energy  deposition  while  the  density  remains  unchanged.  The  fol- 
lowing describes  waves  found  produced  by  these  initial  conditions, 
a.  Pressure  Wave  Form 

The  pressure  wave  produced  by  the  energy  deposition  con- 
sists of  a sharp  pressure  rise  followed  by  a sharp  pressure  drop  and  a 
negative  phase.  This  is  illustrated  in  Figure  1 3 where  the  pressure  shown 
is  positive  in  the  negative  direction  since  this  is  the  way  the  data 
from  the  pressure  transducers  are  recorded  on  the  scope  pictures.  The 
width  of  the  positive  pressure  pulse  is  typically  1 m sec  and  the  nega- 
tive phase  has  about  the  same  width.  The  amplitude  of  the  positive  phase 
was  a function  of  the  energy  deposited  in  the  gas  which  varied  between 
7 and  100  J for  the  15  cm  cavity,  between  120  and  400  J for  the  100  cm 
cavity,  and  between  100  and  205  for  the  45  cm  cavity.  (See  log).  The 
amplitude  of  the  negative  phase  was  generally  proportional  to  the  ampli- 
tude of  the  positive  phase.  Figures  14,  15, and  16show  a typical  pressure 
wave  for  the  15,  45,  and  100  cm  cavity  lengths  respectively.  Figure  14 
for  the  15  cm  cavity  shows  the  wave  for  the  three  different  gases  used 
in  those  tests  (N^,  Ar,  and  He).  Only  Ar  gas  was  used  in  the  45  cm  tests 
and  a typical  wave  form  is  shown  in  Figure  15.  Both  and  Ar  were  used 

in  the  100  cm  tests  and  these  results  are  shown  in  Figure  16.  The  100 
cm  wave  form  shows  two  additional  positive  pressure  pulses  following  the 
negative  phase  of  the  wave.  This  feature  is  not  present  in  the  15  or  45 


cm  results. 
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All  the  results  shown  in  Figures  14,  15,  and  16  are  for  the  case  with 
the  fan  on.  The  character  of  the  wave  generated  by  the  laser  with  the 
fan  off  is  shown  in  Figure  17  for  the  15  cm  case  and  He  gas. 
b.  Thermal  Wave  Form 

The  pressure  waves  described  above  travel  away  from  the 
cavity,  both  upstream  and  downstream,  at  speeds  greater  than  sonic.  The 
pressure  in  the  cavity  region  is  adjusted  back  to  a value  near  ambient 
in  a short  period  of  time.  The  temperature  discontinuities  created  by 
the  energy  deposition  do  not,  however.  A hot  region  of  gas  remains  long 
after  the  pressure  waves  have  traveled  away.  The  center  of  this  hot 
region  of  gas  moves  away  from  the  cavity  at  the  speed  of  flow  of  the  gas 
through  the  cavity.  In  measurements  made  very  near  the  laser  cavity, 
the  hot  region  of  gas  was  found  to  have  the  general  shape  as  sketched  in 
Figure  18. 

Figure  19  shows  the  actual  measured  profiles  for  the  15  cm  cavity 
for  the  N^,  Ar,  and  He  gases.  The  probe  was  located  at  location  # 1 
which  is  143.5  cm  away  from  the  center  line  of  the  discharge  region. 

Figure  20  shows  the  profile  for  the  45  cm  cavity  using  Ar  gas.  Figure  21 
shows  the  measured  temperature  profiles  for  the  100  cm  case  for  and 
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Shot  # 3155 
Date  1/26/77 

Vertical  Scale  .02v/d 
Horizontal  Scale  .5  msec/d 

Notes 

N2 

Fan  On 


Shot  = 3154 
Date  1/25/77 

Vertical  Scale  .02v/d 
Horizontal  Scale  .5  msec/d 
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Arqon 
Fan  On 


Shot  - 31 70 
Date  1/27/77 

Vertical  Scale  .02v/d 
Horizontal  Scale  1 msec/d 

Notes 

He 

Fan  On 


FIGURE  14.  EXAMPLES  FOR  PRESSURE  WAVE  FORM  FOR  THE  15  CM  CAVITY  WITH 
N„.  Ar.  AND  He. 


:CirV 


irrC  1*C 


Shot  # 3227 
Date  2/2/77 

Vertical  Scale  .01v/d 
Horizontal  Scale  1 msec/d 

Notes 

Argon  (45  cm) 

Fan  On 
Corner  Test 


FIGURE  15.  EXAMPLE  OF  PRESSURE  WAVE  FORM  GENERATED  BY  45  CM  CAVITY  USING 
Ar  GAS. 
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FIGURE  16.  EXAMPLES  OF  PRESSURE  WAVE  FORM  GENERATED  BY  100  CM  CAVITY 
USING  N0  AND  Ar  GASES. 
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Vertical  Scale  .02v/d 
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Shot  # 3143 

Vertical  Scale  ,02v/d 
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Notes 
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FIGURE  17. 


EXAMPLE  OF  PRESSURE  WAVE  WHEN  FAN  IS  OFF.  THIS  IS  FOR  15  CM 
CAVITY  AND  He  GAS. 
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Temperature , 


FIGURE  18.  THERMAL  WAVE  FORM  CHARACTERISTICS. 
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FIGURE  19.  EXAMPLES  OF  THERMAL  WAVE  FORM  FOR  THE  15  CM  CAVITY  WITH  Ar, 
N„.  AND  He  GAS. 
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FIGURE  20.  EXAMPLE  OF  THERMAL  WAVE  (UPPER  TRACE)  GENERATED  BV  45  CM 
CAVITY  USING  Ar  GAS. 
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Shot  # 3180 
Date  2/1/77 

Vertical  Scale 
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FIGURE  21.  EXAMPLES  OF  THERMAL  WAVE  (LOWER  TRACE  IN  BOTH  EXAMPLES) 
GENERATED  BY  100  CM  CAVITY  IN  N?  AND  Ar  GAS. 


2.  Wave  Interaction  with  Blower 


The  blower  is  the  first  system  component  the  waves  encounter 
upon  leaving  the  downstream  side  of  the  cavity.  The  experimental  results 
show  that  much  of  the  pressure  wave  is  reflected  back  toward  the  cavity. 
The  thermal  wave  passes  through  the  fan  without  reflection  but  is  re- 
duced in  amplitude  and  spread  over  a larger  distance. 

The  general  characteristic  of  the  pressure  wave  interaction  with 
the  fan  is  illustrated  in  Figure  22.  All  these  features  can  be  seen 
on  a typical  trace  shown  in  Figure  23  where  the  upper  trace  is  from  a 
probe  on  the  upstream  side  of  the  fan  and  the  second  trace  is  from  a 
probe  on  the  downstream  side. 

The  thermal  wave  interaction  with  the  blower  is  illustrated  by  the 
data  shown  in  Figure  24  for  shot  numbers  3217  and  3219.  For  shot  number 
3217,  the  temperature  probe  is  on  the  upstream  side  of  the  blower  while 
for  shot  number  3219  the  probe  is  on  the  downstream  side  of  the  blower. 

The  power  levels  for  the  two  shots  are  322J  and  289J  respectively. 

3.  Wave  Attenuation  with  Respect  to  Distance 

Figure  24  illustrates  the  effect  of  distance  of  propagation 
on  the  thermal  wave  in  the  record  for  shot  number  3217.  The  time  sweep 
was  set  at  .1  s/div  so  that  the  second  passage  of  the  wave  could  be 
observed  at  this  same  location.  The  experiment  shows  that  the  thermal 
wave  took  about  .5  sec  to  travel  around  the  system  and  was  reduced  in 
magnitude  by  a factor  of  about  8 after  a single  pass. 

The  attenuation  of  the  pressure  wave  as  it  travels  around  the  open 

system  is  shown  in  Figure  25  for  shot  number  3079  which  shows  the  pressure 

probe  records  at  the  four  major  locations  in  the  system.  The  effect  of 

distance  is  seen  by  following  the  wave  that  travels  upstream  from  the 
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FIGURE  22.  CHARACTERISTICS  OF  THE  PRESSURE  WAVE  INTERACTION  WITH  THE 
BLOWER. 
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FIGURE  23.  EXAMPLE  OF  PRESSURE  WAVE  INTERACTION  WITH  FAN  (15  CM  CAVITY 
ARGON  GAS). 
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FIGURE  24. 


EXAMPLES  OF  INTERACTION  OF  THERMAL  WAVE  WITH  THE  FAN.  LOWER 
TRACE  OF  *321 7 IS  TEMPERATURE  PROBE  UPSTREAM  OF  FAN  WHILE 
UPPER  TRACE  OF  3219  IS  TEMPERATURE  PROBE  DOWNSTREAM  OF  FAN. 
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FIGURE  25.  EXAMPLE  OF  PRESSURE  WAVE  PROPAGATION  IN  OPEN  SYSTEM  WITH  THE 


cavity  beginning  at  station  4 then  appearing  at  station  3 and  finally 
at  station  2.  Taking  the  amplitude  at  station  4 as  a reference,  the 
positive  phase  amplitude  at  3 and  2 is  .718  and  .684  respectively.  These 
results  are  with  the  fan  on.  For  the  case  when  the  fan  is  off,  refer 
to  Figure  26  which  shows  the  results  for  shot  number  3083.  In  this  case, 
the  amplitude  at  3 and  2 relative  to  4 is  .814  and  .804  respectively. 

This  comparison  shows  that  the  pressure  wave  attenuation  is  much  greater 
when  the  wave  is  traveling  against  the  flow  as  compared  to  when  there 
is  no  flow. 

4.  Pressure  Wave  Interaction  with  Ceramic  Reflectors 

The  ceramic  reflectors  were  installed  between  stations  3 and 
4 in  the  recirculation  system.  The  thin  walled  metal  honeycomb  was  found 
to  have  no  measurable  effect  on  the  wave  propagation.  The  pressure  wave 
was  found  to  be  attenuated  by  the  square  celled  ceramic.  They  were  not 
found  to  cause  any  measurable  reflection  of  the  pressure  wave,  however. 
Examples  of  the  results  are  shown  in  Figure  27  for  shot  number  3097  with 
the  fan  on  and  number  3098  with  the  fan  off.  For  these  shots,  all  four 
sections  (8  cm  in  total  length)  of  the  square  celled  ceramic  were  placed 
between  stations  3 and  4.  Using  the  amplitude  at  station  4 as  the  re- 
ference, the  amplitude  ratio  across  the  ceramic  was  .562  with  the  fan  on 
and  .632  with  the  fan  off.  Compared  with  the  same  ratios  discussed  in 
part  3 above,  the  ceramic  is  found  to  reduce  the  amplitude  of  the  pressure 
wave  by  a factor  of  .78  compared  to  the  open  pipe. 

5.  Pressure  Wave  Interaction  with  Muffler 

The  muffler  was  found  to  have  a dramatic  effect  on  the  pressure 
wave  as  can  be  seen  in  Figure  28  which  shows  the  results  for  shot  numbers 

3102  and  3106  for  fan  and  no  fan  respectively.  The  pressure  wave 
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Date  12/14/76 
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Horizontal  Scale  10  msec/d 
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FIGURE  27.  EXAMPLES  OF  PRESSURE  WAVE  PROPAGATION  THROUGH  CERAMIC  RE- 
FLECTORS WITH  FAN  ON  (*3097)  AND  FAN  OFF  (#3098). 
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Shot  # 3102 
Date  12/17/76 

Vertical  Scale  .02  v/d 
Horizontal  Scale  5 msec/d 

Notes 

Argon 

Muffler 

Fan 

90 

21.8 

Strong 


Shot  # 3106 
Date  12/17/76 

Vertical  Scale  .02v/d 
Horizontal  Scale  5msec/d 

Notes 
Argon 
Muffler 
28.9/90 
No  Fan 
Long  Pulse 


FIGURE  28.  EXAMPLES  OF  PRESSURE  WAVE  PROPAGATION  THROUGH  THE  MUFFLER 
WITH  THE  FAN  ON  (#3102)  AND  FAN  OFF  (#3106). 
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amplitude  ratio  across  the  muffler  is  found  to  be  .143  with  the  fan  on 
(#3102)  and  .181  with  the  fan  off  (#3106). 


Pressure  probes  were  relocated  nearer  the  muffler  and  a more  de- 
tailed study  was  made  of  the  pressure  wave  interaction  with  the  muffler. 

This  study  showed  that  there  is  a negative  phase  pressure  wave  reflected 
from  the  muffler  when  the  pressure  wave  travels  into  the  muffler.  The 
positive  phase  of  the  pressure  wave  is  attenuated  considerably  in  the 
muffler  while  the  negative  phase  is  not  reduced  as  much.  A sketch  of 
the  pressure  wave  Interaction  is  given  in  Figure  29.  Figure  30  shows  data 
which  supports  the  general  view  of  the  muffler  interaction  as  shown  in 
Figure  29.  In  Figure  30,  data  from  shot  numbers  3172  and  3143  are  shown. 

For  #3172,  the  two  pressure  probes  are  placed  at  the  same  side  of  the 
muffler  spaced  about  22  cm  apart.  In  the  lower  trace  of  #3172,  the  re- 
flected negative  wave  is  on  top  of  the  negative  phase  of  the  incident 
wave.  In  the  top  trace,  the  reflected  negative  wave  is  seen  separated 
from  the  incident  pressure  wave.  The  data  on  shot  #3143  was  taken  with 
one  pressure  probe  on  each  side  of  the  muffler.  The  upper  trace  shows 
the  transmitted  wave  and  the  lower  trace  shows  the  incident  wave.  The 
second  negative  phase  wave  on  the  lower  trace  of  #3143  is  the  reflected 
negative  phase  wave  as  has  been  discussed. 

The  performance  of  the  muffler  was  investigated  in  more  detail  with  the 
assistance  of  Mr.  Cason  and  Dr.  Werkheiser  of  MICOM.  They  provided  a 
pressure  spectrum  of  the  pulse  which  is  typical  of  that  generated  in  a 
laser  pulse  shown  in  Figure  30a  . This  figure  shows  that  much  of  the  energy 
In  the  pressure  pulse  is  contained  in  the  100  to  1000  Hz  range.  The  pressure 
wave  data  taken  in  laser  firings  Number  3134,  3137,  3139,  3140,  3142  and 
3143  was  photographically  enlarged  and  over  100  data  points  taken  from 
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FIGURE  29.  CHARACTERISTICS  OF  PRESSURE  WAVE  INTERACTION  WITH  MUFFLER. 
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Shot  # 31 72 
Date  1/28/77 

Vertical  Scale  .02  v/d 
Horizontal  Scale  1 msec/d 

Notes 

Both  delayed  8 msec 
N~,  Fan  On 
MOff  Test 


Shot  « 3143 
Date  12/22/76 

Vertical  Scale  .02  v/d 
Horizontal  Scale  .5  msec/d 

Notes 

He 

Fan  Off 


FIGURE  30. 


EXAMPLE  OF  PRESSURE  WAVE  INTERACTION  WITH  MUFFLER;  #3172 
SHOWS  THE  INCIDENT  AND  REFLECTED  WAVE,  #3143  SHOWS  THE 
TRANSMITTED  WAVE. 
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FIGURE 


each  trace.  These  data  were  input  into  a computer  program  to  obtain  the 
power  spectrum  density  of  the  incident  and  transmitted  waves.  From  this 
data,  the  attenuation  characteristics  of  the  muffler  was  determined  as  a 
function  of  frequency.  The  results  are  shown  in  Figure  30b  for  three 
different  gas  mixtures. 

The  results  show  that  nearly  40  db  attenuation  was  obtained  at  the 
low  frequency  end  and  a 12  db  peak  at  near  500  Hz.  The  first  peak  is 
likely  due  to  the  large  volume  available  for  resonance  and  the  5u0  Hz 
peak  is  likely  resulting  from  the  horn  characteristics.  The  results  show 
a need  to  design  for  the  100  to  200  Hz  frequency  range  since  this  muffler 
had  little  effect  on  energy  in  that  range.  The  range  above  1000  Hz  can 
likely  be  increased  by  lining  the  duct  or  other  standard  approaches  for 
high  frequency  energy. 

6.  Pressure  Wave  Interaction  with  90°  Turns 

Data  was  taken  on  the  pressure  wave  at  a location  before 
and  after  it  passed  a 90°  turn.  No  measurable  change  in  the  wave  form 
could  be  observed  and  no  measurable  reflection  could  be  observed.  Samples 
of  the  data  are  presented  in  Figure  31  which  shows  data  for  shot  numbers 
3211  and  3238  for  the  100  and  45  cm  cavity  cases  with  the  fan  on.  The 
upper  traces  show  the  wave  form  before  the  corner  and  the  lower  traces 
show  the  wave  form  after  the  corner. 
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Shot  # 3211 
Date  2/2/77 

Vertical  Scale  .01  v/d 
Horizontal  Scale  1.  msec/d 

Notes 

Argon 

Fan  On  (100  cm) 

Corner  Test 


Shot  * 3238 
Date  2/2/77 

Vertical  Scale  .01  v/d 
Horizontal  Scale  .5  msec/d 

Notes 

Argon  (45  cm) 

Fan  On 
Corner  i oC 


F 1 0 RE  31.  EXAMPLES  OF  PRESSURE  WAVE  INTERACTION  WITH  CORNER.  TOP  TRACES 
SHOW  WAVE  BEFORE  CORNER  AND  LOWER  TRACE  SHOWS  SAME  WAVE  AFTER 
CORNER. 


76 


7.  Velocity  Wave  Form 


A limited  amount  of  data  was  obtained  from  the  hot  wire  on 
velocity  waves  in  the  system.  An  example  of  the  results  obtained  on  the 
velocity  wave  is  shown  in  Figure  32  for  shot  number  3125.  The  third 
trace  shows  the  velocity  as  measured  at  station  1.  The  upper  trace  is 
the  pressure  record  at  that  same  location.  The  results  show  the  sudden 
change  in  velocity  caused  by  the  shock,  the  decay  in  velocity  caused  by 
the  expansion  wave  following  the  shock,  the  sudden  velocity  increase 
caused  by  the  shock  traveling  back  past  station  1 after  reflecting  off 
the  fan,  and  the  decay  in  velocity  due  to  the  traveling  expansion  fan. 


Chapter  IV 

DESIGN  OF  LASER  CAVITY 


The  UAH  team  assisted  Mr.  H.  C.  Tom  in  the  engineering  design  of 
the  laser  cavity  to  be  incorporated  in  the  laser  system  under  design  and 
construction.  Periodical  meetings  were  held  to  provide  technical  in- 
formation and  design  guidance  in  the  detail  flow  calculations,  thermo 
analysis,  and  configuration  design  of  the  laser  cavity.  Some  laboratory 
analyses  using  a water  table  was  employed  during  this  period  in  obtain- 
ing pertinent  technical  data  for  the  design. 
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Chapter  V 
CONCLUSIONS 

The  theoretical  and  experimental  program  which  was  begun  during  the 
duration  of  this  contract  is  not  complete  at  this  time  but  the  following 
conclusions  have  been  drawn  based  on  the  subscale  testing  and  the  numerical 
model ing  work. 

1.  Impulsive  high  energy  deposition  in  a flowing  gas,  characteristic 
of  a high  energy  laser  pulse,  produces  shock  waves  which  travel  at  super- 
sonic speeds,  both  upstream  and  downstream,  away  from  the  cavity.  The 
energy  pulse  also  produces  a region  of  hot  gas  which  travels  with  the 
flowing  gas. 

2.  The  shock  waves  travel  throughout  the  recirculation  system, 
travel  around  90°  turns  with  little  or  no  attenuation,  reflect  off  a 
compressor  with  little  or  no  attenuation,  and  can  likely  be  absorbed  by 
mufflers  of  suitable  design. 

3.  The  hot  gas  region  (or  thermal  wave)moves  at  the  flow  velocity,  is 
attenuated  by  the  compressor,  and  will  return  to  the  cavity  in  one  re- 
circulation cycle  unless  suitable  heat  transfer  takes  place  in  the  circuit. 

4.  One  dimensional  modeling  of  the  transient  flow  behavior  is  able 
to  predict  the  major  qualitative  characteristics  of  pulsed  operation 
(i.e.,  the  shock  wave  and  thermal  wave  generation)  but  quantitative  com- 
parisons show  a need  for  improved  gas  dynamic  modeling  of  the  transient 
behavior. 

5.  Acoustic  absorbers  of  the  honeycomb  type  tested  in  this  work 

were  found  to  be  of  little  influence  in  the  attenuation  or  reflection  of 

the  shock  wave  nr  the  thermal  wave.  Either  a much  more  closed  honeycomb 
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or  many  more  of  the  ones  employed  would  be  required  to  attenuate  the 
shock  wave.  This  would  be  at  the  expense  of  significantly  higher  pressure 
drop. 

6.  A low-pressure-drop  muffler  can  likely  be  designed  to  eliminate 
che  shock  waves  produced  by  the  pulsed  laser  operation.  The  muffler 
tested  in  this  work  was  successful  in  attenuating  the  shock  wave  by  40  db 
and  better  performance  can  likely  be  obtained  by  redesign  and  further 
testing. 

7.  The  numerical  simulation  of  steady  state  operation  and  transient 
fluid  and  thermal  behavior  needs  further  development  which  could  come 
about  from  imperial  relations  based  on  actual  system  performance.  The 
subscale  system  employed  in  the  work  reported  here  was  not  appropriate 
for  development  of  the  data  needed  on  start  up  transients  and  heat  ex- 
changer performance. 

The  work  reported  here  has  revealed  the  major  transient  fluid  and 

thermal  characteristic  of  high  energy  gas  laser  systems  operated  in  a 

pulsed  mode.  While  the  closed  cycle  EDL  has  been  the  basis  for  the 

study,  the  results  should  be  of  interest  and  application  to  the  design  of 

both  closed  and  open  cycle  systems  including  electric  discharge,  chemical, 

and  gas  dynamic  high  energy  systems.  The  closed  cycly  system  presents  the 

most  severe  environment  in  terms  of  pulsed  operation  since  waves  can  travel 

around  the  system  or  reflect  back  to  the  cavity  and  thus  spoil  beam  quality 

or  otherwise  degrade  the  gas  properties  desired  for  efficient  laser 

operation.  Wave  attenuation  is  required  in  order  to  prevent  waves  from 

returning  to  the  cavity  region  or  from  traveling  into  the  upstream  regions 

of  the  system.  Such  attenuation  can  most  likely  be  performed  using  high 

pressure  drop  screens,  honeycomb  material,  or  orifice  plates.  The  overall 
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efficiency  of  the  laser  system  is,  however,  reduced  by  the  use  of  high 
pressure  drop  devices  and  a better  solution  would  be  to  employ  low  pressure 
drop  muffler  concepts.  The  muffler  concept  employed  in  the  work  re- 
ported here  was  of  the  straight  flow  type  having  no  structure  to  impede 
the  flow  of  laser  gas.  The  waves  are  attenuated  by  designing  resonate 
cavities  and  horn  amplifiers  to  transmit  the  A.  C.  energy  in  the  flow  to  be 
disipated  in  acoustic  absorbing  material  while  not  effecting  the  D.  C. 
energy  contained  in  the  flow.  Improvements  on  this  concept  and  design 
are  confidently  predicted  and  such  devices  should  find  wide  application 
in  pulsed  laser  operation. 
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A DH(3),  DW ( 3 ) , DDHGC3),  DALGC3),  DAF0AC3),  DSIGC3),  DDEl(3), 

B 0 A ( 3 ) , DK(3>,  DSLC3),  DSIC3),  D A L F ( 3 ) , D D H ( 3 ) , C0NV(16),  DELM, 

C R 0 L I M , PLIM,  FG,  TITLF(28),  DP,  DPI,  DPII,  PPR,  PPRM, 

D T A 1 F , RHP,  TIGM,  DKH(3),  ETIG,  E T 02  G , ET06G,  E T07G  , E T 1 , 

E E T 2 , F T 3 , ETA,  E T ft , E T 8 , F T 1 0 , E T 1 A , FT16,  ET18,  ET20,  E T A 2 , 

F EEPR,  C W A , DNP(3)  , GDI,  QDII,  ODA,  EVIM,  PRPM,  VLIM,  R0A1, 

G R0A2.DAFAL  (3)  ,DSLL(3)  , AP(3) 

COMMON  / I V A R / N T AM , NDSGN,  NCS,  KI,  KO,  NLIM,  NP(3),  NFULL, 

1 NeS,  NBPS,  JCPRO 
EQUIVALENCE  C S V P C 1 ) ,T01G) 

2  CALL  ENPUT 
T2G=T02G 
T6G=T06G 
T5G=T02G 
T10G=TC7G 
J V = D 
SOI  11=0 

CALL  GASPCTIG,  FN,  GAMIG,  CPIG,  Z 1 , Z 1 , GMU) 

IF  (JCPRO  .GE.  0) 

1 PIG=R01 *T I G/ ( .0931 78*GMU ) 

IF  (JCPRO  .LT.  0) 

1 R01=.093178*PIG*GMU/TIG 

wG  = P 0 1 *V1 *A  ( 1 ) 

FM1=WG*DSQRT(TIG/(GAMIG*GMU))/(20.777A*PIG*A(1)) 

T0lG=TIG*(1.+FMl*FMl*(GAMIG-1.)/2.) 

T1 1G=T01G 

WRITE(o,19fc)R01,oMU,PIG,TIG 

198  F0PMAT(lH0,3x,'ROl=',El5.7,Sxf'GMU=',ElR.7,SX,'°IG=',El5.7, 

1 f X , 'TIC  = ",F15.7) 

10  DO  3 1 = 1 ,16 
DO  A J = 1 , 1 c 
A AMT ( I , J ) =0  . 
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r : 1 

1 3 6MT(I)=0.0 

WRITE(6,122)  T 0 2 G 
122  FORMAT(1hO,3X,'T02G=',El5.7) 

URITE(t,121)  T1,T2,T3,T4,T6,T8,T10,T14,T16,T18,T20 
121  FORWAT(1h0,5X,'T1,T2,T3,T4,T6,T8,T1P,T14,T16,T1?,T20',/,3X, 

1 9E13.6,/,3X,2E13.6) 

TMI=.5*(T2*T4) 

TMII=.5*(T2*T6) 

TM= .5* (T3*T1  ) 

Z1 =T5G-T4 
Z 2 = TfcG-T2 

IF  (Z1  .LE.  0.  .OR.  12  .LE.  0.)  GO  TO  11 

IF  <Z1  .EQ.  Z 2 ) GO  TO  12 
TGMI=TMIf(Zl-Z2)/(DLOG(Zl)-DLOG(Z2)> 

GO  TO  13 

1 1 TGMI  = . 5*<T5G+T6G) 

GO  TO  13 

1 2 TGWI  = .5* (T5G  + T6G) 

13  Z1=T10G-T6 
Z 2 = T 1 1G-T2 

IF  C Z 1 .LE.  0.  .OR.  Z 2 .LE.  C.)  GO  TO  14 

IF  (Z1  .EQ.  Z 2 ) GO  TO  15 
TGMII=TF"II  + (Z1-Z2)/(OLOG(Z1)-DLOG(Z?)) 

GO  TO  16 

14  TGHII=.5*(T10G*T11G) 

GO  TO  16 

15  TGMII=.5*(T10G+T11G) 

16  Z 1 = T 1 -T  A 1 
Z2=T3-T A2 

IF  <Z1  .LE.  0.  .OR.  12  .LE.  G.)  GO  TO  17 

IF  (Z1  .EQ.  Z 2 ) GO  TO  18 
TMA=TM-(Z1-Z2)/(DL0G(Z1>-DL0G(Z2)) 

GO  TO  19 

17  TMA=.5*(TAl+TA2) 

GO  TO  19 

. 1 £ T* A= . 5* ( T A 1 +T A 2 ) 

19  CALL  GASPCTGMI,  FN,  Z 1 * CPGI,  PRGI,  GIMU,  GMU ) 

CALL  GASPCTGMII,  FN,  Zl,  CPGII,  PRGII,  GIIMU,  GWU) 

WRITE(6,199)GIMU, GIIMU, VI 

199  FORMAT(lHO,3X,'GlWU=',El5.7,5X,'&II"U=',El5.7,5X,'VI=',E15.'n 
CALL  AIPPCTA1,  PAT,  TMA,  RH,  PRA,  AWU  , CPA,  Zl,  AWL) 

CALl  LIQP(TMI  , FG,  ZIMU,  PRI,  CPI) 

CALL  LIQPCTMII,  FG,  ZIIWU,  PRII,  CPII) 

CALL  LIQPCTM,  FG,  ZMU,  PR,  CP) 

AFRGI=DL ( 1 ) *DH ( 1 ) 

AFRGII=DL(2)*DH(2) 

A F R A = DL ( 3 ) *DH  ( 3 ) 

GG  I =WG / ( DS1G ( 1 ) * AFRCI  ) 

GGII=WG/(DSIG(2)*AFRGII) 

GA=WA/  (DSIG (3) *AFR A ) 

AFRI=0W( 1 ) *DH< 1 ) /ONP ( 1 ) 

AFRII=DW(2)*DH(?)/DNP(2) 

AFR  = 0W(3)*DH(?) / 0 N P ( 3 ) 

G I sWI / (OSI  ( 1 ) *A  FRI  ) 

GI  I=WI I / (DSI (2) *AFR  1 1 ) 

W 3 * W I ♦ W I I+Wfl+WM+WVl  ♦ l»V2*WE+WS 
G*W3/  <OSI (3)*AFR) 
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REGI=GGI *DDHG  (1  )/GI*U 

REG1I=GGII*DDHG(2>/GIIMU 

REA=GA*0DHG(3)/AMU 

RE  I =GI *DDH  < 1 ) /Z I MU 

RE  I I = G I I *DDH ( 2) / Z I I VU 

RE=G*DDH(3) /Z*U 

taRITE(6,197)REGI,REGlI,REA,REI,REII,RE 
197  F0RMAT(lH0,3X,''PEGI  = ',El3.6,4Xf'REGII  = ',El3.6,4X,'REA  = 'fE13.6f4X, 

1 'REI=',El3.6,4Xf'REII=',E13.6,4x,'RE=',El3.6> 

Z1=DL0G(REGI ) 

Z2  = -75G.2946+62<>.202fe*Zl-219.5279*Zl**2  + 40.7569*Zl**3-4.253034*Z1* 
1*4  + ,23e>4432*Zl  **5-.005466298*Z1  * *6 
GFI=DEXP(Z2> 

Zl=-360. 4224+335. 276*Zl-129.1436*Z1**2+26. 2 04 69*Z1**3-2.961196*Z1* 
1*4+.1766745*Z1**5-.004345743*Z1**6 
6J  I=0EXP(Z 1 ) 

Z1=DL0G(RE6II ) 

Z2  = -75G. 2946+629. 2028*Zl-219. 52  79 *Zl**2  + 40.7569*Zl**3-4. 253 03 4*Z1* 
1*4+.2364432*Z 1**5-. 005466298*21**6 
GFII=oEXP(Z2) 

Zl  = -360. 4224  + 335. 276*21-129. 143fc*Z1**2+26. 2 04 69*Zl**3-2. 961196*21* 
1*4+.1766745*Z1**5-. 00 4345743*21**6 
G J I I = DE  XP ( Z 1 ) 

Z1=DL0G(REA) 

Z2=-7. 652462+6. 206 772*Zl-1.7151?7*Zl**2+. 17561 96 *Z1**7-. 006 2 07 8 06* 
1 Z 1 * *4 
AF=DEXP(Z2) 

Z 1 = 2 . 123751-1  . 08 3736* Z 1 - . 05 2 05 2 97  * Z 1 * *2 + . 01  64  66  61  *Z1  **3-. 00078  2 2 0? 
1 1 *Z 1 * *4 
AJ=DEXP(Z1 ) 

Zl-DLOG(REI) 

Z1=-22.04452+12.8495b*Zl-3.122052*Zl**2+.3159152*Zl**3-.01161237*Z 

11**4 

ZJI  = DEXP(Z1  ) 

Z 1 =DLOG (RE  1 1 ) 

Zl=-22. 04452+1 2. 84988*Z1-7.122052*Z1**2+. 31 59152*Z1+*7-. 01161237*2 
11**4 

Z J I I = OEXP( Z1  ) 

Z1=D10GCRE) 

7l=-22. 04452+12. 84 98S*Zl-3.122052*Z1**2+. 31 59 152*71**3-. 01 161237*2 
11**4 

Z J=DEXP(Zl  ) 

T03=-2./3. 

HGI=GGI*CPGI*GJI*PRGI**T03 

HGII=GGII*CPGII*GJII*PRGII**T03 

HA=GA*CPA*AJ*PRA**T03 

HI=GI*CPI*ZJI*PPI**T07 

HII=GII*CPII*7JII*PPII**T03 

h=g*cp*zj*pr**to3 

GMI=DSQRT(2.*HGI/(DK(1)*DDEL(1)>) 

GMII=DSQRT(2.*HGII/(DKC2)*DD£L(2))) 

AM=CSQRT  (2  . * H A / (OK  (7)*DDEL  (3)  ) ) 

T03  = GWI*DSL  (1  ) 

A T F I =D  T ANH ( T07 ) /T0  3 
T03=G*I I *DSL ( 2 ) 

ATF I I=DTANH (T03) /T03 
T03=AM*DSL (3) 
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ATFA=0TANH(TO3) /T03 
AT0I=1.-0AF0A(1)*(1.-ATFI) 

A TO  1 1 = 1 ,-DAFOA (2) * ( 1 AT F 1 1 ) 

AT0=1.-DAF0A(,)*(1.-ATFA) 

GMl=DS<dRT(2.*HI/(DK(1)*DA(1))) 

GM2=DSQRT(2.*HII/(DK(2)*DA(2))) 

AMN=DSGRT(2.*H/(DK(3)*DA(3))) 

T03=G*1 *DSLL ( 1 ) 

ATF 1 L=DTANH (T03) /T03 
T03  = GI*2*DSLL(2) 

ATF2L=bTANH (T03) / T 0 3 
T03=AMN*DSLL ( 3 ) 

ATFAL=DTANH(T03)/T03 

AT01  = 1 .-0AFALC1 ) * ( 1 .-ATF1L) 

ATO2=1.-0AFALC2)*C1.-ATF2L) 

ATOL  = 1 ,-DAFAL(3)*(1  .-ATFAL) 

UI=1./(1./(ATOI*HGI)  ♦ AF>(1)/(DK(1)*(1.-DAFOA(1)))  ♦OALG(1)/(DALF(1 
1 ) *HI  *AT01 ) ) 

UII=1./(1./(ATOII*HGII)«-AR(2)/(DK(2)*(1.-DAFOA(2)))*DALG(2)/(DALF( 
12)*HI*AT02) ) 

U=1./(1./(ATO*HA)+AB(3)/(DK(3)*(1.-DAFOA(3)))+DALG(3)/(DALF(3)*HI* 
1 ATOL  ) ) 

CGI=WG*CPGI 
CGI  I=WG*CPGI I 
C A=WA*CPA 
C I=WI*CPI 
C I I=WI I *CPI I 
C =W  3 * C P 

CKINI=OMIN1 (CGI  ,CI  ) 

CMAXI=[,MAX1 (CGI ,CI) 

CMINI I =DMINl (CGI  I ,C I I) 

CWAXI I=DMAXl  ( C G I I * C I I ) 

CMIN  = DMIN1 ( C A , C ) 

C MA  X = Df*A  X 1 ( C A , C ) 

AGI=DL(1 ) * D H ( 1 ) * D W ( 1 )*DALG(1 ) 

AGII=DL(2)*OH(2)*DW(2)*DALG(2) 

AA=DL(3)*DH(3)*DW(3)*DALG(3) 

TUM  = A G I *U  I / C ^ I N I 

TUNII*AGII*UII/CMINII 

TUN=AA*U/CMIN 

TOWI=1.-DEXP(-TUNI*CMIM/(CPAXI*DNP(1))) 

TOWII  = 1.-OEXP(-TUNII*CMIMI/(CMAXII*DNP(2))) 
TOW=1.-DEXP(-TUN*CMIN/(CMAX*DNP(3))) 

ESS  1 = 1 .-DEXP(-TOWI*CMAXI/CMINl> 

ESSII  = 1 .-0  EXP (-TOW  I I*CMAXII/C«INII) 

ESS=1.-OEXP(-TOW*CMAX/CWIN) 

IF  (CM1NI/CMAXI  .GT.  .9999)  GO  TO  31 

EPSI=((1.-ESSI*CMINI/CMAXI)/(1.-ESSI))**NP(1) 

EPSI«(EPSI-1.)/(EPSI-CMINI/CMAXI) 

GO  TO  32 

31  EPSI  = DNP(1 ) * E S S I / (1  •♦ESSI*(DNP(1 ) -1 •) ) 

32  IF  (CWiNII/CWAXII  .GT.  .9999)  GO  TO  33 
EPSII=((1.-ESSI!*CMINII/CMAX1I)/(1.-ESSII))**NP(2) 
fcPSIl  = (EPSII-1.)/(EPSII-CMIMI/CMAXII) 

GO  TO  34 

33  EPSII=DNP(2)*ESSII/(1.*ES$II*(DNP(2)-1.)) 

34  IF  (CMIN/CMAX  .GT.  .9999)  GO  TO  35 
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EPSA=((1.-ESS+CMIN/CMAX)/(1.-ESS))'**NP(3) 

EPSA=<EPSA-1.)/(EPSA-CMIN/CMAX> 

GO  TO  36 

35  EPSA=DNP(?)+ESS/(1.+ESS*(DNP(3)-1.)) 

3 6 CONTINUE 

T03=FMl*FMl 

P01G=PIG*(1.+T03*(GAM1G-1.)/2.)**(GAMIG/(GAMIG-1.>) 
GO  TO  (3G,  40,  50),  NDSGN 
: CAVITY  OPTION  1 CONST.  MACH  NUMBER 

30  F M 2 = f M 1 

A(2)=A(1)*(T02G/T0lG)*+((1.+GAMIG*T03>/2.) 

GO  TO  60 

C CAVITY  OPTION  2 CONST.  AREA 
40  A C 2 ) = A ( 1 ) 

Z1=TClG*(1.+GAMIG*T03)*+2/(2.*(GAMIG+1.)*T03+ 

1 ( 1 . + T03* (GAMIG-1  . ) /2  . ) ) 

21 =T02G/Zl 

IF  ( Z 1 .GT . 1 . ) GO  TO  10000 
1 = 0 

12000  1=1+1 

CALL  GASPCT2G  , FN,  GAM2G,  C P 2G , T03,  T03,  GMU ) 

12  = 2 . *GAM?G*Z 1 -2 . +GAM2G-2  . 
ZJ=GA*2G*GA*2G*(Z1-1.)+1. 

IF  ( Z 3 ,NE . 0.)  GO  TO  41 

FM2  = DS0RT(-Z1 / Z 2 ) 

GO  TO  60 


41  Z4=(-Z2+0SQRT(Z2*Z2-4.*Z1*Z3))/(2.*Z3) 
Z5=(-72-DSQRT(Z2*Z2-a.*71*Z3))/(2.*Z3) 

Z6  = DMIN1  ( Z 4 » Z 5 ) 

IF  ( Z 6 , L T . 0.)  Z6  = DMAXHZ4,  Z5) 

WRITE(6,174)T02G,T03,PR 

174  FORMAT(lH0,3X,'TO2G=',E15.7,5X,'TO3=',E15.7,5X,'PR=',E15.7) 
wRITE(6,173)Z1,72,Z3,Z4,Z5,Z6 

173  FORMAT(1HO,3X,'71=',E15.7,5X,'Z2=',E15.7,5X,'Z3=',E15.7,5X, 
1/,4X,'Z4=',E15.7,5X,'Z5=',E15.7,5X,'Z6=',F15.7) 

F M 2 = D S G R T ( Z 6 ) 

60  TO  60 

C CAVITY  OPTION  3 GENERAL  SHAPE 

50  Z5=(GAMIG+GAM2G)/2. 

Z 6 = F M 1 * F M 1 

N = NC  S-1 
DO  51  J = 1 ,N 

Zl=-2.+(1.+Z6*(Z5-1.)/2.)/(1.-Z6) 

Z 2= ( 1 . + Z5  + Z6)*(-.5*Z1) 

Z3=(ACN(J+1)~ACN(J))/(t.*(ACN(J+1)+ACNCj))) 

Z4= (T02G-T01G) / ( (NCS-1 . )+T0lG-( J-. 5)* (T02G-T01G) > 

Z6=Z6*C1 ,+Z3*Zl+Z2+74) 

IF  (16  .LF . 1 . ) 60  TU  51 
1 0000  WRI TE (KO, 1 05  ) 

105  FORMAT (4lH*+++  CHOKED  CAVITY  * READ  NEW  DATA  **«*♦*) 

GO  TO  2 

51  CONTINUE 

FM2  = DSQPT  (ZM 
A(2)=ACN(NCS) 

60  T260LD=T2G 

T2G=T02G/(1.+FM2*FM2*(GAM2G-1.)/2.) 

T2GD I F=T2G0LD-T2G 
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I F ( DABS ( T2GDI F ) .LE . . 5D 0)  GO  TO  11000 
I F ( I .LE .KLIM)  GO  TO  12000 
PRINT  1 3000  * I , T 2GD  I F 

13G0C  FORMAT (5fcH****FAILURE  TO  CONVERGE  IN  STATIC  TEMP.  LOOP  1 - CYCLES 
' = ,13,'  - ERR0R=',E12.5) 

GO  TO  2 

11000  P2G=WG*DSQRT(T2G/(GAm2G*GMU>)/(2G.7774*A(2)*FM2) 

P02G  = P2G*<1.*FM2*FM2*(GAM?G-1.)/2.)**<GAM!2G/(GAM!2G-1.>) 

: FIRST  DIFFUSER 

P03G=P02G-CF23*(P02G-P2G) 

FM3=FM2*A (2) / A (3) 

GAMi36  = &AM2G 


CALL  MACHfDFLM,  NLIM,  P03G, 
1 A(3),  P3G , T3G, FN.CP3G) 

G AM  3 G , 

FM3  , 

T 02  G , 

WG, 

GMU 

SECOND  DIFFUSER 

P04G=P03G-CF34* (P03G-P3G) 
FM4=FM3*A(3)  / A ( 4 ) 

GAM4G=G AM3G 

CALL  MACHCDELM,  NLIM,  P04G, 

G AM4G  , 

F M 4 , 

T 0 2 G , 

WG, 

GMU 

1 A ( 4 ) , P 4 G , T4G,FN,CP46) 
THIRD  DIFFUSER 

P05G=P04G-C  F 4 S * (P04G-P4G) 

FM5=FM4*A(4)/A(C) 

GAM5G=GAM46 

CALL  MACHCDFLM,  NLIM,  P05G, 

GAM5G , 

F M 5 , 

-1 

o 

ro 

C7> 

WG, 

GMU 

1 A(5),  P 5 G , T5G,FN,CP5G) 

R05=.093178*P5G*GMU/T5G 
R 06  = R 0 5 

FMfc  = WG*t>SQRT(T06G/GAM5G/GMU)/(20.7774*P056*A(6>) 

1=0 

70  Z 1 = RO  6 

P06G  = P05G-6GI  *GG  I * ( ( 1 . +D S I 6 ( 1 ) * D S I G ( 1 ) ) * (R05/R06-1 

1 GFI*DALG(1)*DW(1)*R05/(.5*DSIG(1)*(R06*R05)))/ 

2 C9266.11*R05) 

GAM6G=G AM5G 

CALL  MACH(DELM,NLIM,P06G,GAM6G,FM6,T06G,WG,GMU,A(6),P6G,T6G,FN, 

* CP6G) 

RO6=.093178*P6G*GMU/T66 

1*1*1 

Z3=R06-Z1 

IF  (GAE>S(Z3)  .LF.  ROLIM)  GO  TO  71 
IF  (I  .LT.  NLIM)  GO  TO  70 
WRITE (<0, 101  ) Z 3 

101  F ORMA  T ( 5 3H  * * * * FAILURE  TO  C0NVER6E  IN  DENSITY  LOOP  1 - ERROR  = , 
1 £12.5,6H*****0 

71  CALL  GASPCT07G,  FN,  G AM7G , CP7G,  T03,  T03,  GMU  ) 
wRITE(KO,150)  I , Z 3 

150  F ORMAT < 1 Ox , 'DFNS  IT Y LOOP  1 - CYCLES  =',I3,'  - ERROR  =',F12.5) 
P07G=P06G*BPR 
J =0 

80  FM7=WG*DSQRT(T07G/(GAM7G*GMU))/C20.7774*P07G*A(7)) 

CALL  maCHCDELM,  NLIM,  P 0 7 G , GAM7G,  F M 7 , T07G,  WG  , GMU  , 

1 A ( 7 ) , P7G,  T7G,FN,CP7G) 
fyferf M7*A (7)  / A ( p ) 

P08G  = P07G-C  F7R* (P0  7G-P7G) 

GAM8G=GAM70 

CALL  MACHCDFLM,  N L I m , P08G,  G AMSO  , F Ml  8 , T07G,  WG,  GMU  , 
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1 A ( £ ) , P8G,  T8G,FN,CP£G) 

P09G  = P08G-C F£9* ( P08G-P8G  ) j 

FM9=FMfc*A(8) / A(9) 

GAM9G=GAM8G 

CALL  "ACHCLELM,  NLIM,  P09G,  G A M 9 G , F M 9 , T 0 7 G , W G , GMU, 

1 A ( 9 ) , P 9 G » T9G,FN,CP9G) 

P01 0G  = PO9G-C  F 9 1 0 * (P09G-P9G) 

FM10=FM9*A(9)  / A(10) 

GAM1 0G=GAM9G 

CALL  «ACH(DELM,  NLIM,  P010G,  G A M 1 0 G , FMIO,  T07G,  WG,  GMU, 

1 A ( 1 0 ) , PIP G,  T10G,FN,CP10G) 

ROl0=.G9317R*P1CG*GMU/Tl0G 
ROI 1 = R G 1 C 

F M 1 1=WG*DSURT  CTO1 1 G/GAMlOG/GMU) / (20.7774*P01 OG*  A ( 1 0) ) 

1 = 0 

75  Z 1 = k 0 1 1 

P011G=P0106-GGII*GGI1*((1.*DSIG(2)*DSI6(2))*(R010/R011-1.)* 

1 GFII*DALG(2)*DW(2)*R010/(.5*DSIG(2)*(R010+R011)))/ 

2 <9266.11*R010) 

GAM1 1 G=GAM1 OG 

WRITE(6,183)P011G,GAM11G,FM11,T01G,WG 
183  FORMAT(1h0,3X,'PO11G=',E15.7,5x,'GAM11G=',E15.7»5X,*FM11s', 

1 E15.7-5X,'T01G=',E15.7,5X,'WG=',E15.7) 

CALL  MACH(DELM, NLIM, P0llG,GAMl1G,FMl1,T0lG,WG, GMU, A ( 1 0 ) , P 1 1 G , 

* T11G,FN,CPl1G) 

R011=.093178*P11G*GMu/T11G 
I-I  + 1 

Z3=R011-Z1 

IF  (0ALS(Z3)  .LE.  ROLIK)  GO  to  76 
IF  (I  .LT.  NLIM)  GO  TO  75 
WRITE(KO,120)  Z 3 

120  FORMAT (5 3H****  FAILURE  TO  CONVERGE  IN  DENSITY  LOOP  2 - ERROR  = , 

1 E12.5,6H******) 

76  FM11=WG*DSQ9T(T01G/(GAMIG*GMU))/(20.7774*P011G*A(11)) 

WP I Tt (KO , 1 5 1 ) I,  Z 3 

151  FORMAT (10X , 'DFNSITY  LOOP  2 - CYCLES  =',I3,'  - ERROR  =',E12.5) 

CALL  MACHCDELM,  NLIM,  P011G,GAM11G  ,FM11,  T01G,  WG,  GMU, 

1 A(11),  P 1 1 G , T1 1 G , FN , CPI 1G> 

IF  (JCPRO  .GE.  0) 

1 Z 1=R0l*TIfi/( .093178*GMU)*(1 .♦(GAMIG-1 .)*FM1*FM1  / 2.) 

2 ** (GAMIG / (GAMIG-1  . ) ) 

Z1=P011G-CF111*(P011G-P11G)-Z1 
IF  (JCPRO  .LT.  C) 

1 Z1=P011G-CF111*(P011G-P11G)-P01G 

P01 1G=P01 1G-Z 1 
P01 Ob=POl OG-Z 1 
P09G=P09G-Z 1 
P08G=PG8G-Z 1 
P07G  = P07G-Z  1 

CALL  GASP(T6G,FN,  GAM6G,  C P6G , Z 2 , Z 2 , GMU ) 

BPR=P07G/P06G 

Q6G=60.*A(6)*FM6*DSQRT(GAM6G*T6G*32.174*1545,43/GMU) 

ATAb=. 730-7. C9843* (DABS (&6G/BRPM/NBS-. 5355))**2.0P014 
A TAb  = 0 . 6 

WRITE (6,186)Z1 ,P07G,P06G,Q6G 

186  FORMAT(1HO,3X,'71=',E15.7,5X,'P07G=',E15.7,5X,'P06G=',E15.7,5X, 

1 '!J6G  = ',E1  5.7) 
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WRITECfc,187)  T06G , 9PR , 6AM6G , ATAP 

187  FORMATC1h0,3X,'TO6G=',E15.7,5x,'BPR=',E15.7,5X,'GAM6G=',E15.7,5x, 
I 'ATAB=',E15.7) 

TC7G  = T06G*  C 9 P R * * C CGAM6G-1  .)  / (AT  AB*GAM6G)  ) ) 

WPITE (6, 180)  T06G,BPR,GAM6G,ATAB,T07G 
180  F0RMATClH0,3X,'TO6G=',El5.7,5X,'BPR=',El5.7,5X,'GAM6G=',E15.7,5X, 
1 'ATAB=',E15.7,5X,'T07G=',E15.7) 

J = J + 1 

IF  (DABS(Zl)  .LT.  PLIM)  GO  TO  85 
IF  (J  .LT.  NLIM)  GO  TO  80 
WRITE  CKO, 1 03)  Z 1 

103  FORMAT C52H****  FAILURE  TO  CONVERGE  IN  PRESSURE  LOOP  - ERROR  = , 

1 E12.5,6H ******) 

85  CALL  AIRP(TA1,  PAT,  T A 1 , PH,  Z2,  Z2,  CPA1,  Z2,  AML) 

WRITE  U0 , 1 52)  J,  Z1 

152  FORMAT (1  OX , 'TOTAL  PASSES  THROUGH  PRESSURE  LOOP  =',I3, 

1 ' - tRROR  ='  ,E1 2 .5  ) 

CALL  AIRPCTAl,  PAT,  T A 2 , RH,  Zl,Zl,  C P A 2 , G A M 2 A , AML) 


CALL 

LIQPCT1 , FG, 

Z1  , 

ZI  , 

CPI  ) 

CALL 

L I 0 P C T 2 , F 0 , 

zi . 

Zi  , 

CP2  ) 

CALL 

L I Q P C T 3 , F G , 

Z1  , 

zi , 

C P 3 ) 

CALL 

L10PCT4 , FG , 

zi , 

21  , 

CPA) 

CALL 

L I Q P C T 6 , F G , 

ZI  , 

Zi  , 

CP6) 

CALL 

LIQPCTR,  FG, 

ZI  , 

zi  , 

C P 8 ) 

CALL 

LIQPCT10,  FG 

.Zi  . 

zi , 

CP10) 

CALL 

LIQPCT14,  FG 

.Zi  , 

zi , 

CPU) 

CALL 

LIQPCT16,  FG 

.zi , 

Zi  . 

CPI  6) 

CALL 

LIQPCT18,  FG 

,Z1  , 

zi , 

CP  18) 

CALL 

LIQPCT20,  FG 

tZi  , 

zi , 

C P 2 0 ) 

AMTC 

1,  1 ) =W6* CCPIG+CP 

26 ) 12. 

AMTC  1,  2)=-AMT(1,1) 

8MT ( 1 )=0S-.948*PLKW*( 1 ./EFFL-1 . ) 

AMT(2,2)*CGI 

AMTC2,3)=-AMT(2,2) 

AMTC2,t)=WI*.5*CCP2*CP4) 

AMT (2 ,t ) =-AMT (2 , 6) 

AMT(3,2)=EPSI*CMINI-CGI 
AMT(3,6)=-EPSI*CMINI 
AMTC3,3)=  CGI 
AMT C 4 , 4 ) = 1 . 

AMT  C4,3)=-CP07G/P06G) * * ( (GAMfcG-1 . ) / CATAB*GAM6G)  ) 

AMT  C 5 , 4)  = CGI  I 

AMT (5,1 ) =- AMT (5,4) 

AMT C5,6)=UII*.5*CCP2+CP6) 

AMT(5,9)»-AMT(5,6) 

AMT(6,4)=EPSII*CMINI1-CGII 
A M T ( A,  6)=-EPSII*CMlNlI 
AMT (6 , 1 ) =CGI  I 
AMT(7,7)«W3*.5*(CP1+CP3) 

AMTC7,5)=-AMTC7,7) 

AMTC7,16)*-CA 
BMTC7)=AMTC7,16)*TA1 
AMTC  B,  7)=FPSA*CMIN 
AMT(£,1A)=-CA 
BMT(fe)=(EPSA*CMIN-CA)*TAl 

AMT (9,5)*W3*.5*(CP1*CP2) 

AMTC9,6)*-A*TC9,5) 
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BMT (9 ) =-QP 

AMT  (10,8)=-WI*.5*(CP4+CP3) 

AWT (10,9)=-WII*.5*(CP6*CP3) 

AWT(in,13)=-WM*.5*(CPl6+CP3) 

A*',T(10,6)=-we*.5*(CP2  + CP3) 

AMT(10,1G)=-wVl*.5*(CPa*CP3) 

AMT(10,11)=-WV2*.5*(CP10+CP3) 

AMT(10,14)=-WE*.5*(CP1&*CP3> 

AMT(10,15)=-WS*.5*(CP2G*CP3) 

AWT(10,7)=-AMT(10,S)-AMT(10,9)-AMT(1C,13)-AMT(10,6)-AMT(10,1G)- 
* AMT(1G»11)-AMT(10,14)-AMT(1G,15) 

AMT(11,6)=wM*.5*(CP2+CP14) 

AMT ( 11,12 ) =-AMT (11,6) 

BMT (11 )=-QCM 

AMT(12,12)=WM*.5*(CP14+CP16) 

AMT(1?,13)=-AMT(12»12) 


PMT ( 1 2 ) =-QM 

AMT(1’,6)=wVl*.5*(CP2+CP8) 

AMT<13,10)=-A*T(13,<S) 

BMT (1 3) =-QVl 

AMT (14 ,6)=WV2*.5* (CP2+CP10) 

AMT ( 1 4 , 1 1 ) = - AM  T ( 14,6) 

6MT ( 1 4 ) =-0 V2 

AMT(15,6)=WE*.5*(CP2+CP18) 

AMT(15,14)=-AMT(15,6) 

BMT ( 1 5 ) =-QE 

AMT (16,6)=WS*.5*(CP2+CP2C) 

AMT (16,15)=-AMT  (16,6) 

BMT ( 1 6) =-QS 

CALL  SIMQCAMT,  PMT,  16,  kS) 

WRITE(6,122)  TO  26 
J J = 0 

6 Z 2 = F M 1 

TIG=bMTC1)/(1.+(6AMIG-1.)*FMl*FMl/2.) 

CALL  GASP (TIG,  FN,  G A M I G , CPIG,  Zl,  Z 1 , G MU ) 

IF  (JCPRO  . G E . 0) 

1 PIG=R01*TIG7 ( .0931 78*GMU) 

IF  (JCPRO  .LT.  0) 

1 RO1=.093178*PIG*GMU/TIG 

WG=ROl*VI*A(1  ) 

WRITE (6»160)TIG,GAMI6,6mu,wG 

160  FORMAT(1h0,3X,'T1G=',E14.6,5X,'GAMIG=',E14.6,5X,'GMU=',E14.6,Sx, 

1 'WG* ' , E 1 4 . 6 ) 

WRITE (6,162) 

162  FORMAT(1hO,8X,'SVR',16X,'BMT') 

WRITE (6 , 1 61 ) (SVP  (MP)  ,BMT  (MP)  ,MP  = 1 , 1 6) 

161  FORMAT < 3X , E 1 4 .6 , 5X , E 1 4 . 6) 
FM1=w6*DSQRT(Tir-/(GAMIG*GMU))/(20.7774*PIG*A(1)) 

J J = J J ♦ 1 

Z3=Z2-FM1 

IF  (0ABSCZ3)  .LT.  0 ELM ) GO  TO  7 
IF  (JJ  .LT.NLIM)  GO  TO  6 
WRITE (K0,4C0)  Z’ 

4 OG  FORMAT (SAH****  FAILURE  TO  CONVERGE  IN  FINAL  MACH  LOOP  - EPROR  = , 

1 El2.b,6H ******) 

7 WRITE  (KO  ,301  ) J J , Z 3 

301  FORMATdOX,  'TOTAL  PASSES  IN  FINAL  MACH  LOOP  =',I3,'  - ERROR  =', 
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w 


1 E12.5) 

11=11*1 

IF  (KS  . E Q . 1 ) GOTO  91 
DO  93  KS  = 1 , 16 
NN(KS)=-1 

TD1F(KS)=SVR(KS)-BMT(KS) 

IF  (DABS(TDIF(XS))  .GT.  CONV(KS))  NN(KS)=1 

93  CONTINUE 
WRITE ( K 0 , 201 ) 

201  FORMAT ( 10X , 'STATE  VARIABLE  ERRORS') 

WRITE(K0,300)  TCIF 
300  F0RMAT(1QX,10E1?.4) 

WRITE(K0,153)  II 

153  FORMAT  ( 1 OX , 'END  OF  SIMQ  PASS  NO. ',13/) 

IF  (II  .GT . NLIM)  GO  TO  92 
DO  94  KS  = 1 , 1 6 

IF  (NN(KS)  .GF.  0)  GO  TO  95 

94  CONTINUE 
GO  TO  90 

92  WRITE(KO,104) 

104  FORMAT (48H****  FAILURE  TO  CONVERGE  IN  STATE  VARIABLE  LOOP*, 

1 6h******) 

90  FM6=WG*DSQRT(TOOG/(GAM6G*GMU))/(20.7774*PO6G*A(6)) 

CALL  MACH(DELM,  NLIM,  P06G,  GAM6G,  FM6,  T06C,  WG,  G M U * 

1 A ( 6 ) , P6G,  T6G, FN,CP6G) 

FNBPS  = DFLOAT  (N6PS  ) 

FNBS=DFLOAT(NSS) 

DPSS=((P07G-P06G)*T06G*28.96t>)/(P06G*PRPM*BPPM*516.7*GMU*FNPPS) 
DPSb=. J396-7. 2323* ( DABS (Q6G/CBR PM*  FNbS) -.46*75  ) ) **1.42327 
DPSb=0PSe*1 . E -9 
E RR  = DPSB-DPSS 

GC=11.661E“7/(BRPM*FNBS)*11.589E-5*VIM**.667/(BRPM*BRPM*FNBPS) 

GC=ERR/GC 

JV=J V*1 

WRITE (KO, 204)  JV,  GC 

204  F0RMATC5X, 'PASS  NO. ',13,'  THROUGH  VELOCITY  LOOP,  CORRECTION  =', 

1 E12.5) 

WRITE (fc,175)Q6G,BRPM»FNBS 

175  FORMAT(1hO,3X,'Q6G=',E15.7,5x,'PRPM=',e15.7,5X,'FNBS=',E15.7) 

GC  =0 .0 

IF  (DAf-S  (GC  ) .LF.  VLIM)  C-0  TO  502 
IF  (JV  .GE.  NLIM)  GO  TO  503 
GO  TO  506 

IF  (GC  .LT.  0.  .AND.  G 6G / ( B R P M * F NB S ) .LT.  .46875)  GO  TO  504 
506  CONTINUE 
VIM=VIM*GC 
VI=VIM*3,281 
GO  TO  501 
504  WPITE(KO,205) 

205  FORMAT(50H****  FLOWERS  ARE  CHOKED  - PROCESS  NEXT  CASE  ******) 

GO  TO  2 

503  WRITE (KO, 207)  GC 

203  FORMAT ( 5 2 H * * * * FAILURF  TO  CONVERGE  IN  VELOCITY  LOOP  - ERROR  * , 

1 E12.5.6H******) 

502  CONTINUE 
T03G=T02G 
T 04  G = 7026 
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T05G=T02G 
T0£G=T076 
T09G=T07G 
T01 0G=TO7G 
TO  1 1 G = T 0 1 G 

BLPR  = WG*  .5*  (CP6G+CP7G)  * T 0 6 G * ( P p P * * C (1  . - 2 • / CGAM6G+GAM7G)  ) 

1 /ATAb>-1 ,)*778./550. 

BLPPM  = BLPR*  . 7457 

DPI=4.317E”4*GI*6I*FI*DL(1)*0NP(1)/(R0LI*DDH(1)) 

DPII=4.317E-4*GII*GII*FII*DL(2)*DNP(2)/(R0LII*DDHC2)) 

DP=4.317E-4*G*G*F*DL(3)*DNP(3)/(ROL*DDH(3)) 

R0A1=144.*PAT*Aml/(1545.32*TA1) 

ROA2=POA1 
LL  = 0 

37  Z 1 = R 0 A 2 

PA2=PAT-GA*GA*((1.+DSIG(3)**2)*<ROA1/ROA2-1.)+ 

1 AF*DALG(3)*DW(3)*ROAl/(.5*DSir-  (3)*(R0A1+RCA2>))/ 

2 (ROA  1 *2  . *32  . 1 74*1  44  .) 

ROA2=144.*PA2*AML/(1545.32*TA2) 

L L = L L * 1 

Z3=R0A2-Z1 

IF  (DAPS(ZT)  .IT.  ROHM)  GO  TO  38 
IF  (LL  . L E . M IM)  GO  TO  37 
wRI TE (KO,401  ) Z 3 

401  FORMAT (49H****  FAILURE  TO  CONVERGE  IN  AUX.  DENSITY  LOOP  - *, 

1 'ERROR  =',E12.5,6H******) 

GO  TO  39 

38  WRITE (kO,  202  ) LL,  Z 3 

202  FORMAT  MCX  , 'TOTAL  PASSES  THROUGH  AUX.  DFNSITY  LOOP  =',I3, 

1 ' - ERROR  =',F12.S) 

39  PF=WA*CPA2*TA2*((PAT/PA2)**m.-1./GAM2A)/ATAP)-1.)*776./550. 
PFM= .745 7*PF 

R06G=.C93178*P6G*GMU/TcG 

QCOPR=60.*WG*DSQRT(1.4*GMU*518.7/(GAM6G*28.(57*T06G))/R066 
DO  99  KS  = 1 , 16 
99  SVR ( K S ) =BMT  (KS  ) 

QDI=WI*1800.*(CP2+CP4)*(T4-T2) 

QDII=U 11*1800. *<CP2+CP6)*(T6-T2) 
QDA=W3*180Q.*(CP3*CP1)*(T3-T1) 

CALL  OUTPUT 
GO  TO  2 

95  DO  96  K$  = 1 ,16 

96  SVR(KS)=8MT(K$) 

GO  TO  10 

91  WRITE(K0,102) 

102  FORMAT('  SINGULAR  MATRIX') 

GO  TO  2 
END 
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) • ENPUT 

SUBROUTINE  FNPUT 
IMPLICIT  REAL*?  ( A -H , 0-7  ) 

DIMENSION  SVR ( 1^) 

DIMENSION  fcSVP (16) 

COMMON  /VAR/  TC1G,  T02G,  TOtG,  T07G,  T1,  T2,  T3,  T4,  T6,  T6, 

A TIG,  T 1 4 , TU,  T18,  T2G,  TA2, 

P T0  3G  , TO  4G , TO  5 G , TO  6G  , T0  9G,  TOlPG,  T011G,  T5,  T7,  T9,  Til, 

C T 1 2 , T13,  T 1 5 , T 1 7 , T 1 V , TIG,  T2G,  T3G,  T4G,  T 5 G , T 6 G , T 7G , 

D TfcG,  T9G , T 1 0 G , T11G,  T A 1 , EFFL,  P01G,  P02G,  P03G,  P 04G,  P05G, 

E PC6G.  P07G,  POfcC,  P09G,  P010G,  P011G,  PIG,  P 2 G , P 3G , P4G,  P 5 G , 

F peG,  P7G,  PfiG  , P 9 G , PlGC,  P11G,  P(2D),  AC(20),  VC(20),  F M 1 , 

G F M 2 , F M 3 , F*4,  F M 5 , F M 6 , FM7,  FM8,  F M 9 , FM10,  FM11,  A(11), 

H GAMIG,  G A M 2 G , G AM  3 G , GAM A G , G A M 5 G , G A M 6G  , GAM7G,  GAMfiG,  G A M 96  , 

I G AMI  LG , GAM11G,  T GM I , T 6 M 1 1 , TM  A , GIMU,  GI1MU,  AMU,  PRGI, 

J PRGI 1 , PRA,  CPGI,  CPuII,  CPA,  TM I , TMII,  TM,  7Imu,  ZIIMU,  ZMU, 

K PRI,  PRII,  PR,  CPI,  CPII,  CP,  GGI,  G G I I , 6 A,  Gl,  G I I , G,  PEGI, 

L REGII,  REA,  RFI,  R E I I , PE,  GF1,  G F I I , AF,  FI,  F 1 1 , F,  GJI, 

M GJII,  AJ,  ZJI,  ZJII,  ZJ,  HGI,  HGII,  HA,  HI,  H I I , H,  CGI,  C G I I , 

N CA,  Cl,  C I I , C,  UI,  UII,  U,  A G I , A G I 1 , AA,  TUNI,  TUMI,  TUN, 

0 TORI,  TOWII,  TOW,  EPSI,  EPSII,  EPSA,  ESSI,  ESSII,  ESS,  WG, 

P Q6G,  QCORR,  BPR,  ALT,  PAT,  PA?,  RH,  FN(4),  EFFLP,  VI,  VIM, 

0 PIGA,  HIM,  CL1M,  WT1M,  D C M , bLPR,  BLPRM,  CPP,  CDPM,  PF,  PFM, 

R R 0 1 , P L K W , WI,  WII,  WM,  W V 1 , W V ’ , WE,  WS,  WA,  GP,  QCM, 

S QM,  QV1,  Q V 2 , QE,  GS,  CF23,  CF34,  CF45,  CF78,  CF89,  CF910, 

COMMON  /VAR/  CF111,  ATAP,  Hi,  C L 1 , WT 1 , ACN(20),  DC,  DL(3), 

A D H ( 3 ) , D W ( 3 ) , DDH6(3>,  D A L G ( 3 ) , DAF0A(3)  , DSIG(3),  DDEL(3), 

P D A ( 3 ) , Dt((3),  D S L ( 3 ) , DSIC3),  D A L F ( 7 ) , DDH(3),  C0NVC16),  DELM, 

C ROLIM,  PLIM,  FG,  TITLF(2fc) , DP,  DPI,  DPII,  PPR,  PPRM, 

D T A 1 F , RHP,  TIGM,  DKHC3),  ETIG,  E T 0 2 G , ET06G,  ET07G,  E T 1 , 

E E T 2 , E T 3 , ETA,  E T 6 , E T 6 , ETIG,  ET14,  ET16,  ET16,  ET20,  E T A 2 , 

F E BPR , CWA,  D N P ( 3 ) , GDI,  QD1I,  Q D A , EVIM,  BRPM,  VLIM,  R0A1, 

G R0A2  , D AF AL ( 3 ) , DSLL (3 ) , AP ( 3 ) 

COMMON  /IVAR/  NTAM,  NCSGN,  NCS,  1C  I , K 0 , N L I w , N P < 7 ) , NFULL, 

1 NBS,  NbPS,  JCPRO 

EQUIVALENCE  (FSVR(1),  ETIG) 

EQUIVALENCE  (SVR  (1  ) ,T01G) 

NAMELI ST/FPT/  NTAM,  T A 1 F , PLKW,  ETIG,  ET02G,  ETOfcG,  E T 0 7 G , ET1, 

A E T 2 , E T 3 , E T 4 , E T 6 , ETA?,  LT14,  ETU,  E T 8 , ET10,  ET13,  ET2C, 

P EFFLP,  WI,  wll,  WB,  WM,  WVl,  WV2,  WE,  WS,  CWA, 


C 

QP, 

(j  C M , Q M , 0 V 1 

, Q V2 , QE  , QS  , 

CF23, 

CF34, 

D 

C FAS 

, C F 7 8 , CF89, 

C F91 0 , C F 1 1 1 

* 

E 

EVIM, 

ATAP,  PIGA, 

HIM,  CLl M , WT 1 M 

, F N , 

A CN  , 

F 

DL  , 

DH,  DW,  DDHG 

, dalg,  dafoa. 

DSIG  , 

DDEL  , 

G 

D A , 

DKH,  DSL,  DSI  , 

DALF , DDH  , 

NP  , 

H 

PAT 

, RHP , E b P K , 

A,  CONV,  D F L M , 

N D S G N 

, NCS 

I NLIM,  ROLIM,  PLIM,  FG,  BRPM,  VLIM,  NES,  NBPS,  JCPRO,  ROl, 

J OAFAL.DSLL,  AP 

KI  = 5 
K C = 6 

READ(K'I,2DG)  (T  TTLE  ( J ) , J = 1 , 14) 

READCVl  ,?0G)  (TITLE(j)  , J = 1 3 , 2 p ) 

200  F0RMAT(13A6,A?) 

RE  AD (K I , FPT ) 

T A 1 = T A 1 f+459,69 
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T I G = T A 1 + 20. 

T02G=TA1+210. 

T06G  = TA1 +20  . 

T07G=T Al + 70. 

T 1 = T A 1 +2C  . 

T2  = T Al +20  . 

T3  = TA1 ♦AG  . 

T A = T A 1 +55  . 

T6=TAl+30. 

T P = T A 1 +65  . 

T 1 C = T A 1 +65  . 

T 1 A = T A 1 +30  . 

T1 6 = TA1 +75  . 

T18=TA1+65. 

T2G=T A1 +65  . 

TA2  = TA 1 +30  . 
b PR  = 1 .17 
VIK=1 2C . 

DO  30  <=1,16 

30  IF  (ESVR(K)  .GT.  0.)  SVR  (K)=ESVR(K) 


I F 

( EPPR 

• G T * C . ) 

SPR=EBPR 

I F 

( LV  IM 

.GT.  0.) 

V I M = E V I M 

I F 

UTAH 

.LT.  C) 

call  EXIT 

I F 

( NT  AM 

.LE . 0) 

GO  TO  10 

DO 

20  J«1 

,19 

20  SVR (J )=TA1 
10  CONTINUE 

PIG=PIGA*1 A. 7 
RH=RHP* .01 
E F FL  = F F FLP* .01 
VI=VIM*3.2£1 
DC=DCM+.03281 
Hl=HlM*. 03281 
wT1=wTlM+.032fi1 
CL  1 = C L 1M* .032R1 

CALL  AIRPCTA1,  FAT,  T A 1 , PH,  X,  X,  X,  X,  A V L ) 

R0Al=PAT*AML*.C9317/TAl 

WA=ROAl*CWA*DSOPT(TA1/51fc.7)/6U. 

DK  ( 1 ) = DKH  ( 1 > /'1'600. 

DK(Z)  = DKH (2 ) /7600. 

D k ( 3 ) = dkh(3)/3600. 

DNP ( 1 ) =NP ( 1 > 

dnp<2)=np(2> 

DNP(3)=NP(3) 

A ( 1 ) =H 1 +CL  1 
G6G  = V I * A < 1 ) *6P  . 

WRITE  ( K 0 , 102) 

102  FCRMAT('I') 

WRITE (KO , FPT ) 

101  FORMAT ( 21 X , 1 AA6  } 

WR I TE (KO , 1 01  ) TITLE 

R F T U R N 

END 
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• GASP 

SUBROUTINE  GASP(T,  FN,  GAM,  CP,  PR,  UM,  GMU  ) 

IMPLICIT  R E A L * 6 (A-H.O-Z) 

DIMENSION  F N ( 4 ) , W M ( 4 ) , T B l 4 ) , CCA),  U(4),  F K ( 4 ) , 

1  GCA),  X ( A ) , CN(A),  S C 4 ) 

1 = HE 

2 = N2 

3 = C 0 2 

4 = CO 

DATA  Tb  /4 . 21  6D0, 77 . 4DU,  1 94. 66D0, 81  . 5DO/ 

DATA  WM.  /4.P0D0,28.C16D0,  44.O1CD0,  28.010DC/  , R /1.OA646D0/ 
C ( 1 ) = 1 .2406 
IF  (T  . L E . 650.) 

1 FK(1)=(.8+(T-475.)*(.97-.&)/175.)*.1 

IF  (T  .GT.  650.) 

1 FK(1)=(.97+(T-650.)*(1.13-.97)/200.)*.1 

FK ( 1 ) = FK  ( 1 ) / 3 600 . 

G ( 1 ) = 1.6667 

U(1)  = 1 .12+(T-400.)*(1  . 9 c 5 — 1 .12)/5C0. 

U C 1 ) =U ( 1 ) *1 . E-5 


IF  (T  .LE.  600) 

1 C(2>=3.5G2+(T-400.)*(3. 5065-3. 5C2)/200. 

IF  (T  . G T . 600.  .AND.  T .LE.  800.) 

1 C(2)=3.5G65+(T-600.)*(3. 5315-3. 5065)/200. 

IF  (T  .GT.  POO.) 

1 C(2)=3.5315*(T-800.)*(3. 559-3. 531 5 ) / 1 00. 

C C 2 ) = C (2) *R/WM  (2) 

IF  (T  .LE.  650.) 

1 FK(2)=.925+(T-45G.)*(1.265-.9?5)/200. 

IF  (T  .GT.  650.) 

1 FK(2)=1.265+(T-650.)*(1. 573-1. 265)/?00. 

FK(2)=FK(2)*1 .4F-2/36CG. 

G ( 2) =1  .4 

IF  (T  .CT.  600.) 

1 G(2)=1.4*(T-6C0.)*(1. 391-1. A)/300. 

IF  (T  .LE.  550.) 

1 U(2)  = .933XT-450.)*(1.087-.933)/1  0C. 

IF  (T  .GT.  5 5 n . .AND.  T .LE.  725.) 

1 U(2)=1.0o7+(T-550.)*(1.33-1.U87)/175. 

IF  (T  .GT.  725) 

1 U(2)=1.33*(T-725.)*(1.545-1.33)/175. 

U ( 2 ) =U ( 2 ) * 1 . 1 1 72E-5 


IF  ( T .LE . 600. ) 

1 C ( 3)=4.19S+(T-4?0.)* (4.653-4.195) /I  50. 

IF  (T  .GT.  600.  .AND.  T .LE.  750.) 

1 C (3 )*4.653* (T-600. )* (5.045-4.653) /I 5C. 


102 


IF  (T  . 6 T . 750.) 

1 C(3)  = 5.045«’(T-750.)*(5.?65-5.045)/15G 

C(3)=C(3)*R/W*(3) 


IF  ( T . LE  . 600. ) 

1 FK(3)  = .82-KT-425.)*(1.32-.£?)/175. 

IF(T.LE.140.)Fk(3)=5.714E-3 
IF  (T  , G T . 600.) 

1 FK(3)=1.3?*(T-600.)*C2. 12-1.32)7250. 

Fk(3)=Fk(3)*8.407e-,/360G. 

IF  (T  .LE . 450.) 

1 G < 3 ) = 1 .334+(T-4C0.)*(1  .315-1 .334) /5C. 

IF  (T  .GT.  450.  .AND.  T .LE.  550.) 

1 G(3)=1.315+(T-450.)*C1. 266-1. 315)/1C0. 

IF  (T  .GT.  550.  .AND.  T .LE.  650.) 

1 G(3)  = 1.286*(T-550.)*(1.  2645-1.  2 8 6 ) / 1 00. 

IF  (T  .GT.  650.  .AND.  T .LE.  750.) 

1 G(3)  = 1 .2645+(T-650.)*(1  .248-1  .2645) / 1 0 0 . 

IF  (T  .GT.  750.  .AND.  T .LE.  850.) 

1 G(3)  = 1.248+(T-750.)*C1.  235-1.  248) /I  00. 

IF  (T  .GT.  0 5 0 . ) 

1 G(3)  = 1.235+(T-850.)*(1.  2295-1.  2 3 5 ) / 5 0 . 

IF  (T  .LE.  575.) 

1 U(3)=.9175+(T-450.)*(1.152-. 9175)7125. 

IF  (T  .GT.  575.  .AND.  T .LE.  700.) 

1 u(3)=1.152+(T-575.)*(1. 375-1. 152)/1?5. 

IF  (T  .GT.  700.) 

1 U ( 3 ) = 1 .375+CT-700.)* (1 .541-1 .375) /100. 
U(3)=U(3)*9.2067E-6 

FT=FN(1)+FN(2)+FN(3) 

N G A S = 3 

IF  ( F N ( 4 ) .LE.  0.)  GO  TO  10 
FT=FT+FN(4) 

NG  A S = 4 

IF  (T  .LE.  540.) 

1 C(4)=3.50t+(T-400.)*(!. 506-3. 503) /14C. 

IF  (T  .GT.  54".  .AND.  T .LE.  630.) 

1 C(4)=3.506+(T-540.)*(3. 513-3. 506) / 90. 

IF  (T  .GT.  630.  .AND.  T .LE.  720.) 

1 C(4)=3.513+(T-630.)*(3. 529-3. 513)790. 

IF  (T  .GT.  720.  .AND.  T .LE.  810.) 

1 C(4)=3.529+(T-720.)*(3. 552-3. 529)/90. 

IF  ( T . G T . 810.) 

1 C(4)=3.552+(T-810.)*(3. 583-7. 552)790. 

C ( 4 ) = C C4) *.0708989 

IF  (T  .LE.  540.) 

1 U(4)=.84*(T-400.)*<1.075-.64)/140. 

IF  (T  .GT.  540.  .AND.  T .LE.  720.) 

1 U(4)=1 ,075+(T-540. )*(1 .34-1 .075) /180. 

IF  (T  .GT.  720.) 

1 U(4)  = 1.34+(T-720.)*(1.  575-1.  34  ) /I »0  . 

U (4) =U(4) *1 .1 1 3 2 E - 5 


IF  (T  .Lfc . 530  . ) 

1 6 (4)  =1  .4 

IF  (T  .GT.  530.) 

1 G(4)=1.4*(T-530.)*(1.7&P5-1. 0/370. 

IF  (T  .LE.  540.) 

1 FK(4)  = .?25t-(T-400.)*(1.0e7-.625)/1  40. 

IF  (T  , G T . 540.  .AND.  T .LE.  720.) 

1 FK(4)  = 1.0R7-*-(T-540.)*(1.  365-1. 087  ) /I  80. 

IF  (T  .GT.  7?0.) 

1 FK(4)  = 1.3Q5+(T-720.)*(1.  665-1.  385  ) /I  80. 

FK(4)=FK<4)*. 01342/3600. 


0 CP=0. 

D A N = 0 • 

GAM  = CJ. 

UN  = D . 

G M U = 0 • 

Z3  = C. 

DO  11  1=1  » N G A S 
su)  = i.5*re(i)*i,8 
X ( I ) = FN  ( I ) / FT 

1 GNU=6MU*X < I ) *WW U ) 

DO  50  1=1 fNGAS 

CN ( I ) =X ( I ) *WM ( I ) /GMU 
CP=CP*CN(I)*C(I) 

GA  N = C N ( I ) *C ( I ) +GAM 

DAM=CN(I)*C(I)/G(I)+DAN 

Z1=C. 

Z 2 = 0 . 

DO  30  J=1 .NGAS 

IF  ( J .NE . I ) GO  TO  12 

P = 1 . 

A M = 1 . 

GO  TO  29 

2 D=X ( J ) / X ( I ) 

E = DSGRT (S  C I ) *S ( J ) ) 

UR=U(I)/U(J) 

WRl=DSGRT(6.*(1.-t-WM(I)/WM(j))) 

wR  = (w* ( J ) /w* ( I ) ) **  . 25 

P=D* (1 .♦DSQRT (UR)*WR)**?/WRl 

AR'«.25*D*(1.+DS0RT<UR*WR**3*(T+S(I))/(T+S(J))))**2*(T+E)/(T*S 
9 Z 1 = Z 1 ♦ P 
Z2=Z2+A* 

0 CONTINUE 

Z3  = Z 3+FK  ( I ) /Z  ? 

UM  = UM*U  (I ) n 1 
0 CONTINUE 

GAM=GAM/DAM 

PR=CP*UM/Z3 

RETURN 

END 
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MRP 

SUBROUTINE  AIRP(TAMB,  PAMB,  T,  RH,  P,  AMU,  CPA,  GAM, 
IMPLICIT  REAL*8  (A-H,0-Z> 

T2=DEXP( -3  o.  098  73  ♦.  1C57402MAMB-.  6 74  864  E -4  *TA  *6  *TAMP) 
WF=RH*18.016*T2 / (28 .966* (PAMq-T?) ) 

T?  = T*T 

P=.5552+.1346?3E-2*T-.03lE-4*T2*.20668E-8*T*T2 

AMU=1.E-5*.175E-7*(T-400.) 

CPA=.2421-.0015E-2*T*2.E-8*T2 

GAMA=CPA/(CPA-.C685545> 

CP«= .4509-.0054F-2*T+. G0b8E-4*T2 
GAMW=CPW/(CPW-. 11022) 

A=WF*CPW 
B= C 1 ,-WF) *CPA 
CPA=WF*CPW+(1.-WF)*CPA 
GAM=(A>B)/(A/GAMw  + 5/  GA.MA) 

AML=1P.016*WF+28.966*(1.-WF) 

RETURN 

END 


1 

AML  ) 
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MAC  H 

SUBROUTINE  MACHCDELM,  NLIM,  POG,  GA*,  EM,  TOG,  W6,  GMU,  A,  P6, 
1 TG,EN,CP) 

IMPLICIT  REAL*P  (A-H,0-2) 

DIMENSION  FNU) 

DATA  KC  /fc/ 

1 = 0 

1C  F MG  = F M 

TG=1.+FMG*FMG*(GAM-1.)/2. 

WRITE  (6, 1 1 )P0G,TG,GAM 

11  F0PMAT(1h0,3x,'P0G  = ',E15.7,5x,'TG=',E15.?,5X,'GAM,=  ',E15.7) 

PG=P0G/TG**(GAM/(GAM-1.)) 

TG=T0G/TG 

FM=WG*DSQRT(TG/(GAM*GMU))/(20.7774*PG*A) 

WPITE(6,12)PG,TG,FM 

12  F0RMAT(1HC,3X,'PG=',E15.7,5X,'TG=',E15.7,5X,'FM=',E15.7) 

CALL  GASPCTG,FN,GAM,CP,DUM,DUW,GMU) 

1 = 1 + 1 

FMG1 =FMG-FM 

IF  CDAbS(FMGl)  .LT.  DELM)  RETURN 
IF  (I  .LT.  NLIM)  GO  TO  ID 
WRITE(KO.IOO)  FMG1 

00  FORMAT ('  FAILURE  TO  CONVERGE  IN  MACH  NO.  SU°ROUTINE', 

1 ' - ERROR  =',  F 1 2 . f ) 

RETURN 

ENO 


QP 


r ^ 

LIQP 

SUBROUTINE  LIQP(T,  FG,  UM,  PR,  CP) 

IMPLICIT  REAL*8  (A-H.O-Z) 

C = T / 1 .6-273. 

F = F G * F G 

CP=1 .-.465*FG+C*( .0011 5-.001G5* (1 ,-FG)**3.Qc96) 

S F = 1 301.  / (998. 33+8. 1£55*(C-2u.)  + . 00585*  (C-20. >**2) 

Ul»  = DEXP(-3.01174t-.49287*FG+(2.30854  + 1.1  8556*FG*.8  06*F)*SF) 
FK=137.-82.91667*FG+17.91667*F 

FK=1.*(FK+C*(166.-16u.8833*FG+49.0833*F-FK)/100.) 

UF*  = 0. 01  488  1 

PP  = 1 .0E3*UV*CP/ FK 

U*=UM* .000672 

RETURN 

END 


S 1 M Q 

SUBROUTINE  SI"G(A,  B,  N,  KS) 

IMPLICIT  P E A L * P (A-H,0-Z) 

SIMC.  SOLVEb  A SYSTe"  CE  L I N E A R EQUATIONS  AX  = P 
INPUTS 

A - MATRIX  OF  COEFFICIENTS  STORED  COLUMNWISE.  DESTROYED  IN 
COMPUTATION 

B - VECTOR  OF  ORIGINAL  CONSTANTS.  DESTROYED  IN  COMPUTATION 
SOLUTION  IS  RETURNED  IN  o 
N - NUMBER  OF  EQUATIONS  AND  VARIABLES 

KS  - OUTPUT  FLAG  0=N0R*AL  1=SINGULAR  SET  OF  EQUATIONS 

DIMENSION  A ( 1 ) , B (1  ) 

FORWARD  SOLUTION 

T 0 L = Cl  . 0 
KS  = 0 
J J = -N 

DC  6 i J = 1 , N 
J Y = J + 1 
J J=J J+N+1 
BIGA=0.0 
I T = J J - J 
DO  30  I = J ,N 

SEARCH  FOR  MAXIMUM  COEFFICIENT  IN  COLUMN 
I J = I T ♦ I 

IF  (DABS (B IGA ) -DABS ( A ( 1 J )) ) 20,30,30 
20  BIGA=A(IJ) 

I M a X = I 
30  CONTINUE 

TEST  FOR  PIVOT  LESS  THAN  TOLERANCE  (SINGULAR  MATRIX) 

IF  (DABS(BIGA)-TPL)  35,  35,  AO 
35  KS=1 
RETURN 

INTERCHANGE  ROWS  IF  NECESSARY 

AO  11=J+N*(J-2) 

IT=IMAX-J 
DO  50  K=J,N 
I 1 - I 1 ♦ N 
I 2 = 1 1 ♦ ! T 
SA VE  = A ( I 1 ) 

A(I1)*A(I2) 

A ( I 2) -SAVE 

DIVIDE  EQUATION  BY  LEADING  COEFFICIENT 

50  A(I1)=A(I1)/BTGA 
S A V E = P ( IMAX) 
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B C IMAX  ) =B  ( J ) 
B(J)=SAVE/BIGA 


ELIMINATE  NEXT  VARIABLE 

IF  (J-N)  55,  7C,  55 
55  I0S=N*(J-1) 

DO  65  IX  = J Y , N 
IXJ  = IfiiS+IX 
I T = J -I  X 
DO  60  J X = J Y , N 
IXJX=N*(JX-1)+IX 
J J X = I X J X*I  T 

60  A(IXJX)=A(IXJX)-(A(IXJ)*A(JJX)) 
65  B (I  X ) =e ( IX ) -CB  ( J )*A  ( IX J ) ) 

BACK  SOLUTION 

70  N Y = N - 1 
I T = N * N 

DO  80  J = 1 ,NY 
I A=I T-J 
I 6 = N - J 
I C = N 

DO  80  K = 1 , J 

B(IE)=B(IB)-A(IA)*E)(IC) 

I A = I A-N 
»C  I C = I C - 1 
RETURN 
END 


TPUT 
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.OUTPUT 

SU GROUT  I N E OUTPUT 
IMPLICIT  R £ AL  * 6 (A-H.O-Z) 

DIMENSION  VG(11),  Q G ( 1 1 ) , OC(?U),  GCG(20),  DPC(2U>,  DPCR(2Q) 
DIMENSION  P D G ( 1 1 ) , PP(11),  P C ( 2 P ) , P0(11),  p 0 G ( 1 1 ) 

C OM  M ON  /VAR/  TOIG,  T026  , T06G,  T07G,  T1,  T2,  T3,  T4,  T6,  T8, 

A TIO,  T14,  T 1 6 , T18,  T20,  T A 2 , 

P TO  3G  * T04G  , T0  5G,  T0  8G  , T0  9G  , TOIOG,  TOl 1 G , T5,  T7  , T9,  Til, 

C T 1 2 , T13,  T 1 5 , T 1 7 , T 1 9 , TIG,  T2G,  T 3G , T4G  , T 5 G , T6b,  T7G, 

D TPG,  T9G,  TIOG,  T11G,  T A 1 , F f F L , PQ1G,  P02G,  P07G,  P04G,  P05b, 

E POfcG,  P07G,  P08G,  P 0 9 G , P0106,  P011G,  PIG,  P2G,  P 3G , P 4 G , P 5 G , 

F P6G,  P 7 G , PPG,  P 9 G , PIOG,  P11G,  P(20),  A C ( 2 C ) , VC(20),  FMl, 

G FM2,  FM3,  FM4,  F M 5 , FM6,  FM7,  F M 8 , F M9  , F ^ 1 0 , F * 1 1 , A(11), 

H G A M I C , G AM  2 G , G A M 3 6 , G A M 4 G , 6 A M 5 G , G A m 6 G , G A M 7G  , G A M o C , GAM,  96, 

I GAM1CG,  GAM11G,  TGMI,  TGMII,  TM A , GIMU,  GIIMU,  AMU,  PRGI, 

J PRGI  1 , PRA,  CPGI,  CPGII,  CPA,  T M I , TMI  I , TM  , Z 1WU  , ZIIMU,  ZMu, 

K PRI,  PR  I I , PR,  CPI,  CPII,  CP,  GGI,  GGI  I , GA,  GI,  GII,  G,  &EGI  , 

L PFGII,  R £ A , REI,  REII,  RE,  GFI,  G F I I , AF,  FI,  FII,  F,  6JI, 

M GJII,  AJ,  ZJI,  ZJII,  ZJ,  HGI,  HGII,  HA,  HI,  HII,  H,  CGI,  C6II, 

N CA,  Cl,  CII,  C,  UI,  UII,  U,  A G I , AGII,  AA,  TUMI,  TUNII,  TUN, 

0 TOW  I , TOWII,  TOW,  EPSI,  EPSII,  FPSA,  ESSI,  ESSII,  ESS,  WG, 

P Q6G,  OCORR,  PPR,  AlT,  PAT,  P A 2 , RH,  FN(4),  EFFLP,  VI,  VIM, 

Q PICA,  HIM,  CL1M,  WT1M,  DCM,  BLPR,  BLPRM,  CPP,  CPP*,  PF,  PFM, 

P R 0 1 , PLKW,  WI,  WII,  *0,  WM,  W V 1 , WV?,  WE,  WS,  WA,  QP,  UCM, 

5 QM,  0 V 1 , Q V 2 , GE,  QS,  C F 2 3 , CF34,  CF45,  C F 7 8 , CF89,  CF910, 

COMMON  /VAR/  CF111,  ATAP,  Hi,  C L 1 , WT1,  ACN(20),  DC,  DL(3), 

A D H C 3 ) , 0 W ( 3 ) , D D H G ( 3 ) , DALG(Z),  DAF0A(3),  DSIG(T)  , 0 D E L ( 3 ) , 

B D A ( 3 ) , 0 K ( 3 ) , 0 S L C 3 ) , DPI  (39  , D A L F ( 3 ) , DDH (3  > , CONV ( 1 A) , PELM  , 

C ROLIM,  PLIM,  FG,  TITLE  (28),  OP,  DPI,  0 P I I , PPR,  P P R N , 

0 T A 1 F , RHP,  T I G M , D K H ( 3 ) , FTIG,  F T 02  G , ETCAG,  ET07C-,  E T 1 , 

E E T 2 , E T 3 , ET4,  E T 6 , E T 8 , FTIG,  ET14,  FT16,  ET18,  FT20,  E T A 2 , 

F EBPP,  CWA,  0 N P ( 3 ) , 001,  GCII,  Q 0 A , EVIM,  PRPM,  VLIM,  R0A1, 

6 ROA2,OAFAL(3),DSLL(3),AP(3) 

COMMON  / 1 V A R / N T AM  , NDSGN,  NCS,  KI,  KO,  NLIM,  NPd),  NFULL, 

1 NBS,  NBPS,  JCPRO 
WRITEOcO.IOO) 

WRITE (KO, 101 ) TITLE 
WRITE (K0,400> 

400  FORMATC'O  INPUT  LIST'/'  CAVITY  * LASER  SPECIFICATIONS'/ 

1 ' OF  SCR  I PT ION ' , 24X  , ' NAME  VALUE  UNITS') 

WPITE(KO,4(j3)  FIGA 

403  FORMAT ('  INLET  PRFSSURE  ',12('.  ' ) , ' P I G A ',612.6,'  ATM') 

WRITE(KO,404)  HIM 

4G4  F 0 R M A T ( ' INLET  HEIGHT  ' , 1 3 C ' . '),'HlM  ',612.6,'  CM') 

WRI TE (K0,405)  WT1M 

405  FORMAT ('  WIDTH  IN  FLOW  DIRECTION  ',«('.  '),'wT1M  ' , G 1 2 . b , 

1 ' CM') 

WRITE (KO, 4 06)  CllM 

406  FORMAT  ('  LENGTH  ',1fc('.  '),'CL1M  ',612.6,'  CM') 

WRITE ( K 0 , 4 0 8 ) DCM 

408  FORMAT ('  INLET  HYDRAULIC  DIAMETER  ',7('.  '),'DCM  ',G12.6, 

1 ' CM') 

WPITE(k0,4U7)  PLKW 

407  FORMAK'  LATER  OUTPUT  PC«FR  . . ',«('.  '),'PLKW  ',G12.C, 

1 ' K i»  ' ) 

hRITE  ( K 0 , 4 0 1 ) EFFLP 

401  FORMAT!'  LASER  EFFICIENCY  ' . 1 1 < ' . '), 'EFFLP  ',G12.6, 


no 


1 ' PERCENT') 

WRITE(K0,4G9)  NDSGN 

409  FORMAT('  DESIGN  OPTION  ',13('.  '), 'NDSGN  ',I7,5X,'  -') 

IF  ( N C S 6 N . N E . ’)  GO  TO  5 r 
WRITE  (K0.410)  NCS 

4 1 C FORMAT  ( ' NUMBER  OF  AREAS  (NDSGN=3)  ' ,  1 * * * * *  7 ( ' . '),'NCS  ', 

1 17, 5X,' 

WRITE  (KO, 41 1 ) (ACN(K)  ,K  = 1 ,NCS) 

411  FORM  AT  ( ' ARFAS  ',17('.  '),'ACN  ' , 6 ( G1  2 . fe  , 1 X ) , 

1 3 (51 X ,6(G12.6, 1 X) ) ) 

wRI TE (KO,  41 2) 

4 1?  FORMAT ( 6 5 X , 'FT  * * 2 ' ) 

50  WRITE(KO,413)  FN 

413  FORMAT!'  GAS  MIXTURE  ',14('.  '),'FN  '.4F6.2, 

1 ' HE  TO  N2  TO  C 0 2 TO  CO') 

WRITE(KO,414) 

414  FORMAT ( '01 M T I AL  ESTIMATE  OF  STEADY  STATE  QUANTITIES', 

1 3 6 X , ' I F = 0. .USES') 

WRITE  (KO, 41 5)  ETIG 


415 

FORMAT  ( ' 

GAS 

STATIC 

TEMPERATURE 

- STATION 

1',3 

('. 

') 

, 'ETIG 

9 

1 

012.6,' 

DEG  P 

',  10X,  ' 

T A 1 ♦ 20.') 

WRITE (KO, 

416)  E T 0 2 G 

4 1 6 

FORMAT  ( ' 

GAS 

TOTAL 

TEMPERATURE  - 

STATION 

2 ',3 

('. 

') 

, 'ET026 

» 

1 

C-  1 2 . 6 , ' 

DEG  R 

',10X,  ' 

T A 1 ♦ 2 1 C . ' ) 

WRITE (KO, 

417)  ET06G 

417 

FORMAK  ' 

GAS 

TOTAL 

temperature  - 

STATION 

6 

('  . 

') 

, 'ET06G 

0 

9 

1 

G 1 2 . 6 , ' 

DEG  P 

' , 10X  , ' 

T A 1 + 30.') 

WRI TE (KO, 

41*)  ET07G 

4 18 

FORMAT ( ' 

GAS 

TOTAL 

TEMPERATURE  - 

STATION 

7 ',3 

('. 

') 

, 'ET07G 

0 

9 

1 G 1 2 » 6 , ' DEG  P',10X,'TAl  + 70.') 

WRI TE (K0.419)  ET1 

WRITE (KO, 420)  E T 2 

WRITE (KO, 421  ) E T 3 

WRITE  (KO,  422  ) E T 4 

WR I TE (KO , 42  3 ) E T 6 

WRI TE  (K0,4?4)  E T8 
WRITE (KO,  425  ) E T 1 0 
WRITE (KO,  426  ) E T 1 4 
WRITE (K0,42?)  ETIfc 
WRI TE (kO,42R)  E T 1 8 
WRITE ( K 0 , 4 2 9 ) E T 2 0 


<.19 

F ORMA  T ( 

0 

COOLANT 

Tt*Pf kATURF 

- 

STATION 

1 

0 

, 7 ( 

0 

') 

, ' E T 1 

0 

1 

G12.6 

9 

0 

DEG  P ' , 1 OX 

,'TA1  ♦ 20. 

') 

420 

FORMAT  ( 

0 

coolant 

TEMPER ATURE 

- 

STATION 

2 

0 

. 3 ( 

0 

' ) 

, ' E T 2 

* 

1 

G12.6 

9 

0 

DEG  R ' , 1 OX 

, ' T A 1 ♦ 20. 

') 

421 

FORMAT  ( 

0 

COOLANT 

TEMPERATURE 

- 

STATION 

3 

0 

, 3 ( 

0 

') 

,'ET3 

* 

1 

G1  2 .6 

9 

0 

DEG  P ' , 1 0 X 

, ' T A 1 4 4Q. 

') 

422 

FORMAT  ( 

0 

COOLANT 

TEMPERATURE 

- 

STATION 

4 

0 

,7( 

0 

') 

, ' E T 4 

0 

1 

612.6 

9 

0 

DEG  P ' , 1 OX 

,'TA1  ♦ 55. 

') 

423 

FORMAT  ( 

0 

COOLANT 

TEMPERATURE 

- 

STATION 

6 

* 

, 3 ( 

0 

') 

, 'ETfc 

0 

1 

612.6 

9 

0 

DEG  P ' , 1 OX 

, ' T A 1 ♦ 3U. 

') 

4 24 

FORMAT  ( 

0 

COOLANT 

TEMPERATURE 

- 

STATION 

8 

0 

, 3( 

* 

') 

, ' E T 8 

0 

1 

G 1 2 . 6 

9 

0 

DEG  * ' , 1 0 X 

,'TA1  ♦ 65. 

') 

425 

F C R N A T ( 

0 

coolant 

TEMPERATURE 

- 

STATION 

10 

0 

, 3 ( 

* 

') 

, 'ET10 

0 

1 

G1  2 .6 

9 

0 

DEG  R'.IQX 

, ' T A 1 ♦ 65. 

') 

426 

FORMAT ( 

0 

COOLANT 

TEMPFRATURE 

- 

STATION 

1 4 

0 

, 3 ( 

0 

') 

, ' E T 1 4 

0 

1 

612.6 

9 

0 

DEG  P ' , 1 0 X 

, ' T A 1 ♦ 30. 

') 
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427 

F ORM AT ( ' 

COOLANT  TEMPERATURE 

- 

STATION 

16 

' , 3 ( ' . ') , ' E T 1 fc 

1 

G 1 2 . 6 , ' 

DEG  R',10X,'TA1  ♦ 

75. 

') 

428 

FORMAT  ( ' 

coolant  temperature 

- 

STATION 

18 

' , 3 ( ' . '),'ET18 

1 

G 1 2 . 6 , ' 

DEG  P ' , 1 0 X , ' T A 1 ♦ 

65  . 

') 

429 

FORMAT  ( ' 

COOLANT  TEMPERATURE 

- 

STATION 

2 C 

' , 3 ( ' . '),'ET2C 

1 

C 1 2 . 6 , ' 

DEG  R ' , 1 0 x , ' T A 1 ♦ 

65. 

') 

WRITE (KO, 

430)  E T A 2 

430 

F ORMAT ( ' 

AIR  TEMPERATURE 

AFTER 

RADIATOR 

* 

,4('.  ' ) , ' E T A 2 

1 

612.6,' 

DEG  P',10X,'TA1  ♦ 

30. 

') 

W R I T E ( K 0 , 

431 ) E F P R 

431 

F ORMAT  ( ' 

FLOWER  PRESSURE 

RATIO 

' , 9 ( ' . ' 

),' 

ERPR  ',612.6, 

1 

* 

-' , 14X , '1 .1 7') 

WR I TE (KO , 402  ) EVIM 

402  FORMAT!'  INLET  VELOCITY  ',12(  . '),'FVIM  ',612.6,'  M/SEC 

1 1 OX , ' 1 20  . ' ) 

WRITE  ( K 0 , 4 3 2 ) 

432  FORMAT ( 'OCOOL ING  SYSTEM  FLOWS') 


WR  I TE  (KO  , 4.J3  ) WI 

433  FORMAT ( ' HEAT  EXCHANGER  1 (LIQUID) 

1 G 1 2 . 6 , ' LB/SFC') 

WRITE (KO, 434)  WII 

434  FORMAT ( ' HEAT  EXCHANGER  2 (LIQUID) 

1 G 1 2 • 6 , ' LB/SFC') 

WRITE (KO, 435)  WP 

435  FORMAT ( ' MIRROR  BYPASS  ',13('.  '), 

1 G 1 2 . t , ' LB/SFC') 

WRITE (KO, 436)  WM 

436  FORMAT ('  mirror 

1 612.6,'  LB/SFC') 

WRITE  (KO,  4 37)  WV1 

437  FORMAT ('  VACUUM  PUMP  1 

1 G 1 2 . 6 , ' LB/SEC') 

WRITE (K0,43P ) WV2 

438  FORMAT!'  VACUUM  PUMP  2 

1 G12.6  ,'  LB/SEC') 

WRITE (KO, 430)  WF 

.39  FORMAT ('  F-PEAM 

1 G 1 2 . 6 , ' LB/SEC') 

WRITE(KO,440)  WS 


' , 16  ( ' . ')  , 'WM 

. 1 3 ( ' . '), 
* 1 3 ( ' . ')  , 
, 1 6 ( ' . ' ) , ' W E 


40  FORMAT!'  SUSTAInER  ' , 1 5 ( ' . '),'wS 

1 G 1 2 . 6 , ' LB/SEC') 

WRITfc(K0,441)  CWA 

.41  FORMAT('  CORRFCTED  FAN  FLOW  (AIR) 

1 612.6,'  C F M ' ) 

WR1 TE (KO  , 445  ) 

A 4 5 FORMAT ( 'GSPEC I F ICATIONS  OF  BLOWERS') 
WRITE (KO , 442)  ATAP 


' , 7 ( ' . ') 


wl 


',?('.  '),'WII 

'we  ' , 

% 

'w v 1 ', 

'W  V 2 
* 

» 

* 

t 

',7 ('.  ' ) , ' C W A 


t 


f 


t 


t 
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442  FORMAT ( ' BLOWER  POLYTROPIC  EFFICIENCY  *,M'.  ' ) , ' A T AB 

1 612.6,'  -') 

WRITE(K0,446)  BPPM 

446  FORMAT  ('  BLOWER  SPEED  '),'BRPM  ',612.6,'  RPM') 

WRITE (KO, 447)  N*S 

447  FORMAT  ('  NUMBER  CF  eLOwER  SETS  '),'NBS  ',I6,6X, 

1 ' -') 

WRITE (KO, 448)  NPPS 

448  FORMAT ('  NUMBER  OF  BLOWERS  PER  SET  ',7('.  '),'NEPS  ', 


WRITE(KO,  543  ) VLIM 

543  FORMAT('  INLET  VELOCITY  TOLERANCE  ',7('.  '),'VLIM 

1 612.6,'  1 / S E C ' ) 

WRITE (K0.444) 

444  F OPM AT ( 'OTUNNEL  AREAS') 

WPI TE ( K 0 , 4 4 3 ) A 

443  F ORMAT ( ' STATION  1 TO  STATION  11  ',*('.  '),'A  ', 

1 6(1X,G12.6  )/50X,5(lX,Gl2.6  ),'FT**2') 

WRITE  (KO.IOO) 

WRITE  ( K 0 , 1 G 1 ) TITLE 
WRITE ( K 0 , 5 00 ) 

5 OC  FORMAT ('0  INPUT  LIST  (CONT.)') 

WRI TE (KO, 501 ) 

501  FORMAT ( 'OHEAT  EXCHANGER  S P E C I F I C A T I ON S ' , 2 3 X , ' E X C H AN G E R 1',5X, 

1 'EXCHANGER  2 ' , 7 X , 'R A 0 I A TO R ' , 3 X , ' UN I T S ' ) 

WRITE(K0,5G?) 

502  FORMAT (OX , 'DESCRIPTION', 25X,'NAME',8X,'(1)',13X, '(2)', 14X , '(3) ') 
WRITE (NO, 507)  DL 

503  FORMAT ('  CORE  LENGTH  (GAS)  ',11('.  '),'DL  ', 

1 3(2X,G12.6,2X),'FT') 

WRITE (KO,  504  ) DH 

504  FORMAT ( ' COPE  HEIGHT  ',14('.  '),'DH  ', 

1 3(2X,G12.6,?X),'FT') 

WRITE (KO, 5C5)  OW 

505  FORMAT ('  CORE  WIDTH  (LIQUID)  ',10('.  '),'DW  ', 

1 3(2X,G12.6,2X) , ' F T ' ) 

WRITE(K0,506)  DDHG 

506  FORMAT ('  HYDRAULIC  DIAMETER  (GAS)  ',7('.  ') , 'DDHG  ', 

1 3(2X,G12.6,2X) ,'FT') 

WRITE  (KO, 50?)  DALG 

507  FORMAT ('  SURFACE  AREA  / VOLUME  (GAS)  ',*('.  ' ) , ' D A LG  ', 

1 3 ( 2X ,G1 2 .6, 2X  ) , '1  / FT  ') 

WRITE (K  0 , 5 OR  ) DAFOA 

508  FORMAT ( ' FIN  AREA  / SURFACE  AREA  (GAS)  ',5('.  '), 'DAFOA', 

1 3(2X,G12.6,?X) ,'-') 

WRITE(kO,509)  DSIG 

509  F 0 R M A T ( ' FREE  FLOW  AREA  / FACE  AREA  (GAS)  ' , 3 ( ' . '),'DSIG 

1 3(2X,G12.6,2X) ,'-') 

WRITE  (KO, 510)  DDEL 

510  FORMAT ('  FIN  THICKNESS  (GAS)  ',10('.  ')  ,'DDEL  ', 

1 3(2X,G12.6,2X) ,'FT') 

WRI TE (KO, 51  1 ) DSL 

511  FORMAT ( ' FIN  LENGTH  (GAS)  ' , 1 1 ( ' . '),'DSL  ', 

1 3 (2X,G12 .6, 2X) , 'FT') 

WRITE  (KO, 516)  DDH 

516  FORMAT  ('  HYDRAULIC  DIAMETFR  (LIQUID)  ',6('.  '),'DDH  ", 

1 3(2X,G12.6,?X),'FT') 

WRITE(KO,514)  DALF 

514  FORMAT  ( ' SURFACE  AREA  / VOLUME  (LIQUID)  " , 4 ( ' . '),'DALF 

1 3(2x, 612. 6, 2X) ,'1/FT') 

WRITE (KO, 550)  DAFAL 

550  FORMAT('  FIN  AREA  / SURFACE  AREA  (LIQUID)  ' , 7 ( ' . ' ) , 'DAFAL', 

13(2x,G12.6,2X) ,'-') 
wRITE(KO,515)  D^I 

515  FORMAT ('  FREE  FLOW  AREA  / FACE  AREA  (LIQUID)  ',2('.  '), 

1 ' D S I ',3(2V,G12.6,2X)»'-') 

WRITfc(KO,512)  DA 


P 
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r 


512  FORMAT ( ' FIN  THICKNESS  (LIQUID)  ',«('.  '),'DA  ', 

1 3 ( 2 y ,G1 2 .6, 2X ) , 'FT') 

WRITE (KO, 551  ) D^LL 

551  FORMAT ( ' FIN  LENGTh  (LIQUID)  ',1C('.  '),'DSLL  ', 

13<2X,G12.6,2X)  ,'FT') 

WRITE  ( K 0 , 5 5 2 ) a c- 

552  FORMAT ( ' PAPTING  PLATE  THICKNESS  ',?('.  '),'AE  ', 

13(?X,G12.6,2X) ,'FT') 


WR I TE  CkO , 54 1 ) NP 

541  FORMAT!'  NUMBER  OF  LIQUID  PASSES  ',?('.  '),'NP 

1 3(5X,I4,7X) 

WP I TE (KO , 51 3)  DKh 

513  FORMAT ('  FIN  THERMAL  CONDUCTIVITY  ',7('.  '),'DKH 

1 3(2X,G12.6,?X) ,'BTU/HR-FT-DEG  R') 

wRITt(KO,517) 

517  FCRMAT('CHfcAT  INPUT  RATES', 39X, 'VALUE  UNITS') 

WRITfc  ( K 0 , 5 1 5 ) GP 

518  FORMAT ( ' PUMP  ',17('.  '),'QP  ', 

1 2X,G1?.6,2X,'°TU/SEC') 

WPITE(KO,519)  Cl  C M 

519  FORMAT('  MIRROR  CONTROL  ' , 1 2 ( ' . '),'QCM  ', 

1 2X  , G 1 2 . 6 , ?X , 'PTU/ SE  C ' ) 

WRITE  (KO,  52*^)  GM 

5 20  FORMAT!'  MIPROR  ' , 1 6 ( ' . '),'QM  ', 

1 2X,G12.6,?X,'PTU/SEC') 

WRITE (KO, 521 ) Q V 1 

521  FORMAT!'  VACUUM  PUMP  1 ',13!'.  '),'QVl  ', 

1 2X,G12.6,2X,'BTU/SEC') 

WRITE (KO, 52?)  Q V 2 

522  FORMAT!'  VACUUM  PUMP  2 ',13!'.  '),'QV2  ', 

1 2X,G12.t,?X,'FTU/SEC') 

WRITE (KO,  523  ) Q E 

523  FORMAT!'  F-BEAM  ',16!'.  '),'QE  ', 

1 2X , G 1 2 . 6 , 2X , 'PTU/S  E C ' ) 

WRITE (KO,  524  ) QS 

524  FORMAT!'  SUSTAINER  ',15!'.  '),'QS 

1 2X ,G1 2 .6 , ?X , 'BTU/SEC ') 


TO  STATION  3 


WRITE(KO,525) 

525  FORMAT('CDIFFUSFR  LOSS  COEFFICIENTS') 
WRITE(KO,  526  ) C F 2 3 

WPI TE (KO, 527)  CF34 

526  FORMAT!'  STATION 

1 2X,G12.6,?X,'-') 

wRITE(K0,544)  CF111 

WRITE(KO,531  ) CF910 
WPI TE (KO,  5 30  ) C F 3 9 
WPI TE (KO, 529)  CF78 
WRITE(KO,52p)  C F 4 5 


, 7 ( 


) , 'C  F 2 ^ ', 


527 

FORMAT  ( ' 

STATION 

3 

TO 

STATION 

4 

') 

'C  F34 

» 

1 

2 X , G 1 2 . 6 , 2 X , '-' ) 

528 

FORMAT ( ' 

STATION 

4 

TO 

STATION 

5 

') 

' C F 4 5 

t 

1 

2X,612.<‘,2X,'-') 

529 

FORMAT  ( ' 

STATION 

7 

TO 

STATION 

8 

',7('. 

') 

'C  F7? 

t 

1 

2X,Gl?.<5,2x,'-') 

530 

F OPMAT ( ' 

ST  AT  ION 

E 

TO 

STATION 

9 

* * 7 ( ' . 

') 

' C F 8 9 

* 

t 

1 

2X,G12.6,2x,'-') 

5 31 

F 0 R M A T ( ' 

STATION 

9 

TO 

STATION 

10 

',?('. 

') 

' C F 9 1 

^ t 
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1 2x,G12.6,2x,'-') 

544  FORMAT  ( ' STATION  11  TO  STATION  1 ',7('.  '),'CF111', 

1 2X,G12.6,2X,'-'> 

W R I T E ( K 0 , 5 3 2 ) 

532  FCRMAT('CAMPIENT  CONDITIONS,  DECISION  VARIABLES,  AND  CONVERGENCE', 

1 ' TOLERANCES') 

WRITE (KO, 533)  T A 1 F 

533  FOPMATC'  AMBIENT  TEMPERATURE  ' , 1 0 C ' . ')  , ' T A 1 F ' , 2 X , G 1 2 . 6 , 2 X , 

1 'DEG  F') 

WRITE  (KO,  534  ) PAT 

534  FORMAT ( ' AMBIENT  PRESSURE  ',11('.  '),'PAT  ' , 2 X , G 1 2 . 6 , 2 X , 

1 'PS  I A ' ) 

WRITE(KO,535)  RHP 

535  FORMAT ('  RELATIVE  HUMIDITY  ' , 1 1 C ' . '),'RHP  ' , 2 X , G 1 2 . 6 , 2 X , 

1 'PERCENT') 

WRITE(kO,542)  FG 

542  FORMAT ('  FRACTION  OF  GLYCOL  IN  COOLING  LIQUID  . FG  ', 

1 G12.6,  ' -') 

WRITE (KO, 536)  NLIM 

536  FORMAT ( ' MAXIMUM  NUMBER  OF  LOOP  CYCLES  ',5('.  '),'\LIM  ', 

1 5X,I3,8X,'-') 

WP I TE (KO, 537)  PLIM 

537  FORMAT ('  PRESSURE  TOLERANCE  ' , 1 0 ( ' . '),'PLIM  ' , 2 X , G 1 2 . 6 , 2 X , 

1 'PSIA') 

WRITE (KO,538)  DELM 

538  FORMAT  ( ' MACH  NUMBER  TOLERANCE  ',<?('.  '),'DEL*  ' , 2 X , G 1 2 . 6 , 2 X , 

1 '-') 

WRITE (KO, 539)  ROLIM 

5 39  F ORMA  T ( ' DENSITY  TOLERANCE  ',11('.  ' ) , ' POL IM ' , 2X , G 1 2 . 6 , 2 X , 

1 'Lb/FT**3') 

WRITE(KO,540)  C ON  V 

540  FORMAT ( ' STATE  VARIABLE  TOLERANCES  ',7('.  '),'CONV  ', 

1 5(2X,G12.6,1X),2(/50X,5(2X,G12.6,1X))/50X,2X,G12.6,2X, 

2 'DEG  R') 

WRITE (KO,1QO) 

100  FORMAT ('1') 

WRI TE (KO, 1 01 ) TITLE 

101  F0RMAT(21X,14A6) 


call 

GASP(TIG, 

FN, 

GAMIG, 

2 , 

z. 

2 , 

GMU  ) 

CALL 

GASP(T2G, 

FN, 

G AM2G  , 

2 , 

Z . 

2 , 

GMU) 

CALL 

GASP(T7G, 

FN  , 

GAM3G  , 

2 , 

Z , 

Z , 

GMU) 

CALL 

G A S P ( T 4 G , 

FN  , 

G A M 4 G , 

2 , 

z , 

Z , 

GMU) 

CALL 

GASP(T5G, 

FN, 

G AM  5 G , 

2 , 

z , 

z , 

GMU) 

call 

GASP (T6G , 

FN, 

GAM6G  , 

2 , 

Z , 

z , 

GMU  ) 

CALL 

GASP(T76 , 

FN, 

GAM7G  , 

2, 

Z , 

z , 

GMU) 

CALL 

GASP ( TAG , 

FN 

, G AM  <JG  , 

2, 

z , 

z , 

GMU  ) 

call 

GA$P(T9G, 

FN 

, G AM  9G  , 

2, 

z * 

Z , 

GMU) 

CALL 

GASP(T10G 

, FN 

, gamiog 

t 

2 , 

z , 

2 , GM 

CALL 

GASP(T11G 

, FN 

, G A M 1 1 G 

f 

2 , 

z , 

2 , GM 

C0N=32. 174*1545. 43/GMU 
VG(1)*FM1*DSQRT(GAMIG*C0N*TI6) 
VG(2)=FM2*DSQRT(GAM2G*CON*T2G) 
VO(3)=FM3*DSQRT(GAM‘>G*CON*T3G) 
VG (4)=FM4*DSQPT(GAM4G*C0N*T4G) 
VG(5)=FM5*DSQRT(GAM5G*CON*T5g) 
VG(6)=FM6*DSQRT(GAM6G*C0N*T6G) 
VC.  (7)  = FM7*DSQRT(GAM76*C0N*T7u) 
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VG(8)=FM6*DSQRT(GAMRG*C0N*T?G) 
VG(9)=FM9*DSQRT(GAM9G*C0N*T9G) 
VG(1O)=FMlO*DSQPT(GAMlCG*CON*T10G) 
VG ( 1 1 ) =FM1 1 *D$QPT (G  AMI  1 G*C0N*T1 1 G) 


DO  28  K = 1 ,11 
2t  QG ( K ) =VG (K ) * A (K ) *60 . 


WRI  TE  (K.O,  102)  ALT 

102  FORMAT ( '0  ALTITUDE  ',20('.  ' ) , 1 X , G1  2 . 6 , ' (K  FT)') 

WRITE (K0,103)  PAT 

103  FORMAT ('  AMRIFNT  PRESSURE  ' , 1 6 ( ' . '),1X,G12.6,'  ( P S I A ) ' ) 

WR I TE (K  0 , 1 U4  ) TA1F 

109  FORMAT('  AMR1FNT  TEMPERATURE  ',15('.  ' ) , 1 X , G 1 2 . 6 , ' (DFG  F)') 

WP I TE  (KO  , 1 05 ) RHP 

105  FORMAT ( ' RELATIVE  HUMIDITY  ' , 1 6 ( ' . '),1X,G12.6,'  (PERCENT)') 
WRITE (KO, 107)  P I 6 A 

107  FORMAT ('  L A r E P INLET  PRESSURE  ',14('.  ' ) , 1 X , G 1 2 . 6 , ' (ATM)') 


TIGM=TIG/1 .P 
WRITE (K0,10p) 

1 Cfc  FORMAT ('  LASER  INLET  TEMPERATURE  ' , 1 3 ( ' . '),1X,G12.6,'  (DFG  K ) ' ) 
2= (T02G-T01G) /I  .5 
WR I TE (KO , 1 42 ) Z 

142  FORMAT ( ' TEMPFRATURE  RISE  ACROSS  CAVITY  ',9('.  '),1X,G12.6, 

1 ' (DcG  K ) ' ) 

WRITE ( K 0 , 1 U 9 ) VIM 

1C9  FORMAT ('  LASER  INLET  VELOCITY  ',14('.  ' ) , 1 X , G 1 2 . 6 , ' (M/SEC)') 

WR I TE (HO, 1 1 0)  PLKW 

110  FORMAT ( ' LASER  OUTPUT  POWER  . . ',1 3('.  ' ) , 1 X , G 1 2 . 6 , ' (KW)') 

WPITE(K0,141)  EFFLP 

141  FORMAT ( ' LASER  EFFICIENCY  ' , 1 6 ( ' . '),1X,G12.6,'  (PERCENT)') 

WRITfc(K0,114)  WG 

114  FORMAT ('  LASFR  GAS  MASS  E L 0 * RATE  ',12('.  '),1X,G12.6, 

1 ' (LE/SEC ) ' ) 

WRI TE (KO,  1 1 5 ) Q 6 G 

115  FORMAT ( ' VOLUME  FLOW  RATF  AT  BLOWEP  INLET  ',t('.  ' ) , 1 X , G 1 2 . 6 , 

1 ' ( C F*)  ' ) 

WRITt(KOfllfe)  OCORR 

116  FORMAT ( ' CORRECTED  VOLUMF  FLOW  RATE  AT  BLOWER  INLET  ',3('.  '), 

1 1 X , G 1 2 . 6 , ' ( C F M ) ' ) 

WR I TE (KO , 1 39 ) WA 


1 39 

FORMAT  ( 

' COOLING 

AIR  FLOW 

RATE 

' , 1 4 ( ' . 

') 

, 1 X 

, G 1 2 » 6 , 

1 

' ( LI  /SF  C ) ' ) 

WRITE  (KO 

,140)  C W A 

140 

FORMAT  ( 

' COOLING 

AIR  FLOW 

RATE 

' , 1 4 ( ' . 

') 

, 1 X 

,G1  2.6, 

1 

' ( C F M ) 

') 

W P I T E ( K 0 

,117)  BPR 

1 1 7 

F ORM  AT  ( 

' BLOWER 

PRESSURE 

RATIO 

' , 1 4 ( ' . 

') 

, 1 X 

, G 1 2 . 6 , ' (-)' 

w P I T t ( K 0 

,111)  Ei  L P R M 

, P L P R 

1 1 1 

F 0 P M A T ( 

' BLOWER 

POWER  REQUIRED 

' , 1 4 ( ' . 

') 

, 1 X 

,612.6,'  (KW) 

1 EX, G 12. 6, IX,'  (HP)') 

WPIT£(K0,112)  PPRM,  PPR 

112  FORMAT ( ' COOLING  PUMP  POWER  REQUIRED  ',11('.  '),1X,G12.6, 
1 ' <Kw)',8X,G12.6,1X,'  (HP)') 


WRITE(K0,117)  PFM,  PE 

113  FORMAT ('  COOLING  FAN  POWER  REQUIRED  ' , 1 1 ( ' . '),1X,Ul2.6, 

1 ' (KW)',8X,G12.6,1X,'  (HP)') 

W R I T £ ( K 0,106)  F N 

106  FORMAT ('  LASER  GAS  MIXTURE  ',1d('.  ' ) , 1 X , F 5 . 2 , 3 ( ' TO ' , F 5 . 2 ) , 2 X , 


: 


1 


117 


r 


1 '(HE  TO  N2  TO  C 02  TO  CO)') 

WRI TE (KO, 143)  GMU 

143  F OR  MAT ( ' MOLECULAR  WEIGHT  OF  GAS  ',1 2('.  '),1X,G12.6, 

1 ' (-)') 

WR  I TE<  icO, 1 1 ® ) 

1 1 £ FOR MAT('G',PX, 'HEAT  cXCHANGEK  VARIAPLES',15x, 'EXCHANGER  1 ' , 6 X , 

1 'EXChANGER  2 ' , 7X , 'ft  A D I ATOP ' , 6X , 'UNITS' ) 

WRI TE (KO, 1 1 9)  GDI,  Q 0 I I , GDA 

119  FCRMAT('  HEAT  TRANSFtR  RATE  ' , 1 1 ( ' . ' ) , 2 X , 3 ( G 1 2 . 6 , S x ) , 

1 ' bTU/HP') 

WR 1 T E (K 0 , 1 2°)  DPI,  DPII,  DP 

120  FORMAT('  LIQUID  PRESSUftF  DROP  ' , 1 0 ( ' - ' ) , 2 X , 3 ( G 1 2 . 6 , 5 X ) , 

1 ' PSI A') 

WPITE(KO,121)  NP 

121  FORMAT  ('  NUMPFR  OF  LIQUID  PASSES  ' , 9 C ' • ' ) , 2 X , 3 ( 6 X , I 2 , 9 X ) , ' -') 
N F U L L = a 

IF  ( f,  FULL  . N E . 0)  RETURN 
WRITE  (K 0,122)  T y I , TMI  I , T M 

122  FORMATE'  MEAN  TEMPERATURE  ',9('.')  , '(LIQUID)', 4('.  '),2X, 

1 3(012. 6, 5X) ,'  DEG  ft ' ) 

WPITE(K0,12T)  TGMI,  T G M I I , T M A 

123  FORMAT(32x,'(GAS)',1Gx,3(G12.!S,?x)) 

WRITE(KO,124)  ZIMU,  Z I I M U , ZNU 

1 2 C FORMAT ('  VISCOSITY  ',16  ('.')  , '(LIQUID) ',4('.  ' ) , 2 X , 3 ( G 1 2 . 6 , 5 X ) , 

1 ' LB/FT-SEC') 

WRITE (KO, 1 23)  GIMU,  GIIMU,  AMU 
WRITE(KO,125)  PPI,  PRII,  PR 

125  FORMAT ( ' PRANDTL  NUMBER  ',11  ('.')  , '(LIQUID)', 4('.  '),2X, 

1 3(612. 6, 5X),'  -') 

WR I TE (KO , 1 2*  ) PPG!  , PR  6 I I , PRA 
WPITE(k0,126)  CPI,  CPU,  CP 

126  FORMAT ( ' SPECIFIC  HEAT  ',12('.')  , '(LIQUID)', 4('.  '),2X, 

1 3 (G1  2 .6 , 5 X ) , ' ElTU/LB-DEG  R') 

WRITE (K0,12T)  CPGI,  CPGII,  CPA 
WRITE(k0,12T)  G I , G I I , G 

127  FORMAT ('  *ASS  FLUX  ',16('.')  , '(LIQUID)', 4('.  '),2X, 

1 3(G12.6,5X) ,'  LS/FT**2-SEC ') 

WRITE(K0,123)  GGI,  GGII,  GA 
WRITE (K0,128)  PEI,  RE  I I , RE 

128  FORMAT  ('  REYNOLDS  NUMBER  ',1Cl('.'),'(LICUID)',4('.  '),2X, 

1 3(012. 6,5X) ,'  -') 


129 


130 


1*1 


132 


WRITE  ( K 0 , 1 2 3 ) REGI,  k E G I I , RE» 
WRITE(K0,129)  FI,  F 1 1 , F 
FORMAT ('  FRICTION  FACTOR  ' , 1 0 ( ' 
3(G12.6,5X)  ,'  -') 

WPITE(K0,123)  G F I f 0,  F I I , AF 
WR I Tt (KO , 1 30)  ZJI,  Z J I I , ZJ 
FORMAT('  heat  transfer  FACTOR  ' 
3(612. 6, 5X)  ,'  -') 

WRITE! K 0,123)  GJI,  G J 1 1 , Aj 
WRITE ( K 0 , 1 3 1 ) HI  , HII  , H 
FORMAT ('  HEAT  TRANSFER  COEF.  ', 
3(012. 6, 5X),'  BTU/FT  * * 2 - S E C~DEG 
WRITE ( K 0 , 1 2 3 ) H C I , H 0 I I , H« 
wPITE(fO,132)  C I , C 1 1 , C 
FORMAT ( ' FLUID  CAPACITY  ft  A T F ', 
3 (G  1 2 .6 , 5X  ) , ' EiTU/SFC-DEG  R') 


) , 


, 5 ( ' 


6 ( ' . 

ft') 


6 ( 


(LIQUID) ' , 4 ( ' . ' 


')  , '(LIQUID) ' , 4 ( 


), '(LIQUID)', 4(' 


), '(LIQUID) ',4(' 


2 X , 


') ,2X, 


'),?X, 


') , ?X, 
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WRI TE (KO, 1 23)  CGI,  CGII,  CA 
WRITE (KO, 1 33)  UI,  UII,  U 

133  FORMAT  ( ' OVERALL  HEAT  TRANSFER  COEF.  ' ) , 2 X , 3 ( 6 1 2 . 6 , 5 X ) , 

1 ' bTu/FT**2-SEC-DEG  R') 

WRITE (KO, 1 34)  AG  I , A GI I , AA 

134  FORMAT ( ' GAS  HEAT  TRANSFER  SURFACE  AREA  ' , 5 ( ' . ' ) , 2 X , 

1 3(G12.6, 5X) , ' FT**2') 

WRI  TE  (KO  , 1 35  ) TUNI,  TUMI,  TUN 

135  FORMAT!'  NUMPER  OF  TRANSFER  UNITS  ',?('.  '),2X, 

1 3(612. 6, 5X),'  - ' ) 

wR  I TE  (KO  , 1 36 ) TO W I , T 0 » I I , TOW 

136  FORMAT!'  CAP  C A * M A ',16('.  '),2X,3(G12.6,5X),'  -') 

WRITE (KO, 137)  ESSI,  ESSII,  ESS 

137  FORMAT ( ' SINGLE  PASS  EFFECTIVENESS  ',«('.  '),2X, 

1 3(612. 6, 5X) ,'  -') 

WRITE (KO,13R)  EPSI,  LPSII,  EPSA 

138  FORMAT  ('  TOTAL  EFFECTIVENESS  ' , 1 1 ( ' . '),2X, 

1 3 ( G 1 c . 6 , 5 X ) ,'  -') 

WRITE (KO.IOO) 

WRITE (KO.IOI  ) TITLE 
WRITE (K  0 , 200 ) 

20C  FORMAT('0',16X,'GAS  FLOW  SYSTEM  VAPIAPLFS  BY  STATION') 
WRIT£(K0,2u1 ) 

201  FORMAT ('0  STATION  N UMB E R ' , 1 6 X , ' 1 ' , 1 3 X , ' 2 ' , 1 3 X , ' 3 ' , 1 3 X , ' 4 ' , 

1 13X,'5',13X,'6'/37X,'7',13x,'P',13x,'9',12x,'10',12Xt'11') 

WRITE (K0,2U4)  P01G,  P02G,  P03G,  P04G,  P05G,  P06G,  P07G,  P08G, 

1 P09G,  P010G,  P011G 

204  FORMAT ('0  TOTAL  PRESSURE  (PSIA)', 

1 6(1X,G12.6,1X)/31X,5(1X,G12.6,1X)) 

WRITE (KO,  203  ) PIG,  P 2 G , P 3 G , D4  G , P 5 G , P6G,  P 7 G , P8G,  P9G,  P10G, 
1 PI  1 G 

203  FORMAT ('0  STATIC  PRESSURE  (PSIA)', 

1 6(1X,G12.6,1X)/31X,5(1X,G12.6,1X)> 

WRITE(KO,206)  T 0 1 G , T02G,  T03G,  T04G,  T05G,  T06G,  T07G,  T08G  , 

1 T09G,  T010G,  TC11G 

206  F0RMAT('0  TOTAL  TEMP.  (DEG  R)', 

1 6(1X,G12.6,1X)/31X,5(1X,G12.6,1X)) 

WRITE (KO,  205  ) TIG,  T 2 G , T 36  , T 4 G , T 5 G , T 6 G , T76,  T8G , T9G,  T10G, 
1 T 1 1 G 

205  F ORMA  T ( '0  STATIC  TEMPERATURE  (DEG  P)', 

1 6(1X,G12.6,1X)/31X,5(1X,G12.6,1X)) 

WRITE(K0,207)  FM1,  FM2,  F m 3 , F M 4 , F M 5 , FM6,  F M 7 , FMP,  FM9,  FMIIC, 
1 F M 1 1 

207  FORMAT ( '0  MACH  NUMBER  (-)', 

1 6(1X,G12.6,1X)/31X,5(1X,G12.6,1X)) 

WRITE (KO, 20°)  VG 

208  FORMAT ('0  VELOCITY  (FT/SEC)', 

1 6(1X,G12.6,1X)/31X,5(lX,Gl2.fc,1X)) 

WRITE (K0.2U9)  QG 

209  FORMAT ('0  VOL.  FLOW  PATE  ( C F M ) ' , 

1 6(1X,G12.6,1X)/31X,5(1X,G12.6,1X)) 

Z=GMU*. 093178 
R 06 ( 1 )=PIG*Z/TIG 
R OG ( 2)  = P2G*Z/T2G 
ROG ( 3)=P3G*Z/T3G 
R OG ( 4)=P4G*Z/T4G 
R OG ( 5)=P5G*Z/T5G 
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r 


RCG(  6)=P6G*Z/T6G 
R OG ( 7)=P7G*Z/T7G 
R OG ( 8 ) =P86*Z /TpG 
R OG  ( 9)=P9G*Z/T9& 

ROG ( 1 0) =P1 GG*Z /T1 OG 
ROG ( 1 1 ) =P1 1G*Z/T1  1G 
WRITE (KO,220)  ROG 

220  F0RMAT('0  DENSITY  (lS/FT**3)', 

1 6(1X,G12.6,1X)/31X,5(1X,G12.6,1X)) 

UP  I TE (KO , 2 1 2 ) A 

212  FORMAT('0  A R F A (FT**?)', 

1 6(1X,G12.6,1X)/31X,5(1X,G12.6,1X)) 

WRITE(KO,213)  GAMIG,  6AM2G,  G A M ?G , GAM4G,  GAM5G,  GAM6G,  GAM7G, 

1 GAM PC,  G A M 9G  , G A M 1 GG  , G AM  1 1 G 

213  FORMAT ('0  GAMMA  (-)', 

1 6(1x,G12.6,1X)/31x,5(1X,Gl2.6,1X)) 

P D ( 1 ) = (PC1G-P02G) *27.7 

PD(?)=(P0?G-P03f )*27.7 

PD(3)=(P03G-P04G)*27.7 

PD(4)=(P04G-P05G)*27.7 

PD(5)=(P056-P06G)*27.7 

PD(6)=(P06G-P07G)*27.7 

PD(7)=(P07G-P0RG)*27.7 

PC  (8  ) = ( P08G-P09G ) *27. 7 

PD(9)=(P09G-P010G)*27.7 

PD(10)=(P010G-P011G)*27.7 

PD(11)=(P011G-P01G)*27,7 

DO  27  K =1  , 1 1 

27  PP(K)=FD(K)/(?7.7*P06G) 

WPITE (K0,?10)  PD 

210  FCRMAT('0  PRESSURE  DROP  (IN  H20)',6X, 

1 6<1X,G12.6,1X)/37X,5(1X,G12.6,1X)) 

WRITE (KO, 21 1 ) PP 

211  FORMAT ( '0  PRFSSURE  DPOP  / P06  (-)',6X, 

1 6(1X,G12.6,1X)/37X,5(1X,G12.6,1X)) 

wPITE(kO,214) 

214  FORMAT ( 'O', 1 2X , 'COOL ING  AIR  VARIABLES') 

WRI TE (KO, 21 5) 

215  F0RMAT('0  STATION  N U M B E R ' , 1 8 X , ' 1 A ' , 1 2 X , ' 2 A ' ) 

WRITE (KO, 21 4)  PAT,  PA? 

216  F0PMAT('C  PRFSSUPt ' , 1 2X , '(PS  I A ) ',2 (1 X ,G1 2.6 , 1 X)  ) 

WPITE  ( K 0 , 2 1 7 ) T A 1 , T A 2 

217  FORMAT ( '0  TEMPER ATURE' ,8X ,'( DEG  R ) ' , 2 ( 1 X , G 1 2 . 6 , 1 X ) ) 

WR I TE ( KO  , 2 1 8 ) R 0 A 1 , R 0 A 2 

218  FORMAT('0  DENSITY', 9 X , ' ( L b / F T * * 3 > ' , 2 ( 1 X , G 1 2 . 6,  IX)) 

D P P A = (F  A T -P  A 2 ) *27.7 

WR1TE(K0,219)  DPPA 

219  FORMAT ( '0  PRCSSUPE  DPOP  (IN  M 2 0 ) ' , 7 x , G 1 2 . 6 ) 

WRI T£ (KO, 1 00) 

WRI TE (KO, 1 01 ) TITLE 
WRITE(KO,300) 

300  FORMAT ( 'C' , 16X , 'COOL INu  SySTl*  VARIAELES  PY  STATION') 

WRI TE (KO, 3C1 ) 

301  FORMAT ( 'G  STATION  N o * b F R ' , 1 9 X , ' 1 ' , 1 3 X , ' 2 ' , 1 3 X , ' 3 ' , 1 3 X , ' 4 ' , 1 3 X , 

1 '5',13x,'6',17x,'7'/3  8x,'8',1‘,x,'9',12x,'10',12x,'11',1?x, 

2 '12',12X,'13',12X,'14',/37X,'15',12X,'16',1?X,'17',12X,'1°', 

3 12X,'19',12X,'20') 
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WRITE (k0,302)  P 

302  FORMAT ( '0  PRESSURE  (PSIA)', 

1 7<1X,G12.6,1X)/31X,7(1X,G12.6,1X)/31X,6(1X,G12.6,1X)> 

WRITE (KO, 303)  T 1 , T 2,  T3,  T4,  T5,  T6,  T 7 , T°,  T9,  T10,  Til,  T1 
1 T 1 3 , T 1 4 , T 1 5 , T 1 6 , T17,  TIP,  T19,  T20 

3G3  FORMAT  (■”0  TEMPERATURE  (DEG  R)', 

1 7(lX,G12.6,lX)/3lX,7(lX,G12.fc,lX)/3lX,6(1X,G12.6,lX>) 

WRITE (KO, 304)  VC 

304  F ORMAT ( '0  VELOCITY  (FT/SEC)', 

1 7(1X,G1?.*,1X)/31X,7(1X,G12.6,1X)/31X,6(1X,G1?.6,1X)) 

WRITE  ( K 0 , 3 0 5 ) QC 

305  FORMAT ('0  VOL.  FLOW  PATE  ( F T * * 3 / S E C ) ' , 

1 7(1X,G12.6,1X)/31X,7(1X,G12.6,1X)/31X,6(1Y,G1?.6,1X)) 

WPITE(kO,306)  QCG 

306  FORMATt'O  VOL.  FLOW  RATE  (GAL/MIN)  ' , 

1 7(1X,G12.6,1X)/31X,7(1X,G12.6,1X)/31X,6(1Y,G12.6,1X)) 

WR I TE (KO, 307)  DPC 
WRITE (KO, 309)  AC 

309  FORMAT('0  AREA  (FT**?)', 

1 7(1X,G12.6,1X)/31X,7(1X,G12.6,1X)/31X,6(1X,G12.6,1X)) 

307  FORMAT ('0  PRFSSURE  DROP  (PSIA)', 

1 7(1X,G12.6,1X)/31x,7(lX,G12.o,lX)/31x,6(1X,G12.f,lX)) 

WRITE(KO,30C)  D p C R 

3 Ofc  FORMAT ('0  PRESSURE  DROP  / Pi  (-)', 

1 7(1X,G12.6,1X)/31X,7(1X,G12.6,1X)/31X,6(1X,G12.6,1X)) 

RETURN 

END 
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Table  1 

LASER  TEST  CASE 
DECEMBER  1976 

INPUT  LIST 


CAVITY  * LASER  SPECIFICATIONS 


DESCRIPTION 

NAME 

VALUE 

UNITS 

INLET  PRESSURE  . . . 

• . • • 

• 

PIGA 

.100000*001 

ATM 

INLET  HEIGHT  .... 

• • • . 

• 

Hi  “ 

.500000*001 

CM 

WIDTH  IN  FLOW  DIRECTION  . . . 

• 

W T 1 M 

. 5C0C00  + GP1 

CM 

LENGTH 

• . • • 

CLlM 

.4G00C0+C02 

CM 

INLET  HYDRAULIC  DIAMRTER  . . 

D C M 

. °6?600*001 

CM 

LASER  OUTPUT  POWER 

• * • * 

PLKW 

.500000* 

002 

K w 

LASER  EFFICIENCY  . . 

• • • • 

efflp 

.100000*002 

PERCFNT 

DESIGN  OPTION  .... 

• • • . 

NDSGN 

2 

- 

GAS  WIXTURE  

.... 

FN 

J.00  2 

.00 

.60 

50  HE 

TO  F.2 

TO 

CG2 

TO  CO 

INITIAL  ESTIMATE  OF  STEADY 

STATE  QUANTITIES 

I F 

r . , us  es 

GAS  STATIC  TEMPERATURE 

- STATION  1 

ETIG 

.000000 

DEG 

P 

TA  1 

♦ 

20. 

GAS  TOTAL  TEMPERATURE 

- STATION  2 

ET02G 

.000000 

DEG 

R 

T A 1 

♦ 

210. 

GAS  TOTAL  TEMPERATURE 

- STATION  6 

ET06G 

.000000 

DEG 

P 

T»1 

♦ 

JC  . 

GAS  TOTAL  TEMPERATURE 

- STATION  7 

ET07G 

.ocronc 

DEG 

P 

T » 1 

♦ 

70. 

coolant  TEMPERATURE  - 

STATION 

1 

F T 1 

.OOOOOC 

DEG 

r 

T A 1 

♦ 

20. 

COOLANT  TEMPERATURE  - 

station 

2 

« 

E T 2 

.ncnooo 

DEG 

p 

T A 1 

♦ 

20. 

COOLANT  TEMPERATURE  - 

STATION 

3 

• 

• 

ET  J 

.000000 

DEC- 

p 

T A 1 

♦ 

<.0. 

COOLANT  TEMPERATURE  - 

STATION 

u 

. 

• 

E T 4 

•OCOCOC 

DEG 

p 

T A 1 

♦ 

55  . 

COOLANT  TEMPERATURE  - 

STATION 

6 

. 

• 

E T 6 

.000000 

DEG 

p 

T A 1 

♦ 

JO. 

COOLANT  TEMPERATURE  - 

STATION 

• 

FT' 

. "00000 

DEG 

p 

T A 1 

♦ 

6 5. 

COOLANT  TEMPERATURE  - 

STATION 

10 

E T 1 0 

•OOOOOC 

DEG 

p 

T AT 

♦ 

65. 

COOLANT  TEMPERATURE  - 

STATION 

1 U 

E T 1 4 

."00000 

DEG, 

p 

T A 1 

♦ 

i n 

COOLANT  TEMPERATURE  - 

ST  AT  I ON 

1 6 

ETIo 

•PCOOOO 
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APPENDIX  B 
(David  Washington) 

Phase  Diagram  for  a Mixture  of  (a)  Carbon  Dioxide  and  (b)  Nitrogen, 
at  (P  = 1 atm)  Pressure 


Condition  for  Phase  Equilibrium 


2 


f 


A 


or,  at  equilibrium  the  fugacity  of  a given  component  is  the  same  in  each 
phase.  In  general , 


f ^ = 7^  = fjj  = ...  for  all  phases 

f g = 7^  = f g = . . . for  all  phases 


for  all  components. 

Pure  Carbon  Dioxode  at  1 atm  pressure  will  exist  as  a vapor  or  solid, 
depending  on  the  temperature  (See  P-T  diagram).  Pure  Nitrogen  at  1 atm 
pressure  can  exist  in  all  three  phases  (vapor,  liquid,  solid),  depending 
on  temperature  (See  P-T  diagram).  However,  the  fugacity  of  compressed 
liquids  and  solids  can  be  treated  alike,  as  follows. 

For  a pure  substance  at  constant  temperature, 

R T (d  In  f)T  = v d Py 

integrating  between  saturation  state  on  Pressure  P, 
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In  f)T  = f v d PT 


R T In  a*  v (P  - PSAT) 


Assume  constant 
specific  volume  (v) 
for  liquids  and 
solids. 


f 

f$AT 


v (P-P' 
RT 


Since  v is  small  for  liquids  and  solids,  the  quantity  v (P  - P3A^)  is 
small  for  moderate  changes  of  pressure.  Therefore, 


e = 1 

1 iquid 
and 
solids 


fL  _ fSAT 


fS  fSAT 


For  Nitrogen  at. 


T = 200UR 
pSAT  _ 226.853 

in^ 

v = .02613  ft3/Lbm 


R + 55.15 


ft-Lbf 

Lbm°R 


P = 1 atm 


r 


3ST 


e 


-.072 


.93 


Therefore, 


f « f 


SAT 


with  It  error  at  this  check  point. 


From  the  P-T  chart  for  Nitrogen,  it  is  observed  that  only  a small 
portion  of  the  liquid  phase  is  involved  (since  1 atm  is  near  the  triple 
point).  Also,  since  the  fugacity  of  liquids  and  solids  are  calculated 
alike  (fL  = fS<=.  fSAT)  it  is  felt  that  the  Nitrogen  component  can  be 
considered  in  two  phases  (solid  and  vapor)  for  this  problem.  Consequently, 
for  this  problem  (1  atm),  consider  a Binary  Solution  (vapor  and  solid) 


of  two  components  (A-  Carbon  Dioxide,  and  B-Nitrogen).  Therefore, 

Condition  for  Equilibrium 


f ^ = f v 

tn2  tn2 


Is  = f v 
tC02  tC02 


Additional  Assumptions 

1)  Assume  Ideal  Solution. 


7.*  « y? 

l J-\ 


r. 


where 


4> 

i 


fl1 


fT 


= phase 
= component 

= fugacity  of  component  (i)  in  phase  (<j>) 

= fugacity  of  Pure  (i)  in  the  same  phase  as  the 
mixture  on  the  same  temperature  and  pressure 

= mole  fraction  of  ( i ) 


7 fs 

tA 


ZB  fB 


yA  fA 

*B  fB 
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j 


or 


where. 


mole  fraction  of  (A)  in  Solid  Phase 
mole  fraction  of  (B)  in  Solid  Phase 
mole  fraction  of  (A)  in  Vapor  Phase 
mole  fraction  of  (B)  in  Vapor  Phase 


2)  Assume, 


= f 


SAT 


As  discussed  earlier 


B 


= f, 


SAT 


B 


3)  Assume  that  pure  saturated  vapor  at  T and  PSAT  behave  as  an 
ideal  gas.  Therefore, 


,SAT 


,SAT 


f A = Pft  (Vapor  Pressure)A 

fgA^  = PgAT  (Vapor  Pressure)^ 

4)  For  the  Vapor  Phase,  assume  that  pure  gas  (A)  and  pure  gas  (B) 
behave  as  ideal  gases  at  T and  P.  Therefore, 


'A 

„v 


= P 


= P 


Combining  (1),  (2),  (3),  and  (4), 


Z p ^AT  = y p 

kA  yA  y 


l p ~ yB  P 

^B  B D 
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Known  as  (RAOULT'S  LAW). 
Good  for  systems  at  Low 
Pressure. 


Combined  with, 


ZA+ZB  = 1 


K + yB  = 1 

gives  a system  of  four  equations  and  four  unknowns  (Z^,  Zg,  y^,  yg). 
Combining  final  equations. 


(1) 

ZA 

p SAT 
kA 

= *AP 

(2) 

ZB 

p SAT 

kb 

■ VBP 

(3) 

ZA 

+ ZB 

= 1 

(4) 

yA 

+ yB 

= 1 

SAT 


SAT 


+ Zr 


= 1 


SAT 


A 


LA  pW 
B 


+ (1  - 


ZA>  = 


,SAT 


SAT 


TAT 


- 1 


P 

sKT 


- 1 


B 


P - P 


SAT 


B 


p SAT  SAT 
KA  ' B 


SAT 


1 
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CRITICAL  POINT  AND  TRIPLE  POINT  PROPERTIES 
OF  PURE  NITROGEN  AND  CARBON  DIOXIDE 


MOLECULAR 


CRITICAL 

POINT 


TRIPLE 

POINT 


SATURATION 
DATA  AT 
P = 1 atm 


(B)  (A) 


NITROGEN  CARBON  DIOXIDE 

(n2)  (co2) 


WEIGHT 

/^q  n i o Lb  ms s s 

/i/i  m Lb  mass 

Lb  mole 

Lb  mole 

ZC 

.291 

.275 

Tc 

126.2°K 

304. 2°K 

TC 

227. 1°R 

547. 5°R 

TC 

-232 . 57°F 

87.8°F 

PC 

33.5  atm 

72.9  atm 

pc 

492.  Lbf/in2 

1071.  Lbf/in2 

ttp 

63.1 5°K  ( -210°C) 

216.55°K  (-56.6°C) 

ttp 

1 1 3. 67°R 

389. 79°R 

ttp 

-345. 6°F 

-69.88°F 

PTP 

94.01  rrnihg 

3885.1  mmhg 

PTP 

.1237  atm 

5.112  atm 

PTP 

1.818  Lbf/in2 

75.126  Lbf/in2 

tsat 

1 39 . 255°R  (L-V) 

350. 3°R 

-320.415°F  (-195.786C) 

-109. 3°F  (~78.52°C) 

77. 364°K 

194.63°K 

136 


i A G A A 


M 'xr>M  c 


J rj  ll 


C At ? tfa  w 


_____ 

L j 

i 

~~ 

A/ 1 

r f?o &B/J 

! 

i — f—  1 

Phase  Diagram  for  a Mixture  of  (A)  Nitrogen  and  (B)  helium,  at  (P  = 1 atm) 
Pressure 


Solution 


Assume  a binary  solution  (solid-vapor)  at  two  components.  Assume 


an  ideal  solution  and  using  Raoult's  Law,  as  in  Appendix  A,  we  have, 

(vapor) 


Combine  with. 


yield 


(solid) 

l p S 
‘A  kA 


Z„  P*SAT  = y.  P 


RAOULT'S  LAW 


l p SAT  p 

B yB 


ZA+ZB  = 1 


yA  + ye 


= 1 


p - p 


SAT 


B 


p "SAT  p SAT 

KA  ‘ PB 


SAT 


yA  = ZA 
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Phase  Diagram  for  a Mixture  of  (A)  Carbon  Dioxide  and  (B)  Helium  at 
(P  = 1 atm)  Pressure 


Solution 

Assume  a binary  solution  (solid-vapor)  of  two  components.  Assume 
an  ideal  solution  and  using  Raoult's  Law,  as  in  Appendix  A and  B,  we 
have, 


Combine  with, 


(solid) 

(vapor) 

Z P SAT 
kA 

' yAP 

Z P SAT 
LB  VB 

" yBP 

ZA  + ZB 

= 1 

*A  + yB 

= 1 

7s. 

p p SAT 
K B 

ZA 

tist  rs r 

kA  ' kB 

v.  = 1. 

p SAT 
kA 

RAOULT'S  LAW 
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PA**hETE»  M*  1 2 1 JU241  . iSEGrlRO.  -II=ISE6»1, 
COMMON  /Bod/  6tCP,CVtRFtR,RG tO,A,CA,CB,CC,HXt 


Biim  HOT  1J1H1  ! ffl : a£T*f  Bc  1 • "» 

PARAMETER  NZ:4 
DIMENSION  C SL  ( NZ  ) »C  SR  INZ  ) 

COMMON  /DLL/  CSL f CSR » SL , SR 
DO  50  L = l,2 

C OEFINE  C0A5TA*T5 

IF(L.EO.l)  XNNrl.O 

IF  ( L • EC  » 2 ) XNN--1.0 

SEG-I  SEG 

IC  = II  '1 

M I - 1 1 '2 

MU- JJ  '1 

BCU  = C.O 

BC21- 200 .0 

X=BC2  1-BC 1 1 

DX-X/  SE6 

OlS  = 6 .0 

CP-O.  24 

G = 1 *4 

RF  = 53  .34 

*P=  2 1 16.8 

X0:0.  C0234 

XUrlOC.O 

A Z-0 . C 

ZERO:  C.O 

HX:4.  C 

F=0.02 

XL:1*  C 

CVrCP  /G 

RR:RF  /778.U1 

R = RF*  32. 1739 

XT  = XP/XO/R 

ArSOR  T(G*R*Xt  ) 

U-XU/ A 

CA={6-1. 01/(2. 0*G> 

CB: ( G «1.G)/(2.Q*G) 

CC-2. C/IG*1.0> 

C BOUNDARY  CONDITIONS  (LEFT) 

BC1  I 1 ):  B C 1 1 
BC 1 ( 2 ): - 1 0 . 0 
BC1 1 3 >:U 
BC1 (4  1=1.0 
BC1  (5  ):  0 . 0 
BC1  (6  ): 0 • 0 
BC  1 ( 7 ): 0. 25 
BC1  ( 8 ):XP 
BC1  (9  »:X0 
BCitlfl=XT 
BCi  ( 1 l)  =0.2 

C BOUNDARY  CONDITIONS  (RIGHT) 

BC2  ( 1 )=8C2 1 

BC2<2  »:-10.Q 
BC2  ( 3 ):U 
BC2I4  ):  1 . 0 
BC2  ( 5 1-0.0 


M=1B 

AZ  » ZERO  »F «XL 
,8C2f M) , BIC (M  > 
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ramie* 

BC2<8  1=XP 
BC2  C9  1=XD 
BC2«1C1=XT 
BC2<  1 11=0.2 
C FLO*  DIRECTIONS 

BC1 ( 3 1=BC 1(3) *XNN 
BC2(  3 >=BC2(  3 1 *X  NN 
C L0*0  BoUNDARy  CONDITIONS 

00  15  J=1,MJ 
DO  1*4  K = 1,M 
PU»1  tJ  1 =BC  1 ( K 1 
P ( K r I If  J ) =BC  2 ( K > 

14  CONj  I MJE 

15  CONTINUE 

C INITIAL  CONDITIONS 
DO  1 8 I = 2 1 1C 
I H 1 = 1-1 

P(lfl  ,1>=P<1.IM1,11*DX 
18  CONTINUE 
B IC  ( 2 1=0.0 
B IC  < 3 »=U 
BIC  I 4 1=1 .0 
B 1C  ( S 1=0 . 0 
B IC  ( 6 1=0.0 
B IC ( 7 1=0.25 
B IC  ( 8 »=XP 
B IC ( 9 1=*D 
BIC  ( I C)=xT 
BIC  ( 1.11  =0.2 
C FLOW  DIRECTIONS 

8 IC ( 3 1=8IC ( 3 1 *XNN 
C LOAD  INITIAL  CONDITIONS 

DO  17  1=2, IC 
DO  16  K=2fM 
Pt*,I  fll  =-8IC  4K  1 

16  continuE 

17  CONTINUE 
C SHOCK  WAVE 

SW ( 1 I =100.0 
SW ( 2 1 = 0.0 
SW|31 =1.13 
pRES=  1.323 
DENS=  1.221 
TEHP=  1.084 
SPEED =SQRTiTE HP  1 

0 1 =B  I C(  9 1 
S 1 =8 1 C(  5 1 

T 1 =B  1 C(  1 0 1 

QQ=(2C0.0)*(9.4  7 8E-4)*(28.32  * 

Q=(00/01)*(1. 0/32. 17391 

05=0/  11 

T2=Q/CV*T1 

P2=D1  *R*T2 

A2=S0RT(G*R*T21/A 
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HfR3f*»8hi  -1.0/SWC3) > *B IC ( 3 > 


SW<5»  =BIC<3) 

SW  ( 6 ) =PRES*BIC(8  ) 

SWI7)  -B I C ( 6 ) 

SW  ( 8 ) =DENS*BICC9) 

SW191=BIC(9) 

SW( 10  )rTEMP*BIC< 10> 

SWlll  ) = B IC  ( 1 0 ) 

SW1 12  )=SPEE0*BlC (4) 

SW(  13  l=B I C 4 *»  > 

SH=CP  *AL0G  I TEMP  ) -RR*AL0G  (PR^S) 
SW 1 14  )=SH/G/RR*BIC(S) 

SW( 15  l=BIC(5  » 

SW1 16  )=B1C(6) 

SWt  17  1=B  I C C 7 » 


SW(  18  >- B I C ( 1 1 > 
SL(  1 ) =SW  < 1 )-0IS 


SRI  1 ) =SW 1 1 )40IS 
00  32  IK=2,N 

SUIN  >=SW(IK  > 

SR1IK  )-  S W ( I K ) 

32  CONTINUE 

SL13) =-SW(3) 

SL14»=SW<5)4«2.0/1G»1.0))*(SL(3)-1.0/SL(3|) 

C contact  surface 
TA=PRES-1.0 
VB=2.  C/G 
YC  = (6  •»1.0)*PRES 
YD  = G-  1.0 

YF=YA*SQRT(YB/<YC*YO)) 

YF  = 0. C 
CSL<1  )=SL<1) 

CSLI2  )=SL<2) 


300 


301 


109 


CSL(3  )=U-YF 
CSL  ( 4 )=-l  .0 
C SR  ( 1 >=SRC  1 > 

C SR  ( 2 )r  SR  ( 2 ) 

C SR  1 3 )=U»YF 
C SR  1 4 ) = M.0 
WRITE  (6,300) 

FORMA  T( 8 X, 'LOCATION* ,8X, • TIME* , 8 X, • MACH  * , 1 0 X , *U2 ' , 1 OX , ' U 1 ' , 1 OX , 
l 'P2'110X,'P1',9X,'R02'  ,9X,'R01'I 

WRITE  (6.301) 

FORMA  T)l4Xt'T2'f10X,»Tl',10X  ,»A2'f10X,'Al',10X,*S2',l0x,'Sl«  , 

1 7 X , 'DS/OT ',4X,'D1A ME TER',4X, 'FRICTION') 

WRITE  (6,109)  (SW<I> *1=1,18) 

WRITE  (6, 109)  Si 
WRITi.  (6, 109)  SR 

FORMA  1( ( IX,  ' SWr' , 2X , 9 ( 2X , E 1 0 .4 ) , 2X , • MA I N • ) / ) 

DO  33  1K=1,I1 
HL  = SW  (1  )-OIS 
HR=SW (1 ) ♦01 S 

HHX  -P  (1  , IK,  1 ) 

MK  = 0 

IF(HHX.GT.HL  .AN0.HHX.LT.HR  ) KK=1 
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IMHHX.Co.HL.OR.HHx.EO.MR)  KKr-1 
IF(KK  » 3 5,33,3** 

3*  P(4,IK,1»=A2 

P 15,1 K, 1 )=S2 
P < 8 , I K*  1 )=P2 
Pi  10,  IK,  1 ) = T2 
60  ?0  33 

35  P <4 , 1 K, 1 >=SPEED*BIC < 4 ) 

BH=CP  "ML  06  ( TEMP  I -RR*At06  (PRES) 

P ( 5 , I H, 1 )=BH/G/RR«BIC(5) 

P<8, I Kt 1 >=PRES*BIC(8  > 
P(9,IK,1)=0ENS*BIC<9> 

P < 1 0 , IK,  1 >=TEMP*BIC  < 1 0 » 

33  CONTINUE 
U-XNN  *U 
C SOLUTION 

IS=ISEG/2 
IW=IS  *1 

00  13  J - 1 , H J 
DO  12  IZ=1,MI 

1 = IW-IZ 

IF(IZ.GE.IW)  I = I ^ 

00  10  K = 1,H 
X1(K)  =P(K,I,J> 

IP2  = 1*2 
X 2 ( K ) =P  ( K , I P 2 , j ) 

10  CONTINUE 

00  30  LI=1,N 
ML  ( L I kSLUI  ) 
gP(LI  l=SR(LI  ) 

30  CONTINUE 

CALL  FIG <M,X 1 ,X2 ,X3 ,WL,WR > 

00  31  LI -1 ,N 
SL(L1  )- WL  ( L I I 
SR  < L I »=WR<LI  > 

31  CONTINUE 

00  11  K--  1 , M 

IP1  = 1*1 
JP1  = J*1 
P(K,I  tJ|=Xl  (K  ) 

PtK,I P2, J>=X2(K> 

P(K,1 FI, JP1 ) z X 3 (K ) 

11  CONTINUE 

12  CONTI NUE 

13  CONTI  NUE 
C END  SOLUTION 
C PRINTOUT 

WRITE  (6,1) 

1 FORMA  1ilHl,6X, ’ARRAY  P(K,I,J>*> 
WRITE  *6»200)  XNN 
200  FORMA  1(2X, »XNN  = ’ ,F10.3) 

00  20  J=1,JJ,1Q 
00  19  1 = 1, ll 

WRITE  (6,2)(I,J,(P(K,I,J),K=1  ,M>) 
2 FORMA  T(2X,I3,2x, 13, 12I2X,E8.  31) 
19  CONTI NUE 
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SUBROUTINE  RAN4NN ,BBl ,BB2  ,PR  ,XX  ) 

COMMON  /ODD/  6 , CP , C V , RF » R ,RR  ,U  , A , C A , CB  , C C ,H X , A Z , ZE R 0 ,F , XL 
DIMENSION  BB1 <NN  ) ,BB2 (NN ) 

XA=1.CM iG*1.0>/(G-l.O»)*PR 
XB=I6  -J.0)/(G-1.0)*PR 
662(  1 >=BB  1 C 1 » 

BB2  <2  »=BB1  < 2 ) 

BB2(6  ) = BB 1(b) 

BB2 ( 7 )=BB1(7> 

BB2  < 8 )-PR*BB  1 (8) 

BB2(9  )=BBU9)*txA/xB) 

BB  2 ( 1 C) - BB 1 ( IO)*PR* ( XB/X A ) 

BB2  < 1 l)=BBl( 11) 

Vl=SORT(CA^CB*PR) 

V2=V1*(BB1<9)/BB2‘9>) 

BB2 ( 3 ) = U»XX*( Vl-92) 

BB2(4  )=SQRT ( G+R96B2 ( 1 0) ) /A 
BH=CP  *ALOG (BB2( 10)/ BBi ( 10  > ) -RR*ALOG(PR ) 

BB2 ( b l:BH/G/RR*BBl(5) 

RETUR  A 
END 


SUBROUTINE  SUB(N,A,B,C) 
COMMON  /DLL'  CSL » CSR  , SL  * SR 

DIMENSION  A <N  ) »B  (N  ) ,C  (N  ) 

NM  = 0 

CALL  CCS|Nf  A,B,C,C:>L  ,NM> 

I E ( NM  ) 109,109,107 

109  CONTINUE 
MM-  0 

CALL  CCS(N,A,B,C,CSR,MM) 

IF ( MM  ) 110tll0,107 

110  CONTINUE 

CALL  SSS(N,A,B,C> 

107  CONTINUE 
End 


SKH0Mtc!M«?eSUR?»E3l.iN. 

P = <C2  <1 )-Cl< 1 ) )/  < C3«l  »-Cl  < 1 ) > 
DO  IQ  1=3, NN 

C24I)  =P*C3<I  )*ll.O>P)*Cim 
10  CONTINUE 

RETURN 

END 


SUBROUTINE  PPP ( NN , A A 1 , A A 3 ) 

COMMON  J 000/  GfCP,CV,l$F,RtRR  ,U  , A , C A , CB  , CC  , Hx  , AZ  , ZER  0,F  , XL 


DIMENSION  AAllNNI.AAXiium 
ri=l  2 .0/<G-l  .01  j*  l I A A3  o* > 

= ( A A3  ( 5 ) - A A 1 ( 5 ) >*G 
A A 3 < 9 )=AA1  (9)»(C1-C2M.0> 

A A 3 ( 1 C ) = A *A  * A A 3(<4)*AA3(9)/G/R 
AA34B  )=AA3(9)*R*AA3(  10) 

RETUR  N 

ENO  1ylQ 


rrr~ 

C 2 = 

A A 3 


AA 1 U > > /( < AA3 U >*AA 1|*>  >/2 


3 
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SUtROtTlN£  CHKUNN,  AA1  ,AA2,  A A3  , AA4  I 

COMM 0 A /ODD/  G,CP,CV,*F,R,RR ,U , A , C A , CB , CC ,HX , AZ ♦ ZERO ,F , XL 
PARAMETER  M=12,  N = 18 

01 MEN  SION  B2«M),B3(H),Cl(M),C3(M),Dl(M)f 02<M  ) 

OIMEN SION  AA1(NN),AA2(NN),AA3(NN),AA4|NN),SW(N> 

DIMENSION  E2(M»,E3<M)»61(M),F2(M> 

10=0 
ID«  = 0 
IDL  = 0 
100  = 0 

( 2)-AAl (2 > >/ < AA2< 1 )-A A1 ( 1 ) > 

BK  = AA  1(2)-$K*AA1  < 1) 

CALL  G£0(NN,AA1,AA2,  *A3,  AA4  ) 

sfv^o^lH'lSr:l«?IS62’s“’c3*,ow, 

100  RA  = AA  2(  3 )*AA2(4) 

SA= 1 . C/RA 

B*=Cj  12)  -SA*C3(  1 ) 

C 1(  1 > =(  BA-BK >/ (SK-S*  ) 

Cl  ( 2 A =SK*Cl  ( 1 ) »BK 
CALL  COM (NN  * A A 1 tCl * A A2  ) 

CALL  SUB(NN,C1,AA2»C3) 

SW<5) =C3(3) 

SW  ( 7 ) =C3  ( 8 ) 

SW  ( 9 ) =C  3 ( 9 ) 

SW( 11  )=C  3(10) 

SW  ( 13  )=C  3 ( 4 ) 
sw<  is  )=C 3(G) 

SH=CP AaLoGCSW( 10)/$W( 11 ) ) -RR*AL0G(SW(6)/SW<7  » » 

SWt 14  ) = SH/G/RR  + SW(15) 

WRITE  16,109)  (SU(I) ,1=1,18) 

109  F0RNAT((lX,,SW=•,2X,9(2x,El0.4)f2X,•109L,)/) 

PR=  SW  (6)/SW<7) 

CALL  RANtNN,C3,D2»PR,XN) 

CALL  SlML«NN,AAl, AA3,*W,B3, I OL ) 

IFCIOL)  104,104,103 
103  RB  = A A 1(  3 ) -A  A 1 < 4 ) 

SB  = 1 . C/RB 
BB  = B 3 (2>-SB*B3(l  ) 

B2(l) :(BB-B«)/ISK-SB) 

B2(  2)  =SK*Cl  ( 1 )*BK 
CaLL  C0M(NN,AA1,B2, AA2) 

CALL  SUB (NN,AA1,B2,B3) 

SW(5)  =83(3) 

SW<7)  =B  3 ( 6 ) 

S W | 9 ) =B  3 I 9 ) 

SW(  11  )=B  3 ( 1 0 > 

SW(  13  1=63(4  ) 

SW( IS  )=B 3(5) 

SH=CP  AAl06(SW(10)/SW( 11 ) )-RR*ALOG(SW (6) /SW<  7>  » 

SW( 14  )=SH/G/RR«SW (15) 

WRITE  (6,110)  (SW (I), 1=1, 18) 

110  FORMA  T(  ( IX, * SW=* , ?X ,9 (2x ,E10 .4 ) ,2X , * 110L ' )/ ) 

PR=  SW  (6)/SW  <7 ) 

CALL  BANINN,B3,01,PR,XN) 

GO  TO  104 

101  CALL  SIML (NN , AA1 , AA3,SW,E3, IOD) 

I/dOC)  104,  104,99 
99  F 2 ( 1 > =SW 111 
F 2 ( 2 ) :SW ( 2 ) 

F2(3)  =SW(5) 

F2(4)=SW(13)  149  . 

J 


L 


F 2(5) =SW (15) 

F21 6) =SW(16) 

F2(7»  =SW(17> 

F2(  10  >=SW4 11  ) 

F2(  11  >=SW(  18  ) 

F2(  12  >=-1.0 

RC  = A A 1(  3 ) - A A 1 ( 4 ) 

&CZ1.C/RC 
B C = F 3 (2)  -SC*E3(  1 ) 

SS=(F2|2)-AA1(2))/(F2(1)-AA1(1)) 

B S = A A 1(2 ) ~S$*AA  111) 

E2<1>  =(BS-BC)/(SC-SS> 

E 2 ( 2 ) = SS*E 2 1 1 > ♦ B S 
CALL  cOM(NN,AA1 ,E2,F2  ) 

CALL  SUB(NN,AAl,E2,E3l 
SU  ( 5 ) =€3(3) 

SWI7)  =£3(8) 

SW<9>  =£3(9) 

S«( 11  )=£  3(10) 

SW(  13  ) = E3(4> 

SW( 15  )=E3(5) 

SH=CP 4A LOG ( SW ( 10 )/SU ( ll > > -RR*ALOGI SW (6) /SW  ‘7  ) ) 
SW(  14  )=SH/G/RR  + SW(15> 

WRITE (6,111 >(SW( I), 1=1,18) 

111  FORMA  T(  ( IX,  'SW=*  ,2X,9  (2x,E10  .4)  ,2X  , * U1L  ’ )✓  ) 
PR=SW46)/SW<7> 

CALL  fiAN(NN,E3,Gl,PR,XN) 

104  1 F ( 1 0 G ,LE.  0 .AND.  100  .LE.  0)  GO  TO  105 
IF ( ID  R .GT.  0 .AnD.  IDL  .LE.  01  G®  TO  lOb 
IF ( 10  S • G T . 0 .AND.  IDL  .GT.  0)  GO  TO  107 
IF ( 10  R .LE.  0 . AND • 100  .GT.  0)  GO  TO  108 

106  CALL  SUBINN, AA1.D2.AA3) 

10= 1 0 t 

GO  TO  105 

107  CALL  SUB(NN,D1,D2,AA3) 

10  = 10  7 

GO  TO  105 

ine  CALL  SUB (NN  * G 1 * AA2 , A A3 ) 

10=106 

105  RETURN 

end 


HX  » A Z .ZERO,) 


PARAMETER  M= 1 2 , N=18 

OIMENSION  B2(M),B3(M),Cl(M),C3(H),0l(M),02(M) 

DIMEN  SION  AA1(NN),AA2(NN),AA3(NN),AA4(NN),SW(N) 
01 MEN  SION  El  (M),£3(M),G2(M),F1 (M) 

10=0 


IDL  = 0 


IOR  = U 
100  = 0 
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.XL 


w 


XN- 1 • C 

SK-(AA2I2)-AA1I2))/(AA2(1)-AA1(1)) 

BK  = AA  1(2)  ~SK*AAl  I 1 ) 

CALL  GEOINN, AA1, AA2, AA3, AA4 ) 

CALL  SIMRINN,AA1,AA3»SW,B3,I0L> 

IFIIDL)  101,101,100 

100  R A=  A A 1I3)-AA1I4) 

SA=1.  C/PA 

6 ArB  3 (2 > -S  A*  B 3 1 1 ) 

B 2 ( 1 ) :( B A -BK )/( SK-SA) 

B 2 ( 2 ) =SK*B2(  1)*BK 
CALL  C0MINN,AA1,B2,AA2) 

CALL  SUB(NN,AA1,B2,B3> 

SWI  51  :B3<  3) 

SU( 7)  -B 3 ( 8 ) 

SWI 9) =B3I9) 

SWIll  l=B  3 1 1 0 > 

SW  I 1 3 )-B3  14  ) 

SWI  IS  )=B3(  5) 

SH=CP  *ALoGI SWI 10  >/SW I 1 1 ) ) -RRaALOGI SW  16 ) /SW  <7  ) ) 
SWI 14  ) = SH/G/RR*SWtl5) 

WRITE  16,109)  ISWII) ,1=1,18) 

109  FORMA  Tl I IX, ’SW:’ ,2X,9(2X,E10.9  j ,2x, *109r  • >/) 

PR=SWI6)/SW<7) 

CALL  *ANINN,B3,D1,PR,XN> 

CALL  SlHRINN.AA2.AA2,sW,C3,IDR) 

IF  I ID  fi)  104,104,103 

103  RB=AA 21 3) ,AA2I 4) 

SB  - 1 • C/RB 

BB  = C3  12  ) -SB  *C  3 I 1 » 

Clt 1)  :l BB -BK ) / I Sk -SB  1 
C 1 1 2 ) =S»<  *C  1 ( 1)4BK 
CALL  C0M(NN,AA1,C1,AA2) 

CALL  SUBINN, Cl, AA2, 03) 

SWI  5)  =C3  13) 

S W 1 7 • -C  3 < 8 ) 

SW|«J  ZCU<)> 

SWI 11 ) = C3I 10) 

SW  I 1 3 )-C  314) 

SWI1S  )=C3I5) 

SH=CP  AALoGi SW I 10 )/SW 1 11 ) ) -RRaALOGISW  16 ) /SW  < 7 » » 
SWI 14  )=SH/G/RR4SW 1 1 S ) 

WRITE  lb, 110)  ISW II ) ,1=1,18) 

110  FORMA  T( I IX, ’SW=*,2X , 9 I 2X ,E10.4) ,2X, ’ 110R*)/) 
PR-SW  (6  ) /SW  I 7 ) 

CALL  RaN(NN,C3,02,PR, XN) 

GO  TO  104 

101  CALL  S1MRINN, AA3, AA2,SW,E3, IOD) 

IF  (IOC)  10ll,104,99 

9 9 Fill)  =SW(1> 

F1I2)  :SW(2) 

F 1 1 3 ) =S  W I 5 » 

FI  14)  :$W  113) 

F 1 1 5 ) :sw(15) 

F 1 1 6 ) :SW(lt>) 

F1I7)  =SW<17) 
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nm-iinm 

F 1 ( 10  ) = SW(  1 1 ) 

F 1<  11  >=  SU  ( 1 8 > 

Fl(12)=-1.0 

RC  = AA 2( 3 ) *AA2 (4 ) 

SC=1. C/RC 

BC  = E3  12 ) -SC*E  3 ( 1 > 

SS=(F  1(?I-AA2(2> )/(Fl tl>-*A2  (1 ) ) 

BB-A A 2(2 ) “SS*AA2 ( 1) 

E l‘l)  = ( B B - BC ) / I S C -S S > 

Ell?)  =SS*E  1 1 1 ) ♦BB 
CALL  cOH< NN, f 1 ,£ 1 ,AA2 > 

CALL  SUB(Nn,E1,AA2,E3> 

SU(S)  =£3(3) 

SWI  7)  =E3(8) 

SW<9>  =£3  (9) 

S«(  11  )=E3(  10) 

SUI 1 3 )=£ 3(4) 

SUI  IS  )=E3(5) 

Sh=CP«aCoG(SU(10)/SU(11> )-RR>*AL0G(SU(6)/SW<7)  > 
SUI 14  )-SH/G/RR*SW (1 5) 

URITt(6»Ul>  ISUII)  ,1  = 1,18) 

111  FORMA  Tl I IX, *SW  = * ,2X ,9 (2x,E10  .4 ) t2X , • 11 1R*  >/> 
PR:  SU  (6  ) / SU  ( 7 ) 

CALL  RAN(NN,E3,G2,PR, XN) 

104  IF(IDl.LE.O.AND.IOO.LE.O)  go  To  105 
IF(IUL.GT.O.AnD.IDR.LE.D)  go  to  106 
1F(10L*GT»0.*N0»10R.GT.0)  60  TO  107 
IF ( ID  L. LE .O.AND*IDD*GT *0 ) GO  TO  108 

106  CALL  SUBINN, 01 fAA2,AA3> 

IDzlQ* 

GO  TO  105 

107  CALL  SUBINN, 01, D2.AA3I 
10=107 

GO  TO  105 

106  CALL  SUBINN, AA1,G?,AA3) 

10=108 

105  RETURN 

End 
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SU6R0  LTI  NE  GE0(NN,AA1,AA2»AA  3.  A A <4  ) 

01 MEN  SION  AA1 INN ) * A A2 INN ) »AA3lNNJ , A A 4 I NN ) 

R13  = (AAl(3MAAi(4>'*AA3(3)*AA3t4))/2.0 
S13-1  .0/R13 

R23=(AA2(3>*AA2<4»*AA3m-AA3<4ll/2.0 

S23- 1 .0/623 

B13  = AAH2I-S13»AA1(1) 

B23  = a A2C2>-S23*AA2( 1 ) 

A A 3 I 1 »=(B23-B13)/(S13-S23) 

AA3 ( 2 )*S1 3*AA3  1 1 >*B13 
S12=(AA1(2I-AA2(2)»/IAAH1»-AA2<1)» 

8 1 2 = A A1  4 2J-S12*AA1< 1 ) 

R34  = ( AA3  ( 3 ) ♦ A A4  ( 3 ) ) /2 .0 

S 34  - 1 .0/R34 

B34=A  A3(2»-S34*AA3( 1) 

AA4 ( 1 >=(B12-B34>/<S34-S12» 

AA4I2  )=S34*AA4< 1 ) ♦B34 

PF=  (A  A2I  1)  - AA4I1  >)  / I AA2(  1 l-AAl  I 1 ) > 

AA4IS  )-( 1.0-PF)*AA2(5)»PF*AA1(5) 

C HEATING 

XK  = 0.  C 

EI=1Q  .0 

E0=15  .0 
A A 3 ( 6 ) = 0.0 

A A 4 I 6 >=0.0 

JF(AA  2( 1 ) .GE.EI. AND. A A3  I 1 >.LE.EO>  AA3lfc  >=XK 
IF(AA4(1).GE.EI.AND.AA4(1).LE.E0)  A A4 ( 6 I - XK 
PH=(A A2( 1I-AA3III I / ( A A2 ( 1 » - A A1 { 1) ) 

A A 3 ( 1 1 ) = 1 1 » O-PH )*AA2(11)*PH*AA1(11) 

AA3I7  )=0 . 25 
A A 4 I 7 >=0.25 


syeRgiiUNE  si  ml  inn.al.ar  .sw  . aa  , io> 

COMMON  /ODD/  G , CP , C V * RF , R ,RR ,U , A , C A , CB , C C ,HX , AZ , ZER 0 , F , XL 
PARAMETER  N=18 

01 HEN  SION  AL(NN),AR<wN),sW(N>fAA(NN) 

10  = 0 

IFCAL  ll> .EQ.SUI 1 ) .AND.SWI2) . EQ. AZ > Go  TO  108 

1F( AR  (1) .EO.SU< 1 ) •AN0*AR(2) . EO.SU < 2)  ) 60  TO  108 

IF(AL  II > .GT.SW( 1 > > GO  TO  108 

00X=SU4 Il-AR I n 

IF<0DX.GT.HX>  GO  TO  iOb 

0IA  = 0 .25 

S=l  AL  12  ) -AR  I 2> ) / I AL  I 1 )-AR  1 1 ) ) 

6=AL<  2 >-S*AL41> 

C initialise 

VAS  = SW ( 3 1 
VA  = SW  ( 4 > 


1 
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V 


££  = Sit) 

LL-0 

101  CONTINUE 
VBS  = VAS 
VB  = VA 
P8  - PA 
LL-LL *1 
LHAX  = 5 

LN  = Ll  -LMAX 
IF ( LN > 103.103.105 
103  CONTINUE 

R2=0.  5*{SW(3)+VBS»*SU<5> 

S2  = 1.0/f?2 

B2  = SU  < 2)-S2*SW  l 1) 

A A ( 1 ) -1  B2-B  ) / ( S-S2 ) 

OIS=Sfc( 1 l“AA ( 1) 

IFtDI S.LT.ZERO)  GO  TO  108 
VAS  = ~ (C  A tCB*  (PB/PP 1 )**0.5 
VP  = CC  *(VAS-1.0/VAS) 

VA  = VP  *SU  <5  » 

C CONTINUITY  FOR  A STANDING  SHOCK 
VizVA S 
V 2=  VI  -VP 
D2  = SK  (91+V1/V2 
X1  = SU  I6)/SW(8) 

X2  = 0.  5*  ( S U ( 4 ) * S V (4) * A * A ) 
X3=0.5*(VA*VA*A*A)*«1.0^D1S*F/DIA) 

PA -02 *<X1>X2-X31 
C TEST  ROUTINE 
TX=0.  C2S 
11:0 

RV=  ( V A-VBl/VA 

RP: ( P A-PBl/PA 

IF  (AB  S<RV  ) .GT.TX  > 1 1 = 1 1 ♦ 1 

I F C AB S(RP> .6T .TX » 1 1 = 1 1 ♦ 1 

IF ( 1 1 ) 1 02.  102,101 

102  CONTINUE 
1 05  CONTINUE 

IF<  AA  11)  .GE.ALt  1 l.ANO.AA  <1  ).  LT.ARU  I 1 10=1 
I F ( I 0 » 1 08,1 08,107 

107  A A ( 2 ) =S*  A A < 1 > ♦ B 
c reoefine  SHOCK 

SW(  1 » =AA  m 
SW  ( 2 > =AA(2» 
sum  =V  AS 
sum  =v  a 

S U < 6 I :PA 
S U { 8 1 =02 
SW( 10  )=PA/<0 2*R) 

SU( 12  )=SORT<G*R*SW( 10) >/A 
SU(  lb  >=0.0 
SVC  17  )=0.2S 
SU( 18  >=0.2 

l£»9  FORMA  T(  ( IX, ‘SW:* ,2X,9  <2X  »E10  .4  ) ,2X  , ' SI  ML  * I / ) 
IFjLn  ) 108,108,104 
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10%  cO*TI 

WRITE  (6,106)  TX  ,RV, RP ,SW ( 1 > , SW ( 2 > ,LMAX 
106  F ORMA T(lX,5(2X,E10.4»,3Xf •LMAX=»,I5,2X,*SlML»» 
106  CONTINUE 
RETUR  N 
END 


D I MEN  SION  AA1(NN)»AA2(NN),AA3(NN),AA4(NN),HH3(NN) 

00  10  1=1, NN 

HH  3 ( 1 UAASt1  ) 

10  CONTI  NUE 
S4  =A  A 4(5  ) 

0 S4  3 = UA4 (6 ) ♦ KH3 ( 6 ) >/2.0 
T 3=  HH  2(2) 

T 2=  A A 2(  2 ) 

Pl  = (2  .0/ (G-l .0) >*AA1 (4)*AA1 ( 3) 

0 2= <2  .0/(G-l ,0) )*AA2(4)-AA2(  3) 

A13  = ( AA1 (4MHH3(4  ) 1/2.0 

A 2 3 = ( -AA2  ( 4 ) ♦ H H 3 ( 4 ) ) / 2 • 0 
U 1 3= ( AA 1 ( 3 ) ♦HH3 ( 3) ) /2 . 0 
U23= ( AA2(3)*HH3( 3J1/2.0 
D3  = HH  3(7) 

D 1 = A A 1(7) 

02  = A A 2(  7 ) 

E 3=HH  31 1 ) 

El  = AA  k 1 ) 

mmn  ( 6 ) ♦HH3 ( 6 ) 1/2.0 

QS2  3 = UA2(6)*HH3(6)  )/2.0 
F13=( AA1 ( 11 ) *HM3 ( 11 ) )/2.0 
F23=(  AA 2(11) ♦ HH 3 ( 11  ) >/2.0 

0 1 3=  ( AAl(7MHH3(7))/2»0 
023= ( AA2 ( 7 > ♦HH3 (7>)/2.0 
T3  = HH3(2) 

T 1 = AA  1(2) 

T2=AA  2(2) 

T4=  A A 4(  2 ) 

S 1 = * * 1(5) 

S2=  AA  2(5) 

S 3 = S 4 40  S43*  ( T 3-T4  ) 

P 3=P 1 ♦(  (G-l .Q)*A13*DS13 

1 -(  F 1 3*XL  / (2.0*01 3) >*U13*ABS (U13 ) )*(T3-T1) *A13* (S3-S1 ) 

03=02 ♦ ( (G-l .0)*a23*0S23_, 

1 -( F2 3 *XL/ (2.0*023 ) ) *U 2 3* a BS ( U 2 3 ) ) * ( T 3-T 2 ) ♦ *2 3* ( S 3-52 > 

A A 3 ( 3 )=(  P3-03>/2.0 
A A3 (4 )=(G-1.0)*(P3*03)/4.0 
A A3  ( 5 )=S  3 

20  FORMA T( 2X, 3( 3X.E20. 12  ) ) 

RETURN 

ENO 
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?VMRk?£BEN^li<NNN,AAA1,AA*2  »AAA3,SWL,SMR> 

OIMEN SION  AAA1INNN)  , AAA2 (NNN  >, A A A3 (NNN  > 
PARAMETER  NX  = 15 

DIMENSION  AA1CNX ) .AA2<Nx ) ,AA3(NX ) , AA4CN*) 
DIMENSION  SWL(N>  ,SWRCN> 


NN-NN 

DO  1U 


N 

1 = 1, NN 


A A 1 ( I )=A  AA I ( I ) 

A A 2 ( I ) = A A A 2 ( I 1 
10  CONTI  KJE 

IFC  AA  1(2  ) .Lt.  (O.OJ)  AA1(2)=AA2(2) 
IFIAA 2I21.LT. I 0.0))  A A2 ( 2 ) = A A1 (2) 
A A 3 ( 3 )=  < A A 1 ( 3 ) ♦ A A 2 ( 3)  )/2.0 


A A 3 ( 4 >=(AA1(4)+AA2(4)  ) / 2 • 0 


A A 4 ( 3 ) = A A 3 { 3 ) 

nm=o 

CALL  ChKR(NN,AA1,AA2,AA3,AA4  ,SwR,NMJ 
IFiNM  » 109,109,107 

109  CONTINUE 
MM  = 0 

CALL  CHKL(NN,AAl,AA2,AA3,AA4  ,SML  , MM ) 
IF (MM  ) 110,110,107 

110  CONTINUE 

CALL  SUBINN, AA1.AA2  .AA3) 

107  CONTINUE 


M.lhHhfh. 

A A A 2 ( 1>=AA2( I > 

A A A 3 ( I)  = AA3( I ) 

11  CONTINUE 

WRITE  16,2)  A AA3 
2 FORMA  T‘1X»12<2X,E8.3>) 
RETUR  N 
END 


Wl^A,CB,CC,HX,AZ,ZEROtF,XL 


PARAMETER  N = 18 

DIMENSION  ALCNN>,AR<NN)»SW<N),AA(NN) 

10  = 0 

IF  C ABO)  .EQ.SWC1  ).AN0.SWC2>.EQ.AZ)  60  TO  108 
I^(AL  (l).E0.SUil).AN0.AL(2).E0.SU(2) ) 60  TO  108 
IF(AR  m.LT.SWUI)  DO  TO  108 
QDJUAU  l)-SM(l) 

IF(OD X.GT.HX ) GO  TO  108 
0JA  = 0 .25 

S=  ( AL  (2)-AR(2) ) / ( AL  ( 1 )-AR  « 1 ) ) 

B=AL ( 2>-S*AL(l) 

c initialise 

WAS  = SW(3) 

VA  = SMC  4) 
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Jl 


i 


I 


P A T SW  ( 6 ) 

pp  = sw  (7 1 

LL  = 0 

101  CONTINUE 
VPS  = VAS 
VB  r VA 
PB  = FA 
LL  = LL  ♦! 

L HA  X i 5 
LN-LL  -LMAX 
IF  I LN > 103*103,105 
l03  CONTI  KUE 

R2=0.5*4SWI3)*VBS)*SWlS) 

S 2 = 1.0/r2 

62  = SWI  2>-S2*SW 4 1 > 

A A 4 1 ) -< B2-B )/ IS-S2) 

0 I S - A A ( 1 > -SW ( 1 ) 

IFIDI S.LT.ZERO)  GO  TO  108 
VAS  = 4CA*CB*4PB/PP ) >**0.5 
VP  = CC  •( VAS-1 • 0/ v A S ) 

VA-VP «SW  <5 ) 

C CONTINUITY  for  a standing  shock 
V 1- VA s 
V 2:  V 1 -VP 
02=SW  49  ) *V  1 / V 2 
xl  = SW  16  ) /SW  4 ® 1 
X2  = 0.  5*4  SWI  4 ) *SW  t 4 ) *A*A> 
X3=0.5*IVA*VA*A*A)*I1  .0-*DIS*F/DI  A 1 
PA:D2*(X1*X2-X3) 

C TEST  ROUT  lNE 
TX=0. C25 
11=0 

R V = ( V A*V6>/VA 
RP  = ( P A-P6 ) /PA 
IF  I AB  SIRV  1 *GT  *TX  1 1 1 = 1 1 ♦ 1 
IFCAB SIRP) .GT.TX ) 1 1 = 1 1 ♦ 1 
IF ( 1 1 > l02,  102, 1Q1 

102  CONTINUE 
10S  CONTINUE 

IF  I A A ID .GT . AL ( 1 » .AND. A A ( 1 ) . LE . AR  I 1 1 > 10  = 1 
IF  do  ) I°8,l08,l07 
107  AA  4 2)  =S*  A A I 1 ) *B 
C redefine  shock 
SW4  1 1 =A A ( 1 ) 

SWI2)  =A  A 4 2 ) 

SW«  3)  =VAS 
SW(  4)  :VA 
SW(  6)  =P A 
SW48)  =02 

SW  4 10  )=PA/ID2*R> 

SWl 12  ) = SQRT  I G *R * S W I 10)  )/A 
SWI  lb  >=0.0 
SWl 17  >=0.2S 
SWI  lb  1=0 . 2 

l09  F0RHATlllx,'SW=»,2X,9,2X*El0.4)»2Xt*SIMR,l/) 
IF  4 LN  ) 108,108,109 
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^W8kT5B!)D5C^C^!|ev;^^^ft1,ll!A,CA,CB,CC,HXtA2,ZERO,F 

PARAMETER  NZ  = 4 

parameter  nx-is 

DIMENSION  X (M  > ,Y  (M)  ,Z  <M> 

DIMENSION  FACE (NZ > , Ab (NX > , YY  (NX  ) 

I N-  0 
XHX  -4  .0 

IF«VI  D.LT.FACEU  ))  60  TO  200 
XXL:X  (1  ) -FACE  <1  > 

IF< XX L.6 T.xHX ) GO  TO  200 

IF ( X ( 1) .EO.F ACE ( 1 ) . AND.X ( 2 ) . EO.F ACE ( 2 ) ) GO  TO  200 
IF<y<  1).EQ.FACE(1I.AnO.Y(2).EQ.FACE(2)I  GO  To  200 
Z ( 3) = (XI 3 > ♦ Y ( 31 ) /2.0 
Z ( H ) - (X(4)*Y(41)/2»0 
Ab(  31  :Z(  31 


XL 


RIGEtLtEF8A^♦x’Y’^’8B, 

S2  3 = ( 2t 2 >-Y ( 2) )/ (Z( 1 >-Yl 1 ) ) 

B23  = Y (2  ) ~S  2 3*  Y ( 1 1 
RCS=F  ACE  ( 3) 

SCS=1  .O/RCS 

BCS-F ACE  4 2>-SCS*FACE( 1 ) 
E=(BC5-B23)/(S23-SCS) 

T = S23*E*B23 

IF (E . G£ .Z ( 1 > . AND .E.LE • Y ( 1 I)  GO  TO  101 
GO  TO  102 

101  CONTINUE 
IN=1D  1 

pL=(2.0/(G-1.0n*X(4)*X(3> 

0R=<2  .0/(G-1.0))*Y(4)-Y(3) 

SL-X  ( 5) 

SR  - Y ( -) 

CAL=0.25*(G-1.0)*(SL-SR) 

OL=(P  L*QR )*TANH (CAL ) ♦QR 
PRr ( 0 R-QL ) ♦PL 
V:(PR  -0U/2.D 
C = 0.2  £*<G-1.0»*(PR*QL ) 

Y Y ( 1 ) zt 

Y Y ( 2 ) =T 
YY  ( 3 ) =V 

Y v ( 4 ) =C 

Y Y ( 5 I :X(5  1 

Y Y ( 6 ) =0.0 
YY( 7) =0.25 

Y Y ( 11  ) = 0 . 2 

CALL  SSS(M,X,YY,Z) 

F*Cf  C 1)  zt 
FACE  ( 2)  = T 
FACE  < 2)  =V 
WRITE  (fe. 1 ) FACE 

1 F0RMAT(lX,4(2XtE10.4>,2X,,INTERFACE,> 

102  CONTI  *JE 
200  continue 

RETUR  N 
END 
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104  6rIt£(^106>  tx,rv,rp,sw(1) ,sw(2) ,lmax 

106  FOR  HA  HlX»6(2XfE10»‘4)»3X»,LHAX-,,I5*2X,*SlHR,l 
108  CONTINUE 
RETUR  A 
ENO 


ayR!^iiKEA^^HY;^l^f;AAAM|NN, 

DIMENSION  HH3 ( 15 > ,AA4 (15 > 

AA4(3  )=  A A 3 ( 3 ) 

A A 3 ( 3 1=  ( A A 1 ( 3 > ♦ A A ? < 3)  1/2*0 
A A 3 ( 4 )=  ( AA1(4)*AA2(4)  ) / 2 • 0 
LL-0 

101  CONTINUE 
LL=LL  *1 
I.MA  X = £ 

LN=LL  -LMAX 

I F ( LN > 103*103.105 

1Q3  continue 

CALL  GE0(NN,AA1,AA2 »AA3*AA4  ) 

CALL  «0L(NN,AA1, AA2.AA3, AA4,HH3  » 

C TEST  ROUTINE 

TX  = 0. C01 
11  = 0 

RV=< A A3( I ) -HH  3(3) )/AA3(3> 

R A= ( A A3( 4 1-HW3 (4 ) )/AA3 (4  » 

I F ( A b S(RV1  .GT.TX  > I I =1 1 ♦ 1 

!F(AbS(RA)  .GT.TXj  1 1 ; 1 1 ♦ l 

1F( II  ) 102, 102,101 

102  CONTI AUE 
GO  TO  104 

105  WRITE  (6*106)  TX , R V , R A , AA 3 ( 1 ) , A A 3 ( 2 ) , LM AX 

106  FORMA  T(  1*,5(2X,E10.4  ) »3X,  *LMAX  = *,I5.2X,  * S SS  ’ ) 


to- 

(NN,  AAl , A A 3 ) 


A A 3 ( 1 2)  =LL 
RETUR  A 
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APPENDIX  D 


d 


201  SERIES  MINI-GAGE 


Couplers 

To  complete  the  system  a wide  choice  ot  couplers  is  ottered  in  the  648  and  549 
senes.  Optional  filters  and  several  types  of  power  inputs  are  available  Model  587D 
Coupler  provides  increased  capability.  Please  refer  to  our  Coupler  Data  Sheet  tor 
details 

In  addition  Models  583  5030  and  5040  Laboratory  Amplifiers  are  offered  with 
extensive  versatility  and  many  options  tor  more  complete  pressure  studies 


t 

Specification* 

Model  Variation 

PERFORMANCE 

UNITS 

201B1 

201 B2 

'0164 

201B5 

t 

Pressure  Range,  5V  out 

psr 

5.000 

500 

200 

100 

( 

Overrange  - 

PSI 

7,500 

750 

300 

150 

Reaolution  (noise) 

psi  rms 

0 05 

0 005 

0 002 

0001 

1 

Maximum  Preaeure 

psi 

15.000 

5 000 

2,000 

1,000 

Sensitivity  

mV/psi 

1 

10 

25 

' 50 

Linearity,  B.F.S.L. 

o 

-1 

• 1 

!*  1 

2:1 

i 

Resonant  Frequency,  nom. 

k Hr 

500 

500 

500 

250 

Rise  Time,  10-90%  

n sec 

1 

1 

1 

2 

Time  Constant,  R T. 

sec 

1.500 

400 

200 

100 

♦ 

Low  Frequency  Response,  —5% 

Hr 

0 0003 

0 001 

0 0025 

0005 

High  Frequency  Response,  -t-5% 

Hr 

100,000 

100.000 

100,000 

50,000 

» 

ENVIRONMENTAL 

Common  Specs 

if 

Vibration  Sensitivity,  max. 

psi/g 

0 002 

0.002 

0 002 

0.002 

Shock,  1ms 

g 

5,000 

5.000 

5.000 

5.000 

Vibration  Limit 

q 

500 

500 

500 

500 

t 

Temperature  Range  

F 

-65  to  280 

-65  to  280 

-65  to  280 

-65  to  280 

Temperature  Sensitivity  Shift 

/ F 

0 03% 

0.03% 

003% 

0 03% 

I 

ELECTRICAL 

: 

Output  Current,  min. 

mA 

2 

2 

2 

2 

Polarity,  pressure  Increase 

Negative 

Negative 

Negative 

Negative 

Bias  Voltage 

V 

11  2:2 

11  ±2 

11  ±2 

11  2:2 

f 

Circuit  Return 

Case 

Case 

Case 

Case 

Output  Impedance,  max. 

ohms 

100 

100 

100 

100 

MECHANICAL 

r 

Weight  

gms 

< 10 

< 10 

< 10 

< 10  j 

Case  and  Diaphragm  Material  .. 

Stainless  St. 

Stainless  St. 

Stainless  St. 

Stainless  St. 

Mounting  Torque 

m-lb 

24 

24 

24 

24 

f 

Sealing 

All  Welded 

All  Welded 

All  Welded 

All  Welded 

POWER  SUPPLY 

* 

! 

t 

Constant  Current  Source 

mA 

4 2:1 

4 ±1 

4 2:1 

4 2:1 

* 

Supply  Ripple,  max.  

mV  rms 

25 

25 

25 

25 

Supply  Voltage,  no  load 

VDC 

20-30 

20-30 

20-30 

20-30 

Source  Impedance,  nom.  

ohms 

250  k 

250  k 

250  k 

250  k \ 

i 


